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Pure Yang-Mills instantons are considered on x R'^ - so-called calorons. 
The holonomy - or Polyakov loop around the thermal at spatial infinity - is 
assumed to be a non-centre element of the gauge group SU{n) as most appro- 
priate for QCD applications in the confined phase. It is shown that a charge k 
caloron can be seen as a collection of nk massive magnetic monopoles - com- 
ing in n types - each carrying fractional topological charge. This interpretation 
offers a physically appealing way of introducing monopole degrees of freedom 
into pure gluodynamics: as constituents of finite temperature instantons. 

Using the Nahm transform an elaborate treatment is given for arbitrary topo- 
logical charge and new exact and explicit solutions are found for SU{2) and 
charge 2. The k zero-modes of the Dirac operator in the background of a charge k 
caloron are computed, and are shown to 'hop' between the n types of monopoles 
as a function of the temporal boundary condition. The abelian limit - where it 
is assumed that the massive field components can be dropped - is analysed in 
great detail and is shown that the abelian charge distribution of each monopole 
type coincides with the corresponding fermion zero-mode density. 

The 4n/c dimensional hyperkahler moduli space is identified as an algebraic 
variety, defined by four matrices obeying a constraint, modulo a natural adjoint 
action. This moduli space is proved to be the same as the moduli space of stable 
holomorphic bundles over the complex projective plane which are trivial on 
two complex lines. A description is given for its twistor space which allows for 
the computation of the exact hyperkahler metric - at least in principle. 

Finally, lattice gauge theoretic applications are mentioned and is explicitly 
demonstrated how to obtain calorons on the lattice using the method of cooling. 
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Chapter 1 
Introduction 



More than 30 years after its formulation quantum chromodynamics is still not solved, yet 
there is overwhelming evidence for its correctness. One of the most important phenomena 
of QCD is quark confinement. It is poorly understood in terms of first principles and yet 
this phenomenon is vital for our understanding of basic properties of hadronic matter and 
its interactions. The primary reason for the lack of understanding for this non-perturbative 
phenomenon is the fact that we do not know how to describe the true vacuum of QCD, i.e. 
we do not know which are the "good" degrees of freedom to study the dynamics in the 
strongly coupled infrared regime. 

On the other hand perturbative calculations work reliably for short distances due to 
asymptotic freedom. In fact non-abelian gauge theories are the only known examples of 
asymptotically free quantum field theories in four dimensions. In this regime the free 
fields serve as "good" degrees of freedom and quantum fluctuations around the unique 
perturbative vacuum are under precise calculational control. 

Quark confinement is a phenomenon that is present in QCD for small enough temper- 
atures. As the temperature is increased a phase transition occurs at a critical value which 
separates the confined and deconfined phases. Above the critical temperature the quarks 
are liberated - deconfined - and together with the gluons form a quark-gluon plasma. The 
actual mechanism that causes the quarks to confine below the phase transition is believed 
to be the same as the mechanism at zero temperature. Even though it is not directly obvi- 
ous how the limit of zero temperature affects the microscopic description, robust phenom- 
ena such as confinement ought to survive the limit. For this very reason, if the presence of 
finite temperature in QCD made it easier to address some non-perturbative effects present 
at zero temperature as well, there is no reason not to formulate the theory at non- vanishing 
temperature. 
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1.1 Gauge theories 

Not only QCD but essentially all the fundamental interactions, as we know them at present, 
are successfully described by gauge theories. Despite their simple formulation there are 
several inherent puzzling features. The basic field of a gauge theory is not observable and 
does not have any physical meaning. If one wants to correct for this and use variables 
which have clear physical meaning and in particular are gauge invariant then the theory 
quickly becomes utterly complicated. Being forced to use gauge dependent variables leads 
to a whole zoo of possible computational schemes each corresponding to different gauges. 
However, as gauge invariance is important and any physical answer should refer to only 
gauge invariant quantities, additional computations are neccessary to check the gauge in- 
dependence of the results. 

The seemingly innocent looking Lagrangian - presented below - hides the complexi- 
ties underlying non-abelian gauge theories. The innocent look is partly due to the fact that 
if some reasonable assumptions, such as gauge invariance, locality and Lorentz invariance 
are imposed on a theory describing spin 1 fields then the Lagrangian is essentially unique 
and natural. 

For gauge group SU(n) the dynamical variables are the 4 components of an n x n anti- 
hermitian matrix valued field A^(x), the gauge potential. In terms of the field strength 
jf>v = 9^ Ay — dvAfi + [Afi, Ay] the Lagrangian of pure Yang-Mills theory is 



^=--^TrF2„ (1.1) 



where ^ym is the dimensionless coupling constant and for simplicity we assume a metric 
of Euclidean signature. Note that the minus sign is included in order to make the action 
non-negative. The equations of motion that follow from the Lagrangian are 

Dj,F^y = 0, (1.2) 

where = 3^ + A^ acts in the adjoint representation. Here we are considering a theory 
without matter, which would otherwise give a source term to the right hand side. 

The action, or equivalently the equations of motion, are invariant under gauge transfor- 
mations. A gauge group valued field g{x) acts on the gauge potential and correspondingly 
on the field strength as 

F^v gF^vg-^. (1.3) 

Another invariance of classical Yang-Mills theory in 4 dimensions is conformal sym- 
metry. If the metric is rescaled by an arbitrary spacetime dependent factor, the action does 
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not change. This is because the Lagrangian in the presence of an arbitrary - but still of 
Euclidean signature - metric g^v is 

^ = -;rVv^rt""Trf,,f,,, (1.4) 

2^YM 

and if the metric is rescaled by any local factor A(x) as g^iv — ^ ^gfiv then the inverses change 
according to g^'^ ^^^g^^> which is exactly cancelled by the change in the volume factor 
^/detg A^^detg. 

The above analysis was classical and conformal symmetry is destroyed by quantum 
fluctuations, it is even broken on the perturbative level. Physics is not the same on all 
scales. The coupling constant changes with the energy scale with which the system is 
probed, more specifically with the momentum transfer involved 

dgYM af \ /I c\ 

H^^ = ^(^ym), (1-5) 

where }i is the renormalization scale. For non-abelian gauge theories the /3-f unction is 
negative at small coupling, resulting in a decrease of the coupling constant at large ener- 
gies. This phenomenon is called asymptotic freedom and is a direct consequence of the 
self -interaction of the gluons, that is of the non-abelian nature of the theory IIHEI. 

The Lagrangian of non-abelian gauge theories we have presented is not the full La- 
grangian of QCD, only of its bosonic sector. The Dirac fermion fields j/? are in the fun- 
damental representation of the gauge group and come in Uf flavours. Each flavour / 
transforms as i/^y^ gtpf under gauge transformations and it is easy to see that the full 
Lagrangian 

^ = -^Tr + J^tp^(]l)- mf)tpf (1.6) 
/ 

is also gauge invariant. The parameters m ^ are giving bare masses to each flavour and ^ is 
the hermitian covariant Dirac operator. As long as the number of flavours is small enough 
the anti-screening of charge due to the self -interaction of the gluons is over compensating 
the usual screening present also in abelian theories and the /3-function remains negative 
for small coupling. 

In the massless limit a new symmetry emerges - at least classically. The infinitesimal 
transformation by an x Uj anti-hermitian matrix oj, 

Sxpf = C0ff,y5Xpf, (1.7) 

leaves the action (ll.6l t invariant if m^r = 0. This chiral U{njr) symmetry is, however, broken 
as a result of the quantum dynamics. More precisely, the axial U{1) of U{nf) is broken by 
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instantons through an anomaly as we will see in the next section and the remaining SU{njr) 
group is broken spontaneously. The order parameter of the phase transition associated to 
the spontaneous breaking is the chiral condensate {ipip), where an averaging over flavour is 
implicit. Even though classically and to every flnite order of perturbation theory it remains 
zero, in the full quantum theory {xpxp) 7^ 0. A formula due to Banks and Casher relates the 
chiral condensate to the spectral density |0( A) of the Dirac operator around zero eigenvalue, 

(#) = 7rlim lim (1.8) 

if the theory is formulated in finite volume V and with non-vanishing masses m for each 
flavour ||3l. The order of the two limits is important, first the thermodynamic limit should 
be taken, followed by the chiral limit. The spectral density p{A)dX counts the average 
number of eigenvalues of Ip between A and A + dA and thus the Banks-Casher formula 
relates the chiral condensate to the low-lying spectrum of the Dirac operator. We will see 
in the next section that instantons dramatically affect the spectrum of the Dirac operator, 
in particular they give rise to zero-modes and hence are of significant importance for the 
phenomenology of chiral symmetry breaking. 

Due to the running of the coupling constant a dimensionful parameter, Aqcd/ has to 
emerge in the theory. This new parameter is essentially the constant of integration that 
naturally appears when solving H1.5|l and fixing Aqqd fully specifies the theory with no ad- 
justable parameters. In particular the coupling constant will also be fixed by the j6-function 
equation ()1.5|l . This fashion of trading a dimensionless coupling for a dimensionful one is 
called dimensional transmutation. 

As we will be concerned with certain classical solutions of Yang-Mills theory, dimen- 
sional transmutation does not play a role and we will put gyu = 1- 



1.2 Instantons 

Topological excitations are special gauge configurations in Yang-Mills theory ||4l|5|. They 
are required to have finite action and be stable minima of the action functional. As a result 
they are solutions of the equations of motion. However, the requirement of stability puts 
further constraints on them besides eq. ()1.2|l and a quick inspection of the following trick 
provides us with such a constraint 

S = J d^xTrFf^y = -\ j ^"^^Tr {F^y±P^yf'±^ j d^xTrF^yP^y = 

= -\j d^^Tr {F^y±Pf,yf±8n\, (1.9) 
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where f^^y = ^e^vpaFpa stands for the dual field strength and we have introduced the topo- 
logical charge 

' d^xTrFuyPuy. (1.10) 



1671^ 



For the prototypical example of spacetime being it is an integer once the action is re- 
quired to be finite. In this case the field strength must go to zero at infinity and hence the 
gauge field must be a pure gauge A^ = U^^dfiU, where U is only defined on the bound- 
ary S^. Such SU{n) mappings are classified up to homotopy by an integer which is 
exactly given by ()1.10|l . 

It follows from the above trick that 87r^|fc| < S, and equality is achieved if and only if 

F^y = ±F^,y, (1.11) 

which are the celebrated (anti)self-duality equations depending on the ± sign. They are 
really three equations, 

foi = ±^23 / Fq2 = ifsi / f03 = ±^12 • (1.12) 

If a configuration is (anti)self-dual then it automatically satisfies the equations of motion, 
however the converse is in general not true. Self-dual configurations are called instantons 
if > and anti-tnstantons if < 0. For reviews, see [6n7n8t|9J. 

There is an alternative definition in terms of chiral fermions that is useful. Using the 
representation 

— Z(7j 

for the Dirac 7-matrices the covariant Dirac operator becomes 



7„ = ( " ) . (1.13) 



P = i7,D, = ( ) . (1-14) 

where we have introduced the chiral and anti-chiral Dirac operators (also called Weyl op- 
erators) D = (Tj^iDn and = — o-^D^;. Here the cr^ are the basic quaternions; our notation 
is summarized at the beginning of chapter |2| In this representation the chirality operator 
is 

75 = ( -J ° ) , (1.15) 

where the blocks are 2 x 2. It is easy to check that 

D+D = -D;,D^ - h^vF^,v , (1.16) 
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with the familiar anti-self-dual 't Hooft tensor fj^y = fJliv^i- Since the contraction of a self- 
dual and an anti-self-dual tensor vanishes, we have the following alternative definition: a 
gauge field is self-dual if and only if the corresponding D^D operator is a real quaternion 
and in particular commutes with the quaternions. In this case it equals the negative of the 
covariant Laplacian. For anti-self-dual fields the definition is similar with the role of D and 
interchanged. 

Whether or not a gauge field is (anti)self-dual, its topological charge is always given by 
the formula ()1.10|l and is always an integer as long as the field strength falls off faster than 
for large x. A non-vanishing topological charge has dramatic effect on the spectrum 
of the Dirac operator and will be described below. 

Since the Dirac operator anti-commutes with the chirality operator, {l/),^^} = 0, it fol- 
lows that its real non-zero eigenvalues come in pairs. If A is an eigenvalue with eigenmode 
xp then 'y^ip is an eigenmode with eigenvalue —A. Also we see that if is a zero-mode then 
so is j5Xp and then the combinations 5(1 ± 75)1/' also zero-modes and are eigenmodes 
of 75 as well with eigenvalue ±1. Thus in the space of normalizable zero-modes the basis 
vectors can be chosen with definite chirality. Denote hy n± the number of normalizable 
zero-modes with chirality ±. It follows from the explicit forms I)!. 14111. 15)l that in terms 
of the chiral and anti-chiral Dirac operators, n+(n_) is the number of normalizable zero- 
modes of D(D+). 

We now wish to demonstrate the sum rule fc = n — n+. To this end we consider the 
quantum field theory of a massive fermion coupled to a classical gauge field 1 6 1. The action 
is. 



t>-m)xp, (1.17) 



where the gauge field in P is treated classically, thus only ip{x) and ip{x.) are integrated 
over to compute expectation values. The simplest chiral Ward identity in this theory states 
that 

Mili"")) = -minsH^)) - ^Tr FF(x) , (1.18) 

where = 1/^7^(751/' is the current associated to the axial U{1) transformation, see ()1.7|l . 
The first term on the right hand side is due to explicit chiral symmetry breaking by non- 
zero mass and the second term is the famous Adler-Bell-Jackiw anomaly |10, llj present 
even in the massless limit. Integrating the Ward identity over all of spacetime gives 



m 



{J <fxxp^s^) = -k, (1.19) 



because the left hand side in dl.lSI) is a total derivative and since we are dealing with a 
massive theory no contribution can come from the boundary. The vacuum expectation 
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value - or more precisely the expectation value in the presence of a classical gauge field - 
on the left hand side can be computed using the fact that 

(t/;(x)t/i(y)) = -^^(x-y) (1.20) 

is the known propagator, which gives immediately 

( J d^xtpys^P) = -Tr (75^^) . (1.21) 

Now we have seen that for non-zero eigenvalues the eigenmodes xp and 751/7 belong to 
different eigenvalues, hence are orthogonal. As a result the evaluation of the trace is con- 
veniently done in the basis of eigenmodes and only the zero-modes contribute. Due to 
75 in the trace those with chirality (+) contribute 1, those with chirality ( — ) contribute — 1 
leading to 



d'^xip^s^) = — , (1.22) 

which together with H1.19|l is the desired result k = n — n+, also called the Atiyah-Singer 
index theorem L12, 131. It holds for any gauge field, whether or not it is a solution. We now 
specialize to instantons. 

We have seen that for instantons the covariant Laplacian factorizes, 

D+D = -D^,D^ (1.23) 

and that the number of normalizable zero- modes of D is n_i_ Suppose that n_|_ > 
0. In this case there is a normalizable zero-mode xp for D, Dxp = 0. Applying to both 
sides, multiplying by xp^ and then integrating over spacetime gives 

= J d^xip^D^Dip = J d^xip\-D^D^,)ip = j d^xD^^ip^D^ip, (1.24) 

which is only possible if xp is covariantly constant, contradicting its normalizability. Hence 
n+ = and the index theorem for instantons states that D has no normalizable zero-modes 
whereas has as many as the topological charge of the underlying gauge field. This result 
will be heavily used. 

The Banks-Casher formula ()1.8|l relates the chiral condensate to the low-lying spectrum 
of the Dirac operator hence it is not surprising that instantons play a crucial role in the 
dynamics of chiral symmetry breaking. 

The nature of the fermionic zero-modes can be illustrated by the plots of the exact 
solutions. The field strength square of a generic charge k instanton looks like k lumps. The 
k zero-mode densities xp^xp are such that - after choosing an appropriate basis - they peak 
roughly at the location of the lumps. This harmony between the fermionic and bosonic 
degrees of freedom is expressed by saying that each basic charge 1 instanton carries its 
own zero-mode. 
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1.3 Finite temperature 

In the previous section spacetime was IR^ and had Euclidean signature. This is appropri- 
ate for a Wick rotated Minkowski spacetime or for a finite temperature system in the limit 
of zero temperature. For truly finite temperature one has to consider the imaginary time 
being periodic with period l/k^T where fcg is Boltzmann's constant and T is the tempera- 
ture. Hence the manifold will be x R'^ over which the (anti)self-duality equations will 
be studied and (anti)self-dual configurations will be called calorons. For a comprehen- 
sive discussion of instantons in finite temperature QCD, see IITSll . for a review of caloron 
solutions, see flS]. 

The motivation is the desire to understand or at least be able to say something about 
the phase transition of QCD. The order parameter is the vacuum expectation value (p(x)) 
of the trace of the Polyakov loop 

l/keT 

p(x) = iTrPexp J Ao{t,x)dt. (1.25) 


For high enough temperatures - in the deconfined phase - {p{x)) is close to one of the 
n^^ roots of unity, each representing a vacuum. Clearly, the choice of any specific one out of 
the n possibilities breaks the Z„ symmetry associated to cyclically permuting the n vacua. 
If p{x) is close to a root of unity then its length is fluctuating around its maximal value 
\p{x)\ ~ 1, which is only possible if the eigenvalues of p{x) are all close to being the same. 

On the other hand for low temperatures (p(x)) = and the Z„-symmetry is restored. 
The fact that p{x) is fluctuating around zero means that the length of the Polyakov loop 
is minimal. A simple exercise reveals that this is only possible if the eigenvalues are close 
to being as different as possible. Let us denote the eigenvalues by exp (2m'/i^), ordered as 
^\ < ^2 < ■ ■ ■ < ^n, then 

IpWI^ = 4 E cos27r(^^ -}Ib) + -- (1-26) 

" A<B ^ 

Clearly, two coinciding eigenvalues give a large contribution to |p(x)|^ as cos27r(^y^ — ^g) 
attains its maximum for = ^g. The fluctuations in the Z„-symmetric - or confined - 
phase are such that the eigenvalues repel each other. 

This is our motivation for studying topological excitations - calorons - in a Polyakov 
loop background with no coincident eigenvalues. Since the path ordered exponential 
around a closed loop is also called a holonomy, such a Polyakov loop is also referred to 
as having non-trivial holonomy. 

The mechanism of confinement at zero temperature should be the same as at finite 
temperature T as long as T < Tq. In addition in the real world definitely T > and in 
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most circumstances T < Tc, except maybe at RHIC or in the early universe. Hence under- 
standing confinement for finite temperature is perhaps not enough to collect $1,000,000, 
but is sufficient to understand confinement in the real world [16J. Undoubtedly, studying 
topologically non-trivial solutions of classical Yang-Mills theory at finite temperature will 
not solve or in any way explain this non-perturbative phenomenon as for large coupling 
or low temperature semi-classical arguments are insufficient. Our motivation is solely to 
reveal the true degrees of freedom in the topologically non-trivial sector of QCD at non- 
zero temperature. We find that in the confined phase instantons dissociate into magnetic 
monopoles changing the character of the basic topological object present in the QCD vac- 
uum. 

Just as instantons play a crucial non-perturbative role at zero temperature, we believe 
that a similar role is played by the constituent monopoles that take the place of instantons 
at finite temperature and especially in the confined phase. The true quantum dynamics 
of these monopoles or the quantum dynamics of any degree of freedom for that matter is 
beyond our considerations as we are working in the semi-classical regime, but we would 
like to stress that isolating the "good" variables is the first step in formulating a dynamical 
model. 

Having non-trivial holonomy is essential for arriving at massive constituent monopoles 
and the arguments presented above are in favour of such a scenario. One note, however, 
is in order when discussing the dynamical importance of configurations with non-trivial 
holonomy. It was observed that the one-loop correction to the action of configurations 
with a non-trivial asymptotic value of the Polyakov loop gives rise to an infinite action 
barrier and hence these configurations were considered irrelevant ITU . However, the in- 
finity simply arises due to the integration over the finite energy density induced by the per- 
turbative fluctuations in the background of a non- trivial Polyakov loop Ill7l . The proper 
setting would therefore rather be to calculate the non-perturbative contribution of calorons 
- with a given asymptotic value of the Polyakov loop - to this energy density, as was first 
successfully implemented in supersymmetric theories [18J, where the perturbative contri- 
bution vanishes. The resulting effective potential has a minimum where the trace of the 
Polyakov loop vanishes, i.e. at maximal non-trivial holonomy. 

In a recent study at high temperatures, where one presumably can trust the semi- 
classical approximation, the non-perturbative contribution of the monopole constituents 
was computed IIT9ll . More precisely, the effective potential due to the one-loop determi- 
nant in a caloron background was computed. When added to the perturbative contribu- 
tion with its minima at center elements, a local minimum develops where the trace of the 
Polyakov loop vanishes, deepening further for decreasing temperature. This gives sup- 
port for a phase in which the center symmetry, broken in the high temperature phase, is 
restored and provides an indication that the monopole constituents might be the relevant 
degrees of freedom in the confined phase. 

Monopole based models in the spirit of a dual superconductor that are conjectured to 
lead to confinement were introduced long ago 11201 12l1l . Traditionally, monopoles are static 
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objects in non-abelian Higgs models, but such a Higgs field is absent in QCD. Another 
possibility is abelian projection [22, 23 J but in this approach magnetic monopoles enter es- 
sentially as gauge singularities and their physical interpretation - i.e. gauge independence 
- is not so clear. The alternative we offer to introduce monopoles into QCD, through the 
constituents of finite temperature instantons, is gauge invariant and physically appealing. 

1.4 Collective coordinates 

The set of collective coordinates that may enter the most general instanton or caloron solu- 
tion carries important information about their physical interpretation. Some of the param- 
eters are interpreted as gauge orientations or phases, some others as scales and locations 
Il24il . Exploring the whole moduli space is necessary to identify the role of every parameter. 

What we find is that the 4nA: dimensional moduli space of instantons which consists 
of (4n — 5)k gauge orientations, k 4-dimensional locations and k scales is traded at finite 
temperature for nk 3-dimensional locations and nk phases. The interpreation of these pa- 
rameters is clear, they describe nk magnetic monopoles. A more detailed analysis of the 
moduli space confirms that this is not just an arbitrary juggling with the possible ways of 
factoring the number 4nfc, but really a physically sensible reorganization of the collective 
coordinates takes place. In particular a charge k object is not an approximate superposition 
of k charge 1 basic objects, but rather an approximate superposition of nk objects each with 
fractional topological charge. 

There is an apparent puzzle that seems unavoidable following our discussion of the chi- 
ral fermion zero-modes in the instanton background. We have seen that a generic charge k 
instanton can be thought of as an approximate superposition of k charge 1 instantons each 
carrying a fermion zero-mode. If at finite temperature we have nk basic objects in a charge 
k configuration then how can chiral zero-modes be supported on all of them if the index 
theorem still dictates that only k zero-modes exist? The answer is that the nk monopoles 
come in n distinct types, each corresponding to a U{1) subgroup. In addition, the pres- 
ence of finite temperature necessitates a choice of boundary condition for the zero-modes 
in the compact time direction, they can be chosen to be periodic up to an arbitrary phase, 
exp(— 2m'z). According to the Callias index theorem [25J, for a given choice of this phase, 
say < z < }iA+i> the k zero-modes localize to the k monopoles of type A only. When- 
ever exp(27rzz) passes an eigenvalue of the Polyakov loop, the zero-modes hop from one 
type of monopole to the next, eventually visiting all of them. For z = fij^ the zero-modes 
delocalize or spread over both types A — 1 and A. 

It should also be noted that for finite temperature field theory there is a canonical 
choice, namely that the fermions are anti-periodic. Thus from a physical point of view 
n — 1 out of the n possible boundary conditions are non-physical and perhaps this was 
part of the reason why constituent monopoles were not seen in lattice gauge theoretical 
studies in the past. If one, however, performs simulations with the non-physical boundary 
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conditions as well, the behaviour predicted by the exact solutions is revealed. In this sense 
the zero-modes are used as probes of the underlying gauge configuration rather than as 
dynamical, physical fermions. 

The same comment applies to supersymmetric gauge theory compactified on x R^. 
There is also a canonical choice of boundary condition in this case, the (adjoint or fun- 
damental) fermions should be periodic in order to preserve supersymmetry. This is the 
right choice for computing the caloron contribution to the glutno condensate for instance 
in ^ = 1 super Yang-Mills theory [18J. The non-physical fermions are nevertheless still 
there and can be used for diagnostic purposes but do not play a role dynamically. 

1.5 Outline 

In the following chapter we will present the well-known results on self-dual Yang-Mills 
fields over flat spaces TP x R'? for p + cj < 4. The organizing principle is Nahm's duality 
on the 4-torus that maps U{n) instantons of topological charge k on to U{k) tnstantons 
of charge n on the dual torus T"* with periods inverted. By sending some of the peri- 
ods to infinity or shrinking them to zero this approach puts the ADHM construction, BPS 
monopoles, vortices, calorons, etc. into a common framework. This will be discussed in 
section IZ2l 

Chapter |3l deals with the application of Nahm's transform for calorons with arbitrary 
topological charge. Our rather general results are specialized to various limiting cases in 
section 13341 that include BPS monopoles and the abelian limit when the non-abelian cores of 
the massive monopoles are shrunk to zero size. Both the general formulae and the limiting 
behaviour are exlicitly spelled out for SU{2) and charge 2 in section |331 

The fermionic sector, concretely the zero-modes of the Dirac operator in the caloron 
background, is investigated in chapter lU It is shown how the zero-modes can be used to 
probe the monopole content of the caloron field by varying their boundary condition in 
the compact temperature direction. The abelian limit of the previous chapter is employed 
to achieve maximal localization for the zero-modes. Again, our general results are made 
explicit for charge 2. 

We also show that upon large separation between the constituents the zero-modes 
"see" point-like monopoles. The exact results from the previous chapter are used to re- 
solve the singularity structure of the abelian limit and we obtain the exact zero-modes as 
well. 

Chapter |5l explores the moduli space of multi-calorons. The essential tool is Euclidean 
twistor theory as applied to hyperkahler geometry. We derive an explicit parametrization 
in terms of finite dimensional matrices. Using this a correspondence is established be- 
tween the caloron moduli space and the moduli space of stable holomorphic bundles over 
the projective plane which are trivial on two projective lines. We also construct the cor- 
responding twistor space which encodes the hyperkahler metric of the moduli. We show 
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1.5. Outline 



that upon large separation the moduli space becomes nk copies of S x IR"^ each describing 
a charge 1 BPS monopole. This observation lends support for our constituent monopole 
picture for arbitrary rank and charge from a geometrical point of view. 

Lattice gauge theory is the natural framework to study non-perturbative phenomena 
in QCD and should be decisive on dynamical questions such as the dynamical importance 
of our caloron solutions. Lattice aspects of our analytical work is presented in chapter |6l 
Monte-Carlo simulations are performed to demonstrate the confinement - deconfinement 
phase transition for SU{2) and exploratory investigations are done in the confined phase 
in search of calorons. 

Finally, in chapter |7| we end with a number of concluding remarks. 
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Chapter 2 

Self-dual Yang-Mills fields 



This chapter will outline the general structure of the self-duality equation F^,i/ = F^^v and its 
various limits. Starting with Nahm's tranformation on the 4-torus we show how to obtain 
a whole web of interrelated systems by shrinking some of the periods to zero or by sending 
them to infinity The most general form of Nahm's duality considered here will relate self- 
duality on X Tl and K.^^P^'? x fl for p + q < 4: where T'? is the dual torus with periods 
inverted. In particular this includes the algebraic ADHM construction of instantons on R^, 
BPS monopoles on and their relation to Nahm's equation, calorons on x R^ and their 
relation to Nahm's equation with periodic boundary conditions, vortices on R^, etc. For 
more detailed geometrical aspects see chapter |5l 

The gauge group is limited to be SU{n) (or U{n)), generalization to other classical 
groups are possible. In fact some of the computations will be done in the Sp{n) series 
forSLr(2) = Sp(l). 

Our conventions are that for the quaternions we use (c7"o,(ry) = (1, — ZTy) as well as 
(oq, (Ty) = (1, iTj) where Ty are the usual Pauli matrices, for 't Hooft's self-dual and anti-self- 
dual tensors we define 

V}iv = vlvO-j = \ i^F^v - O-vO-^i) 

Vliv = v'fiv^j = 2 i^ii^v - ^v^ii) ■ (2.1) 

We have the identities 

O-^O-y + (TyO-^ = 25 
O'^O-y + O-yCr^ = 25^v 

cr^o-fy - ayafji = 2f]^^y . (2.2) 

For a quaternion Re/Im will always mean the quaternionic real /imaginary part, that is for 
q = qjid^i we have ^eq = qo and Iraq = qiai. If a quantity, say A^, has 4 components then 



13 



2.1. Nahm duality on the 4- torus 



without the ^ index we will always mean the corresponding quaternion A = A^Vyi. Spatial 
vectors %[, e[, . . . with 3 indices will be bold, x,e, . . . and their dot product will simply be 
written xu, y^, etc. 

The various indices will be such that a, /3, . . . have values 1, 2 and are used for chiral 
spinors, a,b,c, . . . are the dual gauge indices and are running from 1 to k, and A,B,C,. . . 
are the indices for SU{n) and are running from 1 to n. 

2.1 Nahm duality on the 4-torus 

When formulated on the Nahm transform |26U27J assigns to every generic U{n) instan- 
ton of topological charge k, another instanton with topological charge n and gauge group 
U{k) but living on the dual 4-torus whose periods are inverted II28II . A nice feature of 
this duality is the constructive nature of it. Once the U{n) instanton is given with charge 
k, there is a recipe to construct the corresponding U{k) instanton of charge n although the 
actual computations can be cumbersome. Another property is that applying it twice gives 
back the original instanton, in other words the Nahm transform squares to one. In addi- 
tion, since the moduli space of instantons on carries a natural hyperkahler structure, 
one can show that the transformation is a hyperkahler isometry |28|. 

Considering various limits of the periods of one ends up with correspondences be- 
tween objects in a variety of dimensions and it is hoped that the magical properties of the 
Nahm transform survive these limits. Maybe it is worth a note that some of the analytical 
properties have not been rigorously proved for all cases, but from a physical point of view 
the principle is clear. The particular case of the caloron has actually been dealt with in 
a mathematically sound way in 1291 1301 and henceforth we will not bother with rigorous 
proofs. 

Let us start with a U{n) instanton gauge field Aj,i{x) of charge k on the 4-torus with 
4 periods 2nL^. One can modify the gauge field in such a way that self-duality is not 
violated by adding a flat factor, (x) — > (x) — 2niZn where the Zj, are numbers. Indeed, 
such a shift does not affect the curvature. It is possible to change to + / for any 
integers by applying a periodic U{1) gauge transformation, so it is best to think of the 
Zf, variables as parametrizing the dual torus with periods = 1/L^,. 

Now consider the chiral and anti-chiral Dirac operators in the fundamental represen- 
tation Dz = [Dpi — 2niZii) and = — (D^ — Iniz^) with = 9^ + A^. Generically 
Dz will have no normalizable zero-modes, whereas D\ will have k of them according to the 
index theorem. ^ Denote by ^pz{^) the 2n x k matrix of linearly independent orthonormal 

^Naturally, there is an obvious symmetry between Dz and once the gauge field is changed from self- 
dual to anti-self-dual. 
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2.1. Nahm duality on the 4- torus 



zero-modes of D^, 

D,>z = , j xplxp^ = 1 , (2.3) 

with ak X k identity matrix on the right hand side. The ^z{x) parametrically depend on 
the Zji variables and it is possible to define 

A^{z) = J d^xtpl^^tp, (2.4) 

as a X gauge field on the dual torus which will be referred to as the dual gauge field. 
Gauge transformations arise because there is a z-dependent U{k) choice in the xpz matrix 
of zero-modes. The transformation j/7z(x) ipz{x)g^^{z) for unitary g(z) induces 

^ S^i^S-'-§-^g-'- (2-5) 

It is easy to see through integration by parts that is anti-hermitian and it is in fact 
also self-dual. In order to see this notice that 

D+Dz = -(D,, - 2niz^){D^ - Iniz^) - -fj^yF^cv = -{D^ - 2niz^){D^ - Iniz^,) , (2.6) 

because the original field strength is self-dual and drops out when contracted with the 
anti-self-dual 't Hooft tensor. This argument is the same as our fermionic characterisation 
of instantons in H1.16|l . Thus Dz^Dz is a real quaternionic operator and one can introduce 
its Green function or inverse as an n x n matrix /z(x, y) by 

DtDzfzix,y)=5{x-y), (2.7) 

where S{x ~ y) is the periodic Dirac delta on T"^. In terms of the Green function the field 
strength of A ^ is L28J 

F^y{z)=8n^ j d^xd^yiljl{x)fzix,y)r]^yipz{y) , (2.8) 

which is clearly self-dual. 

It can be shown that the topological charge of A^ is n and also that if A^(x) is gauge 
transformed then A^ (z) does not change. This means that the map ^ A^ is a map be- 
tween gauge equivalence classes of U{n) instantons of charge k on and U{k) instantons 
of charge n on T"^, which happens to be a hyperkahler isometry. Also it holds that applying 
it twice gives back the original instanton, hence the Nahm transformation is an involution. 
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2.2. Decompactification and dimensional reduction 

2.2 Decompactification and dimensional reduction 

The Nahm transformation on the 4-torus serves as a rich source for a whole web of in- 
terrelated models. One way of obtaining them is decompactifying some of the periods or 
shrinking them to zero. In the current section we will review self-duality over flat spaces 
that can be obtained as such limits of T"*. 

If a period tends to infinity then the dual period shrinks to zero and the dual torus is 
dimensionally reduced. In addition the dual field strength will not be exactly self-dual, its 
anti-self-dual part will not equal zero but to some source term coming from the boundary 
which will be non-trivial in the presence of non-compact directions. More specifically, 
the formula H2.8|l will have an additional contribution coming from integration over the 
boundary which is absent for the compact case and this extra contribution will violate 
self-duality. One can show that these source terms are singularities as a function of the dual 
variables z^, hence self-duality will hold almost everywhere. These singularities amount 
to special boundary conditions for the dual gauge field. 

On the other hand if some of the periods are reduced to zero, that is the original torus 
is dimensionally reduced then the dual torus will develop non-compact directions. 

From the above it is clear that the Nahm transform relates self-duality on IRP x to 
self-duality on R'* P '? x for p + q < 4:, hence the dimensionality of the problem goes 
from p + q to4: — p, which in our case of the caloron - as well as in some other applications 
- is a considerable simplification. Instead of solving a four dimensional problem directly 
we can achieve the same by solving a problem in one dimension. 

If all four periods are sent to infinity the dual torus is reduced to a single point, and 
self-duality with source terms over this point will give the ADHM equations. In this way 
it is possible to derive the whole ADHM construction from Nahm duality and this point of 
view may help clarify the mysterious fact that self-duality on is solved by an algebraic 
construction [31 , 32J. 

If three periods are sent to infinity and one is reduced to zero, then the original setup 
becomes the BPS monopole problem on 11331 l34l . The dual description is then on R, 
that is Nahm's equation with specific boundary conditions. This is the original Nahm 
construction of magnetic monopoles Ii26k27I. 

If two periods are decompactified and two are reduced to zero, then we obtain an in- 
teresting case where the duality is between the same type of objects, both descriptions 
correspond to vortices. 

If we only send some of the periods to infinity but keep the remaining finite, then we 
obtain flat 4-dimensional spaces. Calorons on x R"^ will be related to Nahm's equa- 
tion on the dual circle with periodic boundary conditions with some singularities, doubly- 
periodic instantons on R^ x will have a description in terms of vortex equations on 
and instantons in a finite box, that is on R x will correspond to singular monopoles on 
f3 f35",|36l|3ZlPl. 

It is amusing to note that the manifolds R^, R x and are mapped topologically 
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2.2. Decompactification and dimensional reduction 



to themselves and if the periods are chosen the self-dual value L = 1 / \J7ji then even the 
metrics stay the same. 

A more detailed presentation of the several variants of Nahm's transform is presented 
in the rest of this section. 

2.2.1 ADHM construction 

Atiyah, Drinfeld, Hitchin and Manin have given a complete recipe to construct all self-dual 
gauge fields on R"^ with gauge group SU{n) and arbitrary topological charge k II31U32I . For 
comprehensive reviews see [39. 40J. Even though this was the first construction of its kind 
we will interpret it in the light of Nahm's duality, which appeared later. 

In our notation we will follow the literature and construct the instanton gauge field 
A^{x) from some auxiliary data. In our exposition of the Nahm transform on - again 
following the literature - we have constructed an auxiliary instanton out of the physical 
Api{x). In this sense the ADHM construction is the analog of the inverse Nahm transform. 
This is the reason for shifting the physical gauge field A^(x) by the auxiliary variables 
for Nahm's transform and - as we will see - shifting the auxiliary gauge field Bfi by the 
physical variables for the ADHM construction. We hope this remark will make it easier 
to relate the two constructions and clarify the logic behind the notation. 

One starts with four k x k hermitian matrices combined into a matrix of quaternions 
B — B^a^ and an n x matrix of 2-component spinors assembled into an n x 2fc matrix A. 
The analog of the Dz and D\ operators is the (n -(- 2k) x 2k matrix 

A(x) = (g^^) (2.9) 

and its adjoint. We see the appearance of A in A(x) due to the boundary which is absent 
for the Nahm transform on T^. For R"^ the A-dependent terms will play the role of a source 
term as already alluded to in the previous section. 

An instanton solution corresponds to the matrices (B, A) if A^(x)A(x) is a real quater- 
nion, compare with 1)2. 6|l . This condition is independent of x and is in fact equivalent to 
B^B + A^A being a real quaternion, 

Im(B+B + AU) =0. (2.10) 

This gives 3 quadratic matrix equations for the algebraic data (B, A), 

[Bo,Bi]-[B2,B3] = -.{x\X2 + ^\^i) 

[Bo,B2]-[B3,Bi] = -(^A+Ai-A+A2) (2.11) 
[Bo,B3]-[Bi,B2] = ^ (a1 Ai - A+A2) . 
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An additional constraint is that B — x = [B^ — x^) should only be degenerate for k 
points, which is an open condition so will generically hold. Once such a data is given, 
the gauge field, field strength and a number of other quantities of physical interest can be 
reconstructed explicitly. 

To this end let us look for the kernel of A^(x) which generically will be n dimensional. 
Choosing n normalized basis vectors gives an {n + 2k) x n matrix v{x) of zero-modes for 
which 

A^{x)vix) = 0, v\x)v{x) = 1 . (2.12) 

The zero-mode matrix v{x) is the analog of ipz as defined by H2.3|l . The gauge field corre- 
sponding to the data (B, A) can be written as 

A^{x) = v\x)d^v{x) . (2.13) 

It is worth pointing out that the above formula defines a gauge field for any (B, A) matrices, 
however it will only be self-dual if (B, A) satisfies the quadratic ADHM equations ()2.10|l . 
Once (B, A) does satisfy ()2.10|l it is clear that the transformation 

B gBg-' , A Ag-i (2.14) 

for g G \l{ki) leads to a new set of data satisfying (12.101 1. The form of the gauge field shows 
that such a transformation does not change A^. Thus it is appropriate to associate a gauge 
field to equivalence classes of ADHM data where equivalence is understood with respect 
to the transformations H2.14|l . We conclude that the moduli space of instantons can be 
explicitly parametrized by matrices (B, A) satisfying (I2.10t modulo the equivalence ()2.14|l . 

Obviously if y(x) is a solution to H2.12|l then so is v{x)g{x)^^ as long as g(x) is unitary 
and this will induce a gauge transformation on the gauge field ()2.13|l . 

It is possible to solve for v{x) explicitly in terms of (B, A). Substituting directly into 
(IZT2)l shows that 



v{x)=\ ,\ w(x)= B^-xM A^ (2.15) 



is a normalized solution, where the square root of the positive n x n matrix, 

^(x) = 1 + u'{x)u{x) (2.16) 

is well defined and u{x) is Ik x n. In terms of these new variables the gauge potential 
becomes 

A^ = j^tg^^ ^-1/2 ^ ^1/2 a^<^-l/2 . (2.17) 
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The analog of the Green function is the inverse of the real quaternion A^(x)A(x) which 
is an ordinary k x k matrix yielding the definition 

= (A+(x)A(x))-i. (2.18) 

Note that on the Green function fz (x, y) is a bona fide Green function for the second 
order differential operator D\Dz, whereas in the ADHM construction it is the ordinary 
inverse of the matrix A^(x) A(x). We will still call it the Green function and in terms of this 
hermitian kxk matrix we have ll39ll4Tll42ll 

(p-^ = 1-A/^A+ 
TrF^ = -DDlogdet/x (2.19) 



where ^{x) is the 2n x k matrix of k normalized fundamental zero-modes of the chiral 
Dirac operator in the instanton background, aj,iDpiXp = 0, whose existence is guaranteed by 
the index theorem, □ is the four dimensional Laplacian, e = (72 is the charge conjugation 
matrix and we have also introduced the n x n matrices 

(Pj{x)=\crjf,\\ (2.20) 

We see that cp and cpj completely determine the instanton gauge field. 

It is a useful excercise to check the value of the total action and normalization of the 
fermion zero-modes. Both the action and zero-mode densities are given as the four dimen- 
sional Laplacian of an expression, hence the integral over 4-space can be evaluated from 
the asymptotics. The definition (12.181 1 yields for the Green function fx = l/x^ + --- which 
indeed leads to S = 8n^k and / d'^xtp^tp = 1 as it should. 

Formulae ()2.19|l show that in order to perform actual calculations the Green function is 
a useful instrument. Its analog for the caloron will be used extensively for finding the new 
exact multi-caloron solutions. 

The above construction is valid for unitary gauge groups. Other classical groups such 
as Sp{n) or 0(n) can be encorporated by considering their embeddings in higher dimen- 
sional unitary groups [40 J. These embeddings are of the form that the generators should 
preserve some additional structure, along with the hermitian metric preserved by the uni- 
tary group. Invoking the ADHM construction and approriately imposing these conditions 
one can arrive at explicit instanton solutions for any compact semisimple Lie group. Since 
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SU{2) = Sp(l) one can apply in this case two forms of the ADHM construction and in 
practice we will find it convenient to use the Sp(l) realization. 

In this variant of the ADHM construction the matrices are taken to be real, sym- 
metric. The initially n x k matrix of chiral spinors, or equivalently the n x 2k matrix A is 
taken to be a quatemionic fc-component row vector with real coefficients. A" = Aj^cr^,. The 
symmetries of the ADHM data are modified accordingly, one has the same transforma- 
tions as in ()2.14|l but only g E 0{k) are allowed. The main advantage of using the Sp(l) 
construction instead of the SU{2) is that some of the formulae simplify considerably. The 
quantity (p in this case is proportional to the identity matrix, hence a scalar function of x. 
This makes it possible to simplify the formula for the gauge field to 

A^ = ^(pf]^J^ydy(Pj. (2.21) 

To summarize, the initial problem of finding solutions to the self-duality equations, 
which are partial differential equations in four variables for the gauge field, is turned into 
an algebraic problem of finding roots in a system of quadratic equations and finding the 
eigenvectors of a matrix corresponding to zero eigenvalue. In this sense the 4 dimensional 
problem is reduced to a zero dimensional one. 



2.2.2 BPS monopoles 

Assuming that the gauge field is invariant under translations in one of the directions of 
one arrives at the Bogomolny equation for magnetic monopoles. Indeed, if ^ = Aq is 
introduced as a Higgs field and assuming that neither (p nor A, for i = 1, 2, 3 depend on xq 
then the self-duality equation will reduce to 

B, = -D,<p, (2.22) 

where D, = 9, + A, acts in the adjoint representation and B/ = ^SijkFjk is the magnetic field. 
The Xo"independent gauge transformations descend to the 3 dimensional gauge symmetry 
of (Izj^ . 

A- gAg-'-d,gg-\ (2.23) 

As BPS monopoles are closely related to our central object - the caloron - we will summa- 
rize some of the well-known facts; for a review and more details see 1431 1441 . All of these 
facts will be rederived in subsequent sections from the caloron point of view. 

The finiteness of the 3-dimensional action implies that at infinity the Higgs field should 
tend to a constant. The appropriate boundary condition at infinity is then 

(pAB{r) =iSAB(}iA-l^ + - ■ ■) , J^}Ia = J^Ia = 0, (2.24) 
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or any gauge transform of the above, where the I a are integers. The magnetic charge 
of the monopole is then {k\,k2, . . with k/^ = h- The numbers ifi^ are the 

eigenvalues of the Higgs field at infinity. 

It turns out that the fields {(p,Ai) are not the most convenient objects to describe a 
monopole. One can define so-called spectral data instead, which are in a one-to-one corre- 
spondence with gauge equivalence classes of monopoles B31II5I . For gauge group SU{2) 
the spectral data consists of a spectral curve which for our present introductory purposes 
is simply a complex polynomial of order k = ki = —k2 > with leading coefficient 
1 and another polynomial q{f]) of order k — 1 such that they have no common root. These 
two polynomials can be combined into a rational function, r{f]) = t^i^) /pi^), and the re- 
markable fact is that the rational function r completely determines the monopole solution 
up to gauge transformations. If it can be written in the form 



r{rj) 



k 

L 

a=\ 



(2.25) 



for non-zero complex numbers Ca = exp(u„ -|- ita) and arbitrary Za then such an r represents 
an approximate superposition of k charge 1 monopoles with U{1) phases exp(z^a) and 
locations {—\va, Zg) G R x C = R^, provided the separation between the locations is large 
enough and the choice fii = —^2 = 1 is made L43J . 

If the polynomials have expansions p{r]) = J/'^ + Efl=o Ps^" and q{r]) = Yj^zlqa^]'^ 
then the complex coefficients po, p\,..., pk-i, qo> qi, ■ ■ ■ , qk-i parametrize the 4fc real dimen- 
sional moduli space. The requirement that the two polynomials should not have a common 
root is expressible as A(p, q) ^ 0, where 



A{p,q) 
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(2.26) 



is a {2k — 1) X (2k — 1) determinant, the so-called resultant of the two polynomials. This 
description of the moduli space of monopoles as the complement of the algebraic variety 
^(P/^) = iri C^*^ is useful as it provides an explicit holomorphic parametrization and 
with some extra work the hyperkahler metric can also be derived II43I . 

If the rank of the gauge group is larger than one then one has a similar picture with 
spectral data for each SU{2) embedding with constraints on the polynomials. We will 
revisit these questions in chapter|5lin more detail and will also show how the spectral data 
of monopoles arises from our construction of caloron solutions and their moduli spaces. 
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2.2.3 Vortices 

If the gauge field on is assumed to depend only on xq and Xi, equations relevant for the 
study of doubly periodic instantons and vortices are obtained 1*461 1471148114911 . In this case 
there are 2 Higgs fields, (pi = Ai, (p2 = A3 and a 2 dimensional gauge field Aq i both of 
which can be combined into complexified fields (p — (pi + i(p2 and A = Aq — iA\. Also, it is 
convenient to work in complex coordinates z = xq + ^^i- The self-duality equations reduce 
to 

B = -{<p,^\ (2.27) 
D(/) = (2.28) 

where B = [D, D] = 9 A - 9 A + [A, A] is the curvature of A and D = 9 + A is in the 
adjoint representation. The symmetry becomes ^ S^S^^ ^^"^ ^ S -^S^^ ^ ^S^S^^ 
for SU{n) valued gauge transformations but note that eq. H2.28|l is invariant under the 
complexified gauge group SL{n, C) whereas H2.27|l only under the original compact group. 
This feature is a general phenomenon also occuring in the other dimensionally reduced 
examples but perhaps is most transparent in the present case. It will be explained and 
exploited in chapter |5l in the general context of hyperkahler geometry. 

2.2.4 Nahm equation 

The most important case for our purposes - for calorons - is the dimensional reduction to 
1 dimension. Assuming that the gauge field only depends onxQ = t self-duality becomes 

A/+[Ao,A,] = le,-^fc[A^-,Afc], (2.29) 

the celebrated Nahm equation, where prime denotes differentiation with respect to t. It is 
an ordinary non-linear differential equation and also plays an important role in the study 
of rotating rigid bodies. Gauge transformations only depend on t and act as ^ 

Ao ^ g^'^-g'g'^ 

^ (2.30) 

The detailed analysis of Nahm's equation will be done in the next chapter where we use it 
to construct new multiply charged caloron solutions. 

■^An interesting observation is that if the range of t is compact and periodic boundary conditions are 
imposed - as for the caloron - then the transformation of Aq is the same as the coadjoint action of the 
centrally extended loop group familiar from WZNW models |50|. 
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2.2.5 Reduction to zero dimension 

There is still a fourth possibility, namely to reduce to zero dimensions and assume that the 
gauge field does not depend on any of the coordinates. In this case only the commutator 
terms in the field strength survive and we obtain 



recognizing immediately the ADHM equations as introduced in section 12.2.11 More pre- 
cisely, the ADHM equations are the above equations in the precence of a source given by 
the A-dependent terms. The conclusion is then clear; the ADHM construction - or from 
our point of view the Nahm transform - relates the four dimensional self-duality equa- 
tion and its moduli space to self-duality in zero dimensions, hence supplying an algebraic 
solution to the former. 

2.3 Existence and obstruction 

Nahm's duality transformation states that if a U{n) instanton of charge k exists on T^, so 
does a U{k) instanton of charge n, if we identify with the toplogically identical T"*. It 
follows then immediately that there can not exist a charge one instanton on the 4-torus 
L28 I. If it existed, the Nahm transform would produce a U{1) instanton of charge n, which 
is clearly impossible as LZ(1) gauge theory is linear. 

One can show that no such obstruction exists for higher charge on |51 1. It is strongly 
believed that on T'^ x R unit charge instantons exist, although it is not proved. For the 
remaining cases, x R^, x R"^ and R'* it is known that instanton solutions exist with 
any topological charge. Those for x R"^ will be discussed in the next chapter, with an 
emphasis on higher topological charge. 



[Ao,Ai] 
[Ao,A2] 
[Ao,A3] 



[A2,A3] 

[A3,Ai] 
[Ai,A2], 



(2.31) 
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Chapter 3 



Multi-caloron solutions 



Nahm duality - see section IZTI - tells us that in order to construct multi-caloron solutions 
one should study Nahm's equation on the dual circle [52J. This method transforms a solu- 
tion of a non-linear 4 dimensional partial differential equation to a solution of an ordinary 
but still non-linear equation. The precise boundary conditions for the dual gauge field, 
formulae for physically interesting quantities and other details of the construction can be 
obtained in a number of ways. One could start from the Nahm transform on T"^, then care- 
fully perform the limit of 3 periods tending to infinity and trace what terms arise as sources 
that violate self-duality, see the comments after eq. H2.8|l . Another possibility is first let 
tend to R** thereby ending up with the ADHM setup and then compactify one direction a 
la Fourier, resulting in x R^. Yet another option is to start from BPS monopoles on R^ 
with corresponding dual discription on R and compactify this R in order to have x R^ 
in the original setup [53 J. The compactification will introduce the time dependence that is 
absent in the pure monopole situation. These approaches are equivalent and we will use 
mainly the second. 

Calorons interpolate between instantons on R** and BPS monopoles on R"^ by varying 
the radius of the circle corresponding to finite temperature. Because of this it is not sur- 
prising that calorons share features with both objects and for the actual construction one 
can use a mixture of the ADHM and BPS monopole methods. 

We will be seeking solutions of multiple topological charge for which the asymptotic 
Polyakov loop, or holonomy, defined as the path ordered exponential at spatial infinity. 



is a generic element with all eigenvalues exp(z/3^^) distinct. In addition, we assume an 
ordering jii < }i2 < ■ ■ ■ < jin- This requirement means that the SU{n) symmetry is 
maximally broken to 11(1)"^ by the holonomy or equivalently that all the monopoles 




(3.1) 
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inside the caloron will have non-vanishing masses, Stt v^^/ /3, where = ^a+i ~ j^A- The 
massless limit giving rise to so-called non-abelian clouds |„54 , 55 J can be taken in a more or 
less straightforward way but we will not be concerned with it here. 

The solutions we obtain will have one special feature though, they will have vanish- 
ing over-all magnetic charge. Including an arbitrary magnetic charge {ki,k2, . . . ,kn-\), as 
mentioned in section l2.2.21 would mean that the rank of the dual gauge field as defined on 
different invervals is not a constant but jumps according to the differences fc^+i — 156J. 
The boundary conditions for these cases are also known but for the sake of simplicity we 
will limit ourselves to vanishing over-all magnetic charge. We will see that the monopole 
constituents come in n distinguished types each being associated with a pair of adjacent 
eigenvalues (f/^, f/^+i). Each of these types has a magnetic charge in the corresponding 
1/(1) subgroup, but the total magnetic charge of the sum of all constituents will be zero 
however. 

Lattice gauge theory considerations also justify the interest in only zero over-all mag- 
netic charge as the simulations are performed in a finite box. Clearly, in finite volume with 
periodic boundary conditions there can be no net magnetic charge. 

Without loss of generality we set the radius of the circle to unity, i.e. /3 = 2/1 which 
results in the period of the dual circle to be 1. 

3.1 Dual description of calorons 

Nahm duality tells us that we should consider the chiral and anti-chiral Dirac operators 
on the dual circle parametrized by z in the background of a ll{k) dual gauge field A, 

D = ^ + cr^A^,{z), = -^-a^A^,{z) (3.2) 

and the requirement of self-duality for A is equivalent to D^D being a real quaternion. 
We have also seen that since x R"^ is not compact, the dual gauge field is only self-dual 
up to singularities and the precise form of the singularities can be obtained by Fourier 
transforming the ADHM equations |57, 58 1. The source term in the ADHM equations was 
Im A^A. Adding the Fourier transform of this term to the self-duality equations for A yields 
the dual description of Sll{n) calorons in terms of the ll{k) dual gauge field A{z) and an 
n xk matrix of 2-component spinors or equivalently a2n x k matrix A, 

A/ + [Ao, Ai\ - -£,yfc [A^-, Afc] = z J] ^(z - ^^)pf , pfd^ = i Im A+P^A , (3.3) 

^ A 

where the prime stands for the derivative with respect to z, the triplet oi k x k hermitian 
matrices at each jumping point z = }Ia are the source terms and is the projection to 
the eigenvector corresponding to the eigenvalue exp(27rz^A) of the holonomy. A factor of 
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i in the definition of is included in order to make them hermitian. It is useful to fix the 
gauge such that the asymptotic Polyakov loop is diagonal, in this case pfai = i Im X^X^ 
with X^ meaning the A^^ row of A, although this will not be always assumed. 

Since 1)3. 3|l is a first order equation, the Dirac deltas on the right hand side give rise to 
finite jumps, 

Aj{}iA + 0) - Aj{iiA - 0) = ipf (3.4) 

in the dual gauge field at z = For this reason the |0-matrices are called the jumps. 

Along with the imaginary part, the real part of A^P^ A will also play a role and for future 
use we define the hermitian kxk matrices 

SA = ReA+PyiA. (3.5) 

The similarity between eq. ()2.11|l and H3.3|l should be clear by now and we recall the 
essential point once more. They express the fact that the dual gauge field, for instan- 
tons and A^ (z) for calorons, satisfies the dimensionally reduced self-duality equation with 
source terms, or equivalently the fact that the appropriate D^D operators are real quater- 
nions. For instantons the reduction leaves only a point and the source is simply the imagi- 
nary part of A^A whereas for calorons the reduction leaves a circle and the source is given 
as the imaginary part of A"*" f*A^(z — }Ia)X. 

It is worth pointing out that the Nahm equation H3.3|l also appears in the study of BPS 
monopoles on R"^ but in that case the range of z is an open interval |26H27l|59|. This follows 
from Nahm's duality as three periods of have to tend to infinity and one has to shrink 
to zero in order to have W?. In the dual description this means that the dual torus reduces 
to R, as we have seen in section IZ2l The finite jumps in the dual gauge field are the same 
for BPS monopoles if the rank of the dual gauge group is the same before and after the 
jumping point. Having only finite jumps for the caloron corresponds to having no over-all 
magnetic charge, but for pure monopoles this is of course impossible, hence in this case 
the boundary conditions are different and in particular involve poles for A;(z). When we 
discuss in sections 13.4.21 and 13.4.31 how BPS monopoles are embedded in calorons we will 
show how these poles arise as boundary conditions from the caloron point of view. 

We have seen that a simplified variant of the ADHM construction exists for the sym- 
plectic series which for Sll{2) = Sp(l) makes practical computations swifter. The require- 
ment of B^i being real and symmetric translates into 

A^(-z) = A,.(z)^ (3.6) 

The A matrix in this case is a fc- vector of quaternions, thus can be written A" = X^^Cpi with 
real coefficients A^. There are only 2 jumping points and we have the following restriction 
on pf and Sa, 

pf = -ph S[ = S2. (3.7) 
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The condition on the jumps is easily seen to be consistent with ()3.6|l . 

The first step in constructing the caloron solution is to solve eq. H3.3|l . The solutions in 
the bulk of the intervals {jiA, are solutions to the homogenous Nahm equation H2.29|) 

and the Dirac deltas give finite jumps across z = First we will investigate the general 
structure of Nahm's equation in the bulk of a fixed interval and then the structure of the 
matching conditions. 



3.1.1 Nahm's equation 

On each interval {}Ia>}^a+i) the dual U{k) gauge field satisfies the homogenous Nahm 
equation Ii27il 

A; + [Ao,A,] = ^e^jk[Aj,Ak], (3.8) 

which is the dimensional reduction of the self-duality equations to 1 dimension as we 
have seen in chapter |2l It follows that Tr A, is a constant. Furthermore, taking derivatives 
explicitly it is easy to see that Tr A,Ay — SijTr is also constant. More generally, 

Tr A/j . . . A/^jj is conserved as long as it is made totally symmetric and traceless in the indices 
z'l, . . . , im, giving rise to constant tensors 

Cii...i„ = ^, TrA,^ . . . A,„, + (permutations) - (traces) , (3.9) 

where the factor of (— z')'" was introduced to make them real. Note that only the first k of 
them are independent. Since they are totally symmetric and traceless, Ci^.„i„ is in the spin 
m irreducible representation of SO (3) and has 2m + 1 independent components. 

Totally symmetric and traceless combinations can be simply encoded using a complex 
null vector y, for which = 0. The conservation laws are then equivalent to the state- 
ment that y/j . . . y!,„Tr A,j . . . A/^^_ is conserved for any null vector y because the monomials 
Vh • • • Vim project on precisely the totally symmetric and traceless part. Obviously scaling y 
with an arbitrary non-zero complex number is irrelevant so it is best to think of y as living 
in a conic of CP^ defined by y\ + y\-\-yl, = 0. Such a complex submanifold is necessarily a 
CP^. The simplest way to see this is by the explicit parametrization of complex null vectors 
- up to an over-all factor - as 

, . (, 1^. M ) . ,3.10) 



Here ^ is a complex number which naturally lives on the Riemann sphere. 

The conservation laws thus can be expressed in a very compact form by saying that 
det (?/ -|- yA) is conserved for any complex ?/ and null vector y since expanding the deter- 
minant in Yj will reproduce the trace of any power of yA up to k. As a result the algebraic 
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curve in the variables {rj, ^) 

aetU^^^a.-e^^yo (3.11) 

is independent of z. This is called the spectral curve and has genus {k — 1)^. 

Let us count the number of degrees of freedom. Locally Aq can always be gauged 
away There are remaining real variables in the 3 anti-hermitian k x k matrices A;. The 
conserved quantities with m indices we have found have 2m + 1 independent components 
and summing them from 1 to gives a total of k^ + 2k conservation laws. After gauging 
away Aq, constant gauge transformations are still allowed, giving k^ — 1 gauge parameters. 
Thus we are left with 3k^ — {k^ + 2k) — {k^ — 1) = {k — 1)^ gauge invariant degrees of 
freedom to be determined. They are related to the {k — 1)^ globally defined holomorphic 
1-forms of the spectral curve [60J. 

We will see that the conserved tensors 1)3. 9|l determine the long-range or abelian prop- 
erties of the caloron and the remaining constants of integration are responsible for the 
short-range or non-abelian behaviour. 

The case = 1 is special because U{1) is abelian. There are no commutator terms and 
the Nahm equation simply says that A, is (covariantly) constant. Based on constant Nahm 
data the complete construction of the most general charge 1 caloron for unitary gauge 
group has been derived in 11571 l6Tll62ll63l . For an extension to arbitrary simple groups, see 

ill. 

Even if the topological charge is greater than unity one can look for special solutions 
which have constant Nahm data. Such an ansatz requires the constant A, to mutually 
commute leading to axially symmetric solutions for any charge [58 J. Here, however, we 
will be concerned with generic solutions with topological charge k > 1. 



3.1.2 Structure of the jumps 

There are n jumps in the dual gauge field for gauge group SU{n) and these will be related 
to locations of the constituent monopoles. More precisely, the jump in Tr A, at z = fi^ is 
zTr pf and since the x-dependence always enters as A, — ixj on both intervals, Tr pf can be 
interpreted as the shift between the center of masses of monopoles of type A — 1 and A. 
For this reason these traces carry a clear physical interpretation. 

Solution of eq. H3.3|l proceeds with solving the homogenous Nahm equation in the bulk 
of the intervals giving independent Nahm data on each. Then the differences at each jump- 
ing point, A;(^^ + 0) — A;(^^ — 0), should equal ipf. However, the pf are not arbitrary 
matrices, but are of a special form 

pfcri = iIm\^PA\, (3.12) 
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and any 3 differences Ai{}iA + 0) — Ai{}iA — 0) in general can not be written in this form. 
This fact is in contrast with the situation for unit topological charge in which case arbitrary 

- and necessarily constant - Nahm solutions on the two sides of the jump can be matched. 

In this section we will analyse in detail what the constraints are on the jumps pf for 
arbitrary topological charge - and consequently on the differences A;(^yi+ 0) — A, (/^^— 0) 

- and what their most general form is for fixed Tr pf. This will be of great practical help for 
the actual construction of the caloron because with prescribed center of mass locations we 
will be in a position to constrain the various moduli present in the most general solution 
of the homogeneous Nahm equation on each interval. 

We will be concerned with a fixed jumping point only, and the A index will be 
dropped in the remainder of this section. A diagonal holonomy will be assumed and the 
A^'' row of A will be denoted by ^. So ^ is a A:-vector of chiral spinors. 

There are Ak real parameters in ^ and 3fc^ real parameters in the three A: x hermitian 
matrices pi. There is a UiV) symmetry ^ e^'^^ that rotates the ^ but does not change 
the jumps pi, hence there are Ak — 1 moduli entering a jump at each jumping point. We 
know already that 3 of these moduli, Tr pi, are associated to constituent locations and our 
analysis will reveal the physical interpretation of the remaining 4fc — 4 parameters. 

Before doing so, let us mention a constraint that holds for the pi matrices as a direct 
consequence of their definition. For a complex null vector (I3.10|) we have det (i//T,) = 0. 
This means that i/;T; is of rank 1, thus can be written as the product of two chiral spinors. 
Concretely, 

y.T, = ^ J^Z] ) ( ^ + h ^-i ) , (3.13) 

where the fact that the two spinors are linear in ^ will become important in chapter|5l Now 
we only wish to point out that as a result, yipi has also rank 1 and can be written as the 
product of two fc-vectors. This constrains the jumps and they must satisfy the following 
quadratic equation, 

iViPif = yiyiaipj, (3.14) 

where we have introduced a, = Tr|0;. To keep flf fixed, we seek the most general form 
allowed by the constraints for the traceless part pi — ai/k and in particular its dependence 

on fl; . 

Spelling out all the indices in (13.12b for diagonal holonomy yields 

S'^H^f^-ptrf = ^^1^^, (3.15) 

where the real part S was introduced in 1)3. 5|l . It is clear that if the above equation is viewed 
as an identity for 2k x 2k matrices, then the requirement is that the left hand side should be 
a hermitian rank 1 projector. Hermiticity is guaranteed by S and pi being hermitian as T/ is 
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hermitian. A necessary and sufficient condition for a hermitian matrix Mpq to be of rank 1 
is that MpqMrs = MpsMyq. Imposing this condition on the left hand side of fl3.15|l and then 
decomposing the answer into its real and imaginary quaternion parts gives 

2s^bgcd _ S'^S'^ + pfpf (3.16) 

cad ^ch I ^adccb n^cdccih ^cb ^ad 

^ Pk +Pk ^ " "^Pk ^ — ''^ijkPi Pj ■ 

Further decomposition of the second equation with respect to symmetric and antisymmet- 
ric (cb) — (ad) indexpairs leads to 

gadpcb ^gcb^ad _ S'^pf + S'^pf (3.17) 
gcd^ab _ gabbed ^ ie^.^pfp^d, (3.18) 



Equations 1)3. 16113. 18)l are the necessary and sufficient conditions for the pi and S matrices 
to have the special form ()3.15|l with some ^. Now taking the trace in H3.18|l with respect to 
indices ad and ab gives two conditions that resemble a u{2) algebra, 

[S,pk] = i^iikpi^j, [Pi>Pj\ = i^ijk {spk - ^kS) , (3.19) 

where s = Tr S. In order to really see this introduce the normalized traceless part of pi and 
S as well as their normalized traces. 



(3.20) 

where a is the length of a,. Then changing variables from Eq, E; to Cq, Cf by 

Co = T^Eo-T^C-Ei (3.21) 
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translates (13.19b to 



[QXj]=ierjkCk, [Co,Q] =0, (3.22) 



which indeed shows that the Cq, C, matrices constitute a k dimensional representation of 
u(2). Contracting the indices ab in (13.151 1 shows that Cq > 1. We have assumed for a 
moment that eg > 1, but we will see that the limit eo ^ 1 is smooth and is relevant for the 
axially symmetric solutions L58J . 
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So far we have imposed condition ()3.18|l but not yet ()3.16|l and ()3.17|) . A lengthy but 
straightforward calculation of what these two conditions mean for the new C; and Cq ma- 
trices leads to the simple result for the Casimir operator 

TrQQ = ^, (3.23) 

which means that the k dimensional representation is the sum of a spin ^ and k — 2 trivial 
representations, hence there exists a unitary matrix U such that 

Q = ^Uf,U-\ (3.24) 

where f; is a A: x A: matrix with the usual Pauli matrices T/ in the upper left 2x2 block and 
zero elsewhere. Because it is so sparse only the first two rows of U will contribute to C,. 
Let us denote these two rows by the two fc-vectors for a = 1,2 and assemble them into 
ak X 2 matrix u. The index a can be thought of as a chiral spinor index. Since U is unitary 
the Ua vectors are orthonormal, u^u = 1. The lengthy but straightforward calculation we 
have referred to above also determines Cq and after putting everything together for the 
original variables pi and S we have 

Pi = - (sinh H{Sij — eiej) + cosh Heiej) UTiU + fl; 

1-1- wu' 
S = -fl/UT/U + coshH— ^ , (3.25) 

where we have introduced the new variable -ff > through cosh H = cq. These are our 
final formulae for the jumps pi and S, they must have these special forms. One can easily 
show by direct substitution that once pi is given, S is completely fixed as well, 

S = -3P2£i±ML=. (3.26) 
V3A/2Tr|0j|0j + 

Let us count the parameters. We have a, and H as 4 real moduli, in the orthonormal 
vectors 2k — 2 complex moduli but multiplication by a phase does not change pi nor S, 
leaving 4k — 5 real parameters, all together 4k — 1 moduli as it should be. 

Let us elaborate on our result. The u matrix is clearly a dual gauge moduli and apart 
from this and the already familiar a, there is only one more parameter, H. This means 
that once flf is fixed by putting the center of mass of the given type of monopoles to a 
prescribed location there is only the H parameter to be determined up to U{k) dual gauge 
transformations . 

Setting -ff = we obtain - after gauging away u - the simple expressions 

= |(^/'f/ + ^o) 
S = ^ (^^i'T; + to) , (3.27) 
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where tq is the 2x2 identity matrix in the upper left block and zero elsewhere, similarly 
to f/. This shows that the matrix structure of pi is the same for all i, namely ejfj + Tq, and in 
particular they mutually commute. This is precisly the situation of the axially symmetric 
solutions found in II58I . However, mutually commuting jumps do not necessarily imply 
axial symmetry. 

We also see that the jump in the dual gauge field A is limited essentially to a u{2) 
subgroup. Most of the components do not change across a jumping point if k is large. 
Another way of seeing this is by noticing that the number of parameters in the jumps 
grows linearly in k whereas the number of components in A grows quadratically. 

Above we have assumed that a ^ 0. The case of vanishing relative separation that 
corresponds to a = can be obtained by taking the limit a ^ and H ^ co while keeping 

exp H = s fixed. In this case we have - again after gauging away u - the simple forms 

s 

Pt = 2'^i 

S = '-TO. (3.28) 

These are the most general forms for the jumps and S up to dual gauge transformations 
for Tr Pi = 0. 

Once A is found over the dual circle, the n normalizable zero-modes of the anti-chiral 
Dirac operator in the background of Api (z) — ix^ have to be determined, from which A^ (x) 
can be computed according to the general prescription of the Nahm transform. Alterna- 
tively, as we have seen for the ADHM construction, the Green function can be used and in 
the next section we proceed with its construction. 



3.2 Green function 



The Green function is the inverse of A^(x)A(x) for the ADHM construction. Compactifica- 
tion from R"^ to x by the Fourier transform turns eq. (12.1811 into [27..57.,.58..63J 



d -^ 



it 



ix) +Y^S{z 

A 



^a)Sa 



f^{z,z')=5{z-z'), (3.29) 



where the k x k matrices Sa are defined by 1)3. 5|l for each jumping point z = The 
above form of the Green function equation should not come as a surprise following our 
extensive discussion of Nahm duality. On the 4-torus eq. 1)2. 7|l defined the Green function 
and we recognize the operator D^Dj of the (inverse) Nahm transform in the first two 
terms of eq. H3.29|l . The S-dependent - and consequently A-dependent - third term is 
the real quaternionic part of the source A"*^ — }Ia)^ and comes from the Fourier 

transformation of Re A"'^A. Naturally, fx{z,z') is periodic in z,z' just as the Dirac delta on 
the right hand side in eq. (13.291) . 
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By applying a dual gauge transformation if necessary, Aq{z) can be assumed to be a 
constant without loss of generality. Using subsequently the non-periodic gauge transfor- 
mation g{z) = exp {Aq — it)z it is possible to gauge away Aq — it completely and we define 

Uz,z')=g{z)f,{z,z')g{z')-\ (3.30) 

however doing so introduces periodicity only up to gauge transformation by the dual 
holonomy h = g{l) = exp (Aq — it), 

f,{z + l,z')=hf,{z,z'), fxiz,z' + l) = f,{z,z')h-\ (3.31) 

Along with fx one has to gauge transform accordingly all other quantities in the Green 
function equation H3.29|l and we introduce 

A,{z) = g{z)Ai{z)g{z)-'^ =exp{Aoz)A,{z) exp{-Aoz) 

Sa = giFA)SAgiFAy^ =exp{Ao}iA)SAexp{-Ao}iA) (3.32) 

thus Aj is only periodic up to gauge transformation by h similarly to fx- Note that the 
above adjoint transformation introduces time dependence neither into A, nor Sa- Clearly, 
since Aj satisfies Nahm's equation (I3.8| l with an Aq term. A; will satisfy one without it, 

A;='^eijk[Aj,Ak\, (3.33) 

which we will call the gauge fixed Nahm equation. 

The construction of fx follows the same methodology as the solution of Nahm's equa- 
tion on the full dual cirle. First we construct solutions on the bulk of each interval and then 
take into account the boundary conditions given by the matching conditions. 

3.2.1 Bulk 

After gauging away Aq — it every column w of fx satisfies the equation 



w = (3.34) 



in the bulk of each interval. As noted earlier, the operator on the left hand side - the 
covariant Laplacian in one dimension - can be written also as — D^D^ with 

Dx = ^ + CTj [Aj - iXj) = ^-Tj {iAj + Xj) 

D+ = -±-aj{Aj-iXj)=-^-Tj{iAj + Xj) , (3.35) 
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again invoking our general argument ()1.16|l . Thus it is natural to look for solutions of 
Dxip = 0. To avoid confusion we note that the global operator Dx including the jumping 
data on the whole circle as boundary conditions has no normalizable zero-modes as stated 
in section IZTI However, now we are seeking local solutions to eq. (13.341) and these are the 
same as local solutions to DxXp = 0. Although it may seem to imply the existence of zero- 
modes with the wrong chirality, these local solutions do not give rise to global, normalizable 
zero-modes. 

To solve Dxtp = we use the ansatz xp{z) = (1 + w,t/)w(z) where Ui is a z-independent 
unit vector and both u and w{z) may depend on x II60II . In fact the spinor part 1 + u,T; is 
a 2 X 2 matrix and hence represents 2 chiral spinors but one may take any column of it as 
they are linearly dependent as a result of |m| = 1. It is straightforward to show that once 
^(z) is annihilated by Dx then the corresponding iv{z) in the ansatz will satisfy the bulk 
Green function equation ()3.34|l . 

The Dx^p = equation for the ansatz leads to 



[lA + x)u 

dz 



w = (3.36) 
y{iA + x)w = 0, (3.37) 



where we have introduced the complex vector y = n^^^ + m^^^ with n^^\ n^'^^ and n'^'^^ = u 
forming a right-handed orthonormal basis of W?, in other words n^'^ x n^j"! = tij^n'^^'' and 

„(o„a) = 3^., 

Obviously there is a ambiguity in defining y from u but this ambiguity y e^'^y 
leaves H3.37|l invariant. Scaling y with an arbitrary non-zero complex number also leaves 
J3.37I I invariant, so it is best to think of y as an element in CP^. Once y is given, u can 
be reconstructed and the fact that and u are orthonormal is translated into the 

properties 

y X y J 
u = i- — r^, uxy = —ty, uy = , y =0, (3.38) 

so y is null. It already appeared in the context of the spectral curve H3.11|l . Its reappearance 
here is not an accident, eq. H3.37|l is an algebraic equation for iv and implies 

dety(zA-fx) = (3.39) 

for non-trivial solutions, which is exactly the equation of the spectral curve (13.11b with 
the identification t] = —iyx. It shows that the ansatz is self -consistent, the assumed z- 
independence of u and consequently of y is justified, since we have shown that the spectral 
curve is z-independent. 

The algebraic constraint (13.391 1 is a polynomial equation of order 2k for ^ if we use the 
parametrization (I3.10t and will generically have 2k roots. Taking its complex conjugate we 
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see immediately that if t/ is a solution, so is y. This means that the 2k roots for ^ naturally 
come in pairs. The transformation y ^ y induces the anti-podal map ^ —1/^ on CP^ 
and the transformation u ^ —u for u. We note in passing that the anti-podal map will be 
an important ingredient in the twistor construction of the moduli space in chapter |5l 

The polynomial (13.391 1 depends on x and this will specify the x-dependence of the roots 
^ and the corresponding y and u. Let us label the roots as ^''"K [^(''+^) = —i/^i") and 
accordingly y^''\ y^"^^^ = y'^"^ as well as u'''^\ u^"^^^ = —u^^h Restricting u and y to these 
special values makes eq. ()3.37|l solvable for non-zero w which as a result has to lie in the 
kernel of y{iA + x). The matrix 

M(z) = adj y{iA{z) + x) (3.40) 

projects exactly to the kernel, where adj stands for the adjoint - or matrix of subdeter- 
minants - defined by Madj M = det M. Thus the solution for iv is proportional to any 
column of M, 

Wa{z) = cp{z)Miz)ad, (3.41) 

for arbitrary d and with (p{z) a scalar function that is to be determined. This completely 
solves eq. H3.37|l and now we turn to eq. (|3.36|l . 

As a result of A satisfying the gauge fixed Nahm equation, it follows that^ 

M'(z) = [u{iA{z) + x),M{z)] . (3.42) 

Using this result in the substitution of ()3.41|l into H3.36|l gives for the function (p{z) the 
differential equation 

(p'Mcd = cp {Mu{iA + x))^^, (3.43) 
where on both sides arbitrary cd components can be taken as M has rank 1. Now using 

1 . 1 . 

Mu{iA + x) = -[M,u{iA + x)] + -{M,u{iA + x)} ^ 

= ~M' + ^{M,iuA} + uxM (3.44) 

where {, } is the anti-commutator, we see that (p{z) can be factorized as 

eMzux + xiz)l^ (3.45) 



^In fact for any holomorphic function h, the matrix H(z) = h{yA{z)) satisfies the differential equation 

H'{z) = [mA(z),H(z)]. 
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This defines the function ;^(z) which satisfies the rather simple equation 

X' = i^^, (3.46) 

which encodes all the information about the non-abelian cores as we will see shortly. 

The parity transformation ^ —1/^ or equivalently y ^ y or u ^ —u can be char- 
acterized by inspecting the large-x behaviour of the roots. For |x| ^ oo the polynomial 
equation (13.391) becomes (yx)^ = which has a fc-fold degeneracy and has roots y that are 
orthogonal to x. The vector u is orthogonal to y by definition thus we conclude that the 
solution for u in this limit is 

M = ±A/ (3-47) 

\x\ 

and due to the parity u —u exactly k roots are u = x/\x\ and the k others are u = —x/\x\. 
This analysis of the asymptotics gives us a natural characterization of the grouping of the 
2k roots into two sets of k elements, one corresponding to ux > 0, the other to ux < 0, at 
least for large enough \x\. We will arrange the 2k roots in such a way that m^^^x > and 
ii{b+k)^ < for large enough \x\. 

This asymptotic region has a transparent meaning also on the level of the bulk Green 
function equation ()3.34|l . For large \x\ the equation and its solutions tend to 

+ x^)w = 0, w{z) = e^\^\^wo , (3.48) 



dz2 

whereas eq. (13.361 1 and its solution tend to 

- MX ) w = , w{z) = e^^^wo , (3.49) 
az J 

for some constant fc- vector wq. This serves as a cross check that the roots u = ±x/|x| are 
indeed correct for large Also, we have established that the roots u = x/\x\ correspond 
to exponentially growing and the roots u = — x / |x| to exponentially decaying solutions of 
the bulk Green function equation in the large |x| limit. 

To package all the solutions, introduce the two k x k matrices (z) with containing 
as columns the k exponentially decaying and containing as columns the exponentially 
growing solutions w. Explicitly, 

/+(z) = (pW(z)M(?(z) 

fabi^) = <p('+'Hz)Mi?+'^)(z), (3.50) 

where the column c can be arbitrary and the superscripts (b) and {b + k) mean that the 
corresponding roots should be plugged into cp{z) and M(z). 
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Let us summarize what has been done. We started off looking for solutions of the 
bulk Green function equation ()3.34|l which is a second order differential equation for k 
variables so we expect 2k solutions. We have produced exactly 2k solutions packaged in 
f^{z), where cp satisfies (13.43b , M is defined by (13.40b and the number 2k comes from the 
2k possible choices for y and u. In particular the z-dependence of the solutions are quite 
similar but one should resist the temptation to think that they are - as solutions - the same 
or linearly dependent. Formulae ()3.50|l represent 2k linearly independent, and as such the 
full set of solutions. The problem of solving ()3.34|l in full generality is thus reduced to 
finding roots of a polynomial and solving the first order differential equation 0.46b in a 
single variable. 

For the subsequent discussion on various asymptotic regions it will be useful to intro- 
duce IIS3 

R^{z)=±f^'{z)f^{z)-\ (3.51) 

which is easily seen to be independent of the initial condition for As a result of f"^ 
satisfying ()3.34|l a Riccati type of equation holds for R^, 

±R^'(z) +R±(z)2 = (zA(z) (3.52) 

From the large |x| behaviour of the columns of given in (I3.48t one can deduce that 
R'^(z) — *• |x| for large \x\, which is indeed consistent with the above Riccati equation. 

The reason the quantities are useful is that they only contain algebraic dependence 
on X and when isolating exponential terms from algebraic ones they will naturally arise. 
The roots y are certainly algebraic in x, so is the matrix M = adj y{iA + x) and all the 
exponential dependence is contained in cp. 

Let us introduce yet two more matrices f ^ (z) that take this into account by leaving out 
the cp factors from thus are simply columns of M*^^) (z), 

Fabi^) = Mi?(z) 

f-,(z) = mS+'^^Iz), (3.53) 

which are clearly algebraic in x |65|. Again, the index c is arbitrary. The definitions (I3.51t , 
the fact that the cp factors drop out and the labelling convention lead to 

^.1 = + x),,f ±f ± (3.54) 

(note that the index b appears 3 times and is summed over from 1 to k) which shows 
directly that R^ has only an algebraic dependence on x as well as on A. This means that 
once A is known no further integration is needed to determine R^. It is also clear that as a 
result of u^^^x \x\ for all b in the |x| —>■ oo limit the above formula leads to R^ \x\ as it 
should. 

After this short intermezzo on R^, to be used later in section l3!4l we continue with the 
construction of the Green function. 
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3.2.2 Matching at the jumping points 

The boundary conditions - periodicity in z, jumps atz = ji^ and the behaviour of fx {z,z') 
as z — > z' dictated hy S{z — z') in H3.29|l - have not been taken into account yet. 

In order to incorporate them properly note first that the terms S{z — }iA)SA give jumps 
in the derivative of fx{z,z') as the equation is second order. It will prove useful to switch 
to a first order equation in the standard way by assembling fx{z,z') and ^/x(Z/Z') into a 
2k X k matrix. It satisfies the equation 



rf_ _ / 1 

dz \ -{Aj-ixjf + ZA^iz-}iA)SA 



Mz,z') \^_f 
ifx{z,z') J ySiz-z') 



(3.55) 



and as a result has jumps in its second component. As this equation is first order, the solu- 
tion is a path ordered exponential with known evolution between the jumps, this we have 
computed in the previous section giving matrices (z) for each interval. This evolution 
for Ha <z,z' < }Ia+i is given by the 2k x 2k matrix 



For the full evolution from jia to jia+i we write Wa = W(^^_|_i,^^). Due to fx being 
only periodic up to a dual gauge transformation, we have Wa+u = hWAh^^- Here and 
henceforth the products of the type HWa where ak x k matrix multiplies a 2k x 2k matrix 
is understood k x k blockwise. 

As z reaches any of the jumping points J4-a one has to insert into the evolution the 
2k X 2k matrix 

^^ = "^p(|, o) = (/^ ?)' ^^-^^^ 

responsible for the jump in ■^fx{z,z'). Thus the full evolution from < z' < piB+i to 
jiA < z < }iA+l> which range may now contain any number of jumping points, is given by 

W(z,z') - W(z, liA) }aWa-i }a-i ■ ■ ■ Wb+i/b+iW(^b+i,z') . (3.58) 

Incorporating the inhomogenous right hand side proportional to S{z — z') is easy, it gives 

( ) = W(.,.„) C(.„,.') - W(z.z-) ( « ) , (3.59) 

with an arbitrary 2k x k matrix C(zo,z'), arbitrary initial point Zq and the step function 
6{z — z') whose derivative is ^(z — z'). 
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All what is left is taking into account the periodicity of fx{z,z') under z — z + 1 and 
z' ^ z' + 1. This requires a little care as we have gauged away Aq — it completely which 
introduced into fx{z,z') periodicity only up to gauge transformation by the dual holon- 
omy h, see (13.311) . This requirement determines C{zq,z') to be 

Substitution into H3.30|l and H3.59|l leads to the final solution for the full, periodic and con- 
tinous Green function 

Mz,z') = g{z)-' [-W{z,z') + e{z' - z)W{z,z')]^^g{z') , (3.61) 

where g{z) = exp {Aq — it)z, the subscript 12 means that from the 2k x 2k matrix the upper 
right kxk block should be taken and we have defined the important contribution, 

r(z,z') = W(z,z„) w(z„,z') . (3.62) 

The choice of zq is arbitrary, different choices are convenient for different applications of 
formula ()3.61|l . The factor W(zo + l,zo) is the evolution operator over a full period and, 
just as any other W factor, is known explicitly in terms of (z)'s which we constructed in 
the previous section for each interval. All the ^-dependence of the Green function is in the 
dual holonomy h = exp(Ao — it) and the factors g{z)^'^ ,g{z'). 

Generally the Green function is hermitian, however, if we use the variant of the con- 
struction specific to the symplectic series for SU{2) — Sp(l) the restriction 1)3. 6|l implies in 
addition that 

fx{z,z'Y =fx{-z',-z), (3.63) 

which, since }i2 = —^\, means 

fx{m,mY = fx{}i\,}i\) , fx{H,mY = fx{H,m) , fx{.m,^\f = fx{m,H)- (3.64) 



3.3 Gauge field 

One of the advantages of calculating the Green function is that the gauge field and its 
field strength can be expressed entirely in terms of fx and A. The general Sll{n) ADHM 
formulae (l2.19l - IZ20t can be directly used for the caloron and the Fourier transform gives 

(P{X)-^ = 1 -Y^PAXfxillA,m)^^PB, (pj{x) = Y^PA^(rjfxi}lA'FB)^^PB, (3.65) 
A,B A,B 
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which shows that we only need to know the Green function evaluated at the jumping 
points. Fourier transforming the ADHM formula ()2.19|l of the action density yields L58J 



-lTrf2^ = -inniogdet #-(20,20)"', ^(zo,Zo)"' = 1 -h-^]N{zQ + l,ZQ) , (3.66) 

which expression is easily seen to be independent of 2o due to the determinant. 

We have noted that for Sp(l) the matrix (p is proportional to the identity hence a scalar 
function of x and that A in this case is a fc-vector of quaternions, thus can be written A" = 
A^cr^ with real coefficient s A"^. The quantities (pj are arbitrary quaternions. Furthermore, 
we have properties H3.63l - l3^64)l . Now the following combinations appear in cp and (pj, 

A«A^ = AJA^, + (AJA^)^,, 
AVyA^ = AjA|;//j^ + ^(AVyA^+AV^A''), (3.67) 

where the identity 1)2. 2|l has been used. The decomposition into symmetric and anti-symmetric 
parts with respect to ah is useful because fxi^i, y-i) and fxiV-i, V-i) are symmetric and thus 
give zero when contracted with the anti-symmetric parts. If the asymptotic holonomy is 
parametrized by coi as P = exp(2m'(x',T;), then we have ^2 = = \co\ = co and we 
obtain 

(p-^ = l-Y^Tr M}Ia,Fa)Sa = 1 - 2Tr fxicv,cv)S2, (3.68) 

A 

where we have used the definition (13.51 1 of and the fact that Sj = S2 as well as (I3.64|t . 

Similarly we can calculate (pj using the decompositions J3.67II together with (I3.7|l and 
the orthogonal projections 

to obtain 

4>i = \ff{-^rCo)Pi{X''ajX^ -f aV,A'')P2 - h.c. -I- 2^Tr/;,(a;,a;)p2 . (3.70) 
2 CO ' 

The usefulness of these quantities is that the gauge field assumes the simple form, 

A^ = ^(Pf]'^ydy(Pj, (3.71) 

in exactly the same way as ()2.21|l for the ADHM construction. 
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3.4 Asymptotic regions 

In order to demonstrate the physical content of formula H3.61|l we define 4 sensible lim- 
its. First, the zero-mode limit in which all but one type of monopoles are far away from 
X. In this limit we drop all exponential corrections coming from monopoles not of the pre- 
ferred type but still keep algebraic corrections. All contributions from the preferred type of 
monopoles will be kept fully. Second, we define the SU{2) monopole limit in which case 
also the algebraic corrections from the non-preferred types will be dropped. This leads 
to quantities being sensitive exclusively to the preferred type, thus we should recover the 
known results for SU{2) BPS monopoles. Third, the SU{n) monopole limit, when only 
one type of monopoles is assumed to be far and all contributions - both algebraic and 
exponential - are dropped, but everything is kept from the remaining n — 1 types. This 
setting is the same as for SU{n) BPS monopoles. The forth kind of limit we consider is the 
abelian - or far field - limit in which case x is assumed to be far from every monopole. 
Exponential corrections coming from any of the non-abelian cores will be dropped and 
only algebraic tails will be kept. 

We emphasize that maximal symmetry breaking leads to all monopoles being massive 
and as a result to exponentially decaying cores thus justifying our procedures regarding 
these contributions. 

The structure of the Nahm formalism tells us that there is a notion of a nearest neigh- 
bour for the various types of monopoles. This is given by the ordering in group space 
originating from the ordering of the eigenvalues fii < fi2 < ■ ■ ■ < }in- Both in the Nahm 
and Green function equation only quantities from neighbouring intervals are directly in- 
teracting through the jumping conditions. Hence one may suspect that monopoles of type 
A are interacting with types A — 1 and A -|- 1 differently than with all other types. Since 
in our analysis we carefully separate algebraic and exponential contributions coming from 
the various types of monopoles, we will be able to determine in a precise way how the 
interation between the monopole types varies. 

3.4.1 Zero-mode limit 

We have chosen the name 'zero-mode limit' because this case will be relevant for the be- 
haviour of fermion zero-modes which are only sensitive to one type of monopoles, see 
chapter 01 In this limit all types but one are assumed to be far away and we neglect all 
exponentially small contributions coming from these, but still keep all algebraic contribu- 
tions |66|. Let the preferred type be associated to the interval (f/^, }Ia+i)- 

The basic ingredient is the evolution operator W(z, z ') defined by ()3.56|l . It can be writ- 
ten 
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for Ha < z,z' < }iA+i> where we have introduced 

Ui{z,z') = fi{z)fi{z')-' ■ (3.73) 

Since we have seen in (I3.54t that is algebraic, it is clear that only U^{z,z') contains 
exponential terms in W(z,z'). In particular U^{z,z') ~ exp ±|x|(z — z') for large \x\. Let 
us introduce furthermore LT^ without arguments to mean LT^ = U'^diA+l/j^A) and also 
Ra{z) = j{R^{z) + R^{z)). For the Green function we will need the full evolution 
operator W(zo + 1/Zo) where we are free to make the choice zq = + 0. In this case it 
becomes 

W{^A + 1/ Ha) = }A+n^A+n-l}A+n-l ' ' ' Ja+I^^A+iJa+I^A , (3-74) 

where /a+h = hj^h^^. The crucial simplification comes from the fact that for B 7^ A we 
drop all LZg terms as they are exponentially small, which leads to 

^« = K R+{Ib+i) ) ^BRBifiB)-' ( RbM, 1 ) , (3.75) 

WbJbWb-1 = l(^j^+^^^^^^^U+RB{}iBr^T.BU+_,RB-l{}iB-l^^^ 1), 

where we have introduced Eg = + Rg{}iB) + ^b_i(Hb)- Note that while Ra{z) only 
carries information about the interval {}iA,}iA+l)> the term Eg involves quantities from 
both (^B_i,f/B) and (;Mb,;Mb+i). 

Upon multiplying more and more terms in order to form ()3.74|l the pattern in H3.75|l 
continues with more and more U^'s in the middle. Eventually for the Green function 
H3.61|l the term (1 — h^^W{}iA + 1, ^a))^ is needed in the limit of large matrices for 
B ^ A. Evaluation of the limit is done using the following formula which can be checked 
directly, 

lhn(^l- (^l^K{c d)^ ' ^l-(^iyca + db)-'{c d) (3.76) 



for k X k matrices a, b, c, d in the limit of a large and generic k x k matrix K. Applying this 
formula for the evaluation of (1 — h^^W{jiA + 1, }Ia))^^ with all LTg for B ^ A contained 



in K and 



( b ) = ( + RUm^) )' ( ) = ( s'-^+i + - 1 ) (3.77) 



leads to 



Wiz,z') = (3.78) 
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1 



1 

A-1 



Here we introduced fc^-i = + (Ha) arid c^+i = S^+i + (Ha+i) which contain 
data from the neighbouring intervals. Using another identity that can be checked directly 
to hold for any matrices involved. 



1 



1)W,( 



-1 



-1 

ca+\ 



-1 ) 



-1 



ca+i 



H -1 



( 

we arrive at the final expressions for f/^ < < z < }Ia+\ 
f,{z,z') = (1, 0)W(z,^^+i)f 



( ca+1 1 ) 

( hA^l -1 ) W 



-1 

A ■ 



X 



X 



X 



-1 



( -1 ) 



X 



or writing out explicitly, 

/x(z,z') = 



X 



X 



(U^(z,^A+l) -JiA(2'?^A+l)ZA+l) X 



-1 



X 



where 



(3.79) 



(3.80) 



(3.81) 



(3.82) 



Let us analyse the final formulae (l3.80H3.8i1 l in more detail. We have started off with 
the approximation that any exponentially small contribution from monopoles not of type 
A are neglected. This means that the non-abelian cores of all monopoles except the pre- 
ferred type are pushed far away and may only contribute through their algebraic - or 
abelian - tails. The non-abelian core of the preferred type should of course be visible 
as no approximation was made regarding that. Hence the expectation could have been 
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that the monopoles of type A will be present fully in the final formula whereas all other 
monopoles contribute algebraically. The terms W{z,}iA+i), W(^^, and W{}iA,z') 
in H3.80|l contain data only from the monopoles of type A and the interaction with neigh- 
bouring monopoles of type A — 1 and A + 1 enter only through Ba-i and Ca in (13.801 1 or 
through and Z^^^ in (13.811 1 and are algebraic. Thus monopoles not of type A — 1, A 
or A + 1 contribute nothing at all in this limit, their effect is purely exponential, once these 
are dropped not even algebraic tails survive. 

3.4.2 SU(2) monopole limit 

In the SU{2) monopole limit we neglect all contributions coming from monopoles not 
of type A. In the zero-mode limit we have neglected the exponential corrections, now 
we have to drop the algebraic ones as well. These - as we have seen - only come from 
monopoles of type A — 1 and A + 1. We have also seen that these contributions are con- 
trolledby Ba-i = Sa + Ra-i (Fa) and c^+i = Sa+i + Ra+i (Fa+i) in (EHQJ and Zj^_ Z^^^ 
in (13.811 1. We have also established that R'^iz) \x\ for large |x| and thus -^(z) —>■ 
|y^_xl large center of mass location |yA-i| of monopoles of type A — 1, and similarly 
R^^-^(z) — > |yA+i| for monopoles of type A + 1. Hence pushing the two types A — 1 and 
A + 1 to infinity means that 

bA-i — > \yA-i\ — ^ oo, z\_^ — > 1 

— ^ Iya+iI — ^ ^\+\ — ^ 1- (3.83) 

Inspection of (13.801 1 or (13.8111 gives in this limit for /^a < z' < z < f/A+l 

1 

fx{z,z') = -UA{z,}lA+\)UA{}iA,^A+\y^UA{}iA,z') (3.84) 

^^'^a[z ) 

where we have introduced Ua{z,z') = U\{z,z') — U^{z,z'). This Green function can 
be computed from data on the interval {iiA,}iA+i) only as appropriate for an SU{2) BPS 
monopole. 

We note that the quantity Ua{z,z') in terms of which the Green function is expressed 
may also be defined as the solution to the bulk Green function equation ()3.34|l with the 
boundary conditions 

= Ra{z). (3.85) 

z'^z 

The boundary conditions for the monopole Green function can easily be obtained as 
well. Since Ua{z,z) ~ for coinciding arguments, formula (I3.84|l implies 

fxi^A+i,z') = fx{z,}lA) = 0, (3.86) 



Ua{z,z) = 0, 



d_ 

dz 



Ua{z,z') 



dz 



-U{z,z') 
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that is, the Green function vanishes at both ends of the interval {}Iai}^a+\)i which are 
precisely the right conditions for BPS monopoles [I26IJ. 

Let us now derive the energy density, eA{x), of monopoles. This will be calculated 
from the caloron action density (13.661) in a similar manner as the Green function above. In 
approximating det (1 — h^^W{}iA + 1, /^a)) we first keep exponentially large terms to 
highest order only and drop exponentially small Ug terms for B ^ A. Similarly to ()3.76|) 
the following formula may be checked directly to hold to highest order in a large matrix K, 

det (^1 - ^ ^ ]^ ( c d)^ — .det {-K{ac + bd)) , (3.87) 

which we will use with the same substitution as in (13.771) . This leads to 

detW{iiA,}iA)-'' det(-K)det (^( c^+i 1 ) , (3.88) 

which, upon applying the logarithm, shows that the action density is a sum of two terms, 
the energy density of monopoles of type A and another term that tends to zero alge- 
braically as the monopoles not of type A are pushed to infinity. Thus we arrive at 

CAix) = --AAlogdet[WA]i2 = --AAlogdet (LT^lf^A+i, ?^A)i^A(?^A)"^) , (3.89) 

where again the subscript 12 indicates that the upper right kx k block should be taken from 
the 2k X 2k matrix Wa- Since no approximation is made with respect to the monopoles of 
type A, formula (|3.89|l is the exact energy density of SU{2) BPS monopoles of charge k. It 
is a generalization of the well-known formula e{x) = — 1/2 AAlog (sinh(v|x|) /\x\) for a 
charge 1 monopole of mass 4:7tv, located at the origin. 
The asymptotic behaviour 

LTaIHa+I/Ha) = U\{^A+i,}iA) -U^{}iA+i>}iA) ~expvA|x|, (3.90) 

for large \x\ assures the mass formula / d'^xeA{x) = AnkvA, which, when summed over 
all monopoles and integrated over the of unit radius gives the correct action for the 
caloron, S = Sn-^kJ^VA = Sn-^k. 

We note in passing that it follows from our discussion of the general structure of the 
jumps in section 13.1.21 that once |yA±i| tends to infinity, the corresponding pf,pf^^ ma- 
trices also tend to infinity. And since they determine the boundary conditions for the dual 
gauge field on the interval {}iA>FA+i) ^ diverging jump will ultimately mean a pole for 
Af (z) atz = fiA and z = }Ia+i- Dimensional analysis shows that since pf and pf^^ diverge 
linearly in |yA±i I the divergence in A,(z) will be linear in (z — ^a)^ arid (z — }Ia+i)^^ thus 
it must develop simple poles at both ends of the interval |26|. This is how the appropriate 
boundary conditions for Nahm's equation describing BPS monopoles are recovered from 
the caloron. In this sense a full caloron configuration acts as a regulator for monopoles 
because with every type at a finite distance no poles are present anywhere in the compu- 
tations, only in the limit of pushing some of them to infinity. 
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3.4.3 SU(n) monopole limit 

In the previous section we made a detour to 51/(2) monopoles because they were swiftly 
obtained from the zero-mode limit which was actually our primary interest. Making a 
similar detour to general SU{n) monopoles is also possible by moving the n^^' type to 
infinity resulting in rather compact and explicit expressions for the energy density. 

A caloron has net vanishing magnetic charge and removing the n'^^ type - which, 
essentially, was neutralizing the caloron - leaves a configuration with magnetic charge 
{k,k, . . . ,k). The moduli space of such BPS monopoles and its metric for k = 1 has been 
fully described in [671 . 

Since the n'^^ type of monopoles to be pushed to infinity are associated to the interval 
{}in,}ii + 1) it is convenient to choose zq = }i\ for the evaluation of the Green function 
H3.61|l . Exponentially small LT^ terms are dropped as usual, as well as algebraic ones in 
order to have no effect at all from the removed n^^ type. 

Going through the same steps as before - using the key formula (13.761 1 to perform 
the limit of dropping exponentially small term with algebraic corrections still kept, then 
rewriting fx{z,z') analogously to (13.801) and finally dropping the algebraic and 
terms as well - we conclude that the Green function for SU{n) BPS monopoles of arbi- 
trary charge is 

Mz,z')=Wiz,lll)l2[W{fli,fln)l2r'W{fln.z')i2. (3.91) 

where }ii < z' < z < }in and again the 12 subscript means that the upper right x ?c 
block should be taken from the 2k x 2k matrices W. The term W{}i\,}in) = W(^„,/ii)^^ 
is the inverse of the full evolution operator over the full range of z from pii to fin- In this 
way we see precisely how the dual circle of calorons opens up into R appropriate to BPS 
monopoles. In the original setup this of course means that x M.^ reduces to IR'^ and our 
Green function (13.911 1 is indeed static. 

Regarding the boundary conditions for Nahm's equation and the Green function the 
same comment applies as for the SU{2) case. The jumps pj and pf at and }in diverge as 
a result of bn = §1 + R+(^i) |y„| — > 00 and c„ = S„ + R„ (^„) \yn\ 00, where y„ 
is the center of mass location associated to the pushed away monopoles of type n leading 
to a simple pole for the dual gauge field, Aj(z) ~ (z — f/i)^^ and A;(z) ~ (z — ^„)^^ as 
z ^ }i\ and z ^ Jin- The jumps at ^2/ • • • / Mh-i remain finite. 

Since W(^i,^i) = W{}in,}in) = 1/ the resulting Green function (13.911 ) has the property 
fx{}i\,z') = fx{z,}in) = as appropriate for magnetic monopoles. Now this Green func- 
tion is essentially given in terms of the f^{z) matrices computed in section 15 . 2 . 1 1 with the 
matching at z = ^2/ • • • / /^n 1 being taken care of by our general formalism. 

The energy density of a general SU{n) monopole of charge k can be obtained analo- 
gously to the SU{2) case in the previous section. We arrive at 

e{x) = -^AA log detW(^„, ^1)12. (3.92) 
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3.4. Asymptotic regions 



This concludes our detour to BPS monopoles and we continue with the abelian limit 
of the caloron Green function that will be used in the context of fermion zero-modes in 
chapter m 



f,(z,z')=U^{z,z'),^^^y:j^. (3.93) 



3.4.4 Abelian limit 

In this limit - also called far field limit - all exponential contributions are neglected and 
only algebraic tails are kept. Since in the zero-mode limit we have dropped already all but 
one type of exponential contributions, only one type is left corresponding to the interval 

For < z' < z < }iA+i the terms U^{z,iiA+i) and U^{}iA,z') are both exponentially 
decaying in ()3.81|l . leading to 

1 

'2R7{z') 

If z ' < z this result itself is exponentially decaying so we obtain a non-vanishing Green 
function only if z = z ' in which case it is 

Recall that Ra{z) = j{R^{z) + R^(z)) is purely algebraic in terms of solutions of Nahm's 
equation. Once A is known so is the Green function, no further integration is needed. 

It is also possible to determine the Green function for z = z' = jia or z = z' = }iA+i 
in the abelian limit. Following the same logic as before, U^{}iA,}iA+i) is exponentially 
decaying hence can be dropped, U^{}Ia>}^a) = U^{}Ia+i>'H-a+i) = 1/ and the remaining 
growing terms cancel identically in (I3.81t , which leads to the simple results 

11- 1 

fxiM'M) = i'^-ZA-i) ^j, ^ =^' fxiM+i'FA+i) = ^ , (3.95) 



2Ra{f. 



with Ha = Sa + Ra{^a) + Note that the Green function is not continuous 

at the jumping points. We have seen that fx{}iA — 0, /i^^ — 0) = (2R^_i(^^ — 0))^^ and 
fx{}iA,^A) = ^^4^ whereas fx{}iA + 0,^a + 0) = (2Ra+i(?^a + 0))"^ Naturally, the full 
Green function is continuous but in the limit we are taking - which can be considered 
to be the high temperature limit as well - it develops a discontinuity. The width of the 
transition from one side of }Ia to the other is inversely proportional to the temperature and 
in the limit of infinite temperature it becomes a finite discontinuity. 

In order to construct the gauge field the Green function evaluated at different jumping 
points is also needed. For z' = }Ia and z = }Ia+1/ however, the only remaining exponential 
term is {}iAr which is exponentially small, hence 

= 0. (3.96) 
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3.5. Solutions for SU(2) and charge 2 



As a result, the matrices (p, (pj in H3.65|l simplify considerably in the abelian limit, 

= 1-EP^AE;,1a+P^, cp^ = J^PA\crjT.-/\'PA. (3.97) 

A A 

These can be used directly to compute the gauge field as in (|2.19|l . 

The action density ()3.66|) simplifies as well. For this the choice zq = /^i is a convenient 
one and then we need to approximate det#'(^i, ^i)^^. Using the expressions (13.741 - 13^51) 
leads to 

W{m,m)-^= (3.98) 

where all terms are exponentially growing. Approximation of the determinant is done 
in the same way as in formula (I3.87| l using the substitutions. 



(3.99) 



yielding 



det#^(^i,^i)-i =ndetL/+;^-^%-, (3.100) 



^ det2J^^(^^)' 
finally leading to the following formula for the action density in the abelian limit. 



lTrF,^„ = -lDniogndetUjg^^^. (3.101) 



3.5 Solutions for SU(2) and charge 2 

For k = 2 the spectral curve is a torus and as a result elliptic functions appear in the 
solution of Nahm's equation, which are quite manageable for practical calculations II65I . 
The complication for > 2 is partly due to the fact that the genus of the spectral curve 
is greater than 4 and the generalization of elliptic functions to higher genus are not so 
elementary ||60|. It should be stressed that conceptually - and qualitatively - there is no 
major difference between k = 2 and k > 2, except for writing explicit formulae. Thus our 
detailed analysis of double topological charge is illustrative of the general multi-charge 
case and henceforth we set k = 2. 

The general strategy has been outlined before and we will use the symplectic formula- 
tion. First we solve Nahm's equation on a fixed interval, then deal with the matching at 
the jumping points and finally solve for the Green function. As the gauge group in this 
section is SU{2), the asymptotic Polyakov loop will be parametrized as P = exp(27rza;iT/). 
There will be two intervals and we have }i\ = -co and ^2 = where co = \w\. 
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3.5. Solutions for SU(2) and charge 2 



3.5.1 Bulk 

We have seen that the traces Tr A/ are conserved and essentially decouple from the su{2) 
traceless part. They parametrize the center of mass and henceforth Tr A, = will be as- 
sumed without loss of generality. We have also seen that the traceless part (in the i,j 
indices) of Tr A/Ay is conserved, so introducing f (z) = — ^Tr A;(z)A;(z) for the trace part 
we can write 

- ^Tr A,(z)A^(z) = Qj + F{z)S,j, (3.102) 

with the constant, symmetric and traceless 3x3 matrix C. Using the basis cry for su{2) the 
dual gauge field can be written as A;(z) = ~Xji{z)aj and the above relation translates into 
X^X = C + F. It follows from the gauge fixed Nahm equation ()3.33|l that 

1 » » » 2 

F ' = gSi/fcTr [Ai, Ay] Afc = -e^y^ Smnp XmtXnjXpk = 4 det X, (3.103) 

which together with X^X = C + F leads to 

f ' 2 = I6(f + ci) (f + C2) (f + C3) , (3.104) 

where the c, are the eigenvalues of C ordered as C2 < Ci < C3. In addition we have C1+C2 + 
C3 = 0. Now recall that the Weierstrass elliptic function ^(z; g2, gs) satisfies the differential 
equation L68J 

^'^ = 4{^ -ei){^ -e2)i^ -63) , (3.105) 

where the elliptic invariants are g2 = — 4(eie2 + ^i^s + £2^3) and ^3 = 4eie2S3- Thus the 
most general solution for F(z) is 

1 

f(z) = -^(z-zo;g2,g3), g2 = -64(ciC2 + C1C3 + C2C3) , g3 = -256C1C2C3, (3.106) 

with an arbitrary zq. This completely determines the z-dependence of the dual gauge field 
which we now reconstruct. 

It is possible to diagonalize C by an orthogonal matrix R and then it follows from 
X^X = C + F that 

X(z) = R diag (v/f(z)+ci, Vf(z)+C2, Vf(z)+C3) R'^ , (3.107) 

with another orthogonal matrix R. Clearly, R corresponds to spatial rotations while R is 
coming from a global 51/(2) dual gauge rotation U. Let us parametrize the eigenvalues of 
the traceless C as 

01 = 0^1^, C2 = d2^, C3 = d2^, (3.108) 
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3.5. Solutions for SU(2) and charge 2 



with a so-called shape parameter < k < 1 (not to be confused with the topological 
charge k which is 2 in this section) and scale parameter D > 0. In terms of these the elliptic 
invariants are 

g2 = ^D\l - k2 + k4) , g3 = ^D6(2k2 - l)(k2 - 2) (k^ + 1) . (3.109) 
It is useful to express the result in terms of Jacobi elliptic functions, 

m = ^. /2(2) = k'^, Mz) = ^, (3.110) 

■'^^ ' cnk(z) ' cnk(z) ' cry^{z) 

which were used in II69II in the BPS monopole context. Note that f\ and /s are symmetric 
functions of z whereas /2 is anti-symmetric in accordance with the symmetry properties of 
(Ji, hence the restriction 1)3. 6|l is fulfilled. We refer to standard definitions [68J, 

<P{z) nil 

dt X f dt 



1 - k2 sin^ i J Vl- k2 sin2 1 



K{k)= / ^ k' = Vl-k2 



snk(z) = sin^(z) , cnk(z) = cos^(z) , dni,(z) = y 1 — k2sn^(z) . (3.111) 

A useful identity is /■[ = /2/3 and its cyclic permutations. Then - using the relations 
F(z) -|- Cf = D^ffiDz) 1 4, which follow from an identity between the Jacobi and Weierstrass 
elliptic functions L68J - we obtain 

A,(z) = -ia, + ^Rijfj{Diz~zo))UajU\ (3.112) 

where we have added the arbitrary trace part a, . The arbitrary constant Zq will be fixed by 
the condition (13.61 1. We see that in order to fully specify a solution to Nahm's equation on a 
fixed interval we have to specify its location a, spatial orientation R, dual gauge orientation 
U, scale D and shape k, as well as zq. The condition 1)3. 6|l further restricts U to be generated 
by the symmetric 02. 

A special case is of interest. Taking k = 1 we see from (13.110b that /i(z) = 0,/2(z) = 
and /3(z) = 1 which gives, up to spatial and gauge rotations, 

—iui 

Mz) = I -ia2 I , (3.113) 

§0-3 - ias 

i.e. constant Nahm data used in [58] to construct axially symmetric solutions. Thus we see 
that the shape parameter k controls how much the configuration deviates from the axially 
symmetric arrangement. 
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3.5. Solutions for SU(2) and charge 2 



3.5.2 Matching at the jumps 

Having found the solutions for Nahm's equation in the bulk of both intervals one has to 
impose the matching conditions at the jumping points z = ±co. For definiteness we take 
the choice of equal mass constituents and set a; = 1 /4 as most appropriate for the confined 
phase where on average the trace of the Polyakov loop vanishes and we have 1 /2 Tr P = 
cos 2nco = 0. The jumping conditions will restrict the pair of moduli (fl,, R, U, D, k, zq) 1,2 
found in the previous section for each interval. For charge 2 and gauge group SU{2) we 
expect 4 constituent monopoles with arbitrary locations. 

We have found two sets of non-trivial solutions that interpolate between overlapping 
and well-separated constituents. It should be stressed that solving the matching conditions 
for the various moduli is, from a numerical point of view, a not very difficult task. One 
has to solve a finite number of equations - although transcendental - for a finite number 
of variables. In the present section we present exact solutions in order to illustrate in a 
controlled way how typical configurations behave. This investigation makes it clear how 
the jumping conditions restrict the moduli, which parameter controls which behaviour of 
the solution, etc. After all this is understood one can generate any number of numerical 
solutions essentially in a straightforward way. 

We will see in section I^TI when we discuss the behaviour of fermionic zero-modes that 
in the abelian limit they localize to ellipses. For large separation between the constituents 
these 2 ellipses collapse to 2 finite segments and the support for the zero-modes become 
4 points at the ends of these 2 segments. Since in this case the singularity structure given 
by the 4 Dirac deltas clearly signals 4 point-like objects with well-defined locations we will 
introduce the extremal points of the major axis of the ellipses as approximate constituent 
locations for arbitrary finite separation. 

For the first particularly simple parametrization - which we call 'rectangular' - the 
two disks are parallel and separated in height by a distance d. The asymptotic Polyakov 
loop is chosen not to be diagonal but rather P = exp {2nicoT2). 

The parameters entering the Nahm solutions on the two intervals are Di = D2 = D, 
ki = k2 = k, R^ = R2 = 1, Ui = U2 = 1, Aui =a]- a] = (0, -d,0) and zj = 0,zl = 1/2. 
We take Aq = 0, thus A/(z) = A,(z). Then the Nahm data on the first interval (-1/4, 1/4) 



IS 




(3.114) 



whereas on the second interval (1 /4, 3/4) we have 



Mz) = - D/2 ^5 fz- i n C72 + Zd 



(3.115) 
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3.5. Solutions for SU(2) and charge 2 



Now the matching condition at z = 1/4 gives the following jumps (remembering that for 
SU{2) the requirements for z = —cv follow from z = co, see 1)3. 7|l ) 



pj = d{ 1 + T2 ] , pl=rf| -1 + T2 I, (3.116) 





provided we impose d = D/2(D/4). This is an example of how the matching conditions 
constrain the various Nahm moduli entering the two intervals. The jumps above are seen 
to be of the allowed general form derived in section 13. 1.21 In fact the constraints discussed 
in that section were used to find these - and subsequent - solutions. All parameters, and 
in particular the separations between constituents, are determined by D and k. 
The configuration is characterized by the 4 approximate constituent locations 




In order to see how the point-like constituents of the axially symmetric solutions are 
deformed by the generic non-constant dual gauge fields, we have considered a second 
family that interpolate between axially symmetric and generic solutions. These we called 
'crossed' configurations because in the abelian limit the two ellipses are in the same plane 
and form a shape of an X as shown in figure 14.11 The interpolation is indicated by the 
arrows from the configuration of 4 constituents forming a cross to all being alined along 
the z axis. 

The moduli corresponding to this family of solutions is D2 = Di = D, k2 = ki = k, 

Afl; = of — a] = (0,0, — rfcos a), Zq = 0, Zq = 1/2, Aq = and 



JJ2 = LTf 1 = exp ^ , R2 = R-^=\ 1 1. (3.118) 



2 

The jumps are taken to be 



cos (p sin (p 
10 
sin (p cos (p 



/ — Tssina \ / Tssino; 

= d \ — T2sina: , = d — T2sina: | , (3.119) 

\ Ti + cos a J \ — Ti — cos a 

again in accordance with our results in section 13.1.21 for the most general allowed form. 
These data are subject to the matching condition A, (1/4 + 0) — A, (1/4 — 0) = ipj which 
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3.5. Solutions for SU(2) and charge 2 



means that the following constraints should hold. 



D sm{9 - cp) (^f3 (^j^ - (^^^^ = d(l-sina) 
D sm{e + cp)(^f 3 (^j^+f^(^^^^ = d(l + sina) (3.120) 

D/2 = dsina. 

The 4 constituent monopoles for this family of solutions are located at the approximate 
positions 




/1^ = ^\ I . (3.121) 





-l)^+''Dcos(/)- (-l)^dcosa 

In order to have exactly point-like constituents in the abelian limit we need to impose 
k = 1, implying sin a = and cos sin ^ = 0. The first possibility is cos 9 = for which 
I cos (p\D = d. Such parameter values lead to two constituents of opposite charge to coin- 
cide as can be read off from ()3.121|l . A solution of this type describes a singular (zero-size) 
instanton on top of a smooth caloron. Excluding this singular case we are left with the 
choice sin^ = for which | sin0|D = d, implying D > d. We now find axially symmetric 
solutions with constituent locations at 

y(j) = ±1 ((-l)^d + i-iro) 1^ j , (3.122) 

where the overall sign comes from the fact that cos a = ±1. For cos 7^ all constituents 
are now separated from each other, giving a regular solution. The jumps H3.119|l become 
in this case one dimensional and it can be shown that for SU{2) exactly point-like con- 
stituents, that is ki = k2 = 1, forces the jumps to be one dimensional for any choice of the 
remaining parameters. 

Nevertheless for k 7^ 1, or sin a ^ 0, insisting as before that equal charge constituents 
are well separated while keeping the centers of mass of these pairs at a fixed distance 
d cos oi, forces k — >^ 1 for increasing D. Thus this setting also tends to point-like constituents. 
We will illustrate this behaviour for 9 = 7r/4. In figure l3Jl we plot for a typical value of 
d the constituent locations as given by (13.1211 1, varying a from n to 0. Note that for a 
given d, a and 9 the constraints (|3.120|l can be used to solve for (p, D and k. The asymptotic 
behaviour for a. = 7r/2 is determined by 



k' = i|^t|^(l + 0(k'^)), D = 2y2d(l + 0(k'^)) 



^'1 / 

(3.123) 
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-0.2 -0.1 

-0J_ 

/\^^^'"^.2- 


0.1 0.2 







Figure 3.1: The locations of the monopoles and antimonopoles (fat vs. thin curves) in the 1-3 plane for 
what we called 'crossed' configuration. The arrows indicate that by varying a from tt to the configuration 
interpolates between the 4 constituents forming a cross and all being alined along the z axis. 

It is also interesting to inspect Ai (z) in the limit k ^ 1 (or D — > oo) in order to under- 
stand to which extent we retrieve the piecewise constant solutions of Nahm's equation as 
these give rise to exactly point-like constituents. To this end we plot /,(D(k)z) in figure 
13.21 which up to an overall scale and rotation represent the constituent locations. The plot- 
ted cases are for k = 1 — 10 and k = 1 — 10^, clearly showing how in the bulk of the 
interval /i 2 and fs^l which is the characteristic behaviour of the axially symmet- 
ric and hence exactly point-like constituents. We also see, however, that at the jumping 
points z = ±1/4 all three functions deviate from their near constant bulk value and de- 
velop a pole that scales as D^. This behaviour is in fact in perfect agreement with the 
general construction of BPS monopoles in which case poles are imposed on the Nahm data 
as opposed to finite discontinuities. 

3.5.3 Green function 

Having found complete solutions for Nahm's equation including the jumping conditions, 
the next step in constructing SU{2) calorons for topological charge 2 is to obtain the Green 
function. The general method was described in section I3.2[ and now we will apply the 
results to the present case. 

One of the ingredients were the zeros of the spectral curve. In terms of the parametriza- 
tion ()3.10|l one has to solve a 4^^' order polynomial equation for ^. However - using a 
simple trick that is specific to charge 2 - one can solve directly for y. The spectral curve 
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3.5. Solutions for SU(2) and charge 2 




Figure 3.2: The three functions /;(D(k)z) for -1/4 < z < 1/4, at k = 1 - IQ-"^ (left) and k = 1 - 10"^ 
(right), illustrating the approach to the point-like limit k 1. 

condition dety{iA + x) = and the fact that y is null yields 

ViVi = , (xiXj - -^x^ - Cij)yiyj = , (3.124) 

where the symmetric, traceless, conserved tensor Qy = — jTr A/Ay + -^Tr Aj^A]^ surfaced 
again. 

Let us denote the eigenvalues of the symmetric, traceless matrix X/Xy — Sij — by 
A; = A,(x) ordered as Ai < A3 < A2 and the orthogonal transformation that brings it into 
diagonal form by C = 0{x). Introducing y = O^y and u = O^u for the corresponding 
quantities in the new frame, the equations (|3.124|l reduce to 

y? + y2 + ys = , ^iy\ + ^lyl + Agyi = o . (3.125) 

These are easily solved and we obtain 

(V A2 - A3 \ 
(_l)«+i^^A^^ , y(«+2)^y(«), (3.126) 
zVA2 - Ai / 

in accordance with our general discussion on the parity transformation y ^ y m section 
13.2.11 For the subsequent formulae we set Zq = 0, a, = 0, R = 1, LT = 1 and D = 1 all 
of which can be reinstated at the end by appropriate scaling, spatial rotation and gauge 
rotation. In this case we have M(z) = adj y{iA{z) -\- x) = y{x — iA{z)), from which the 
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3.5. Solutions for SU(2) and charge 2 



explicit expressions for the matrices in H3.53|l are 

2[ 2xyW-y«/3(z) 2xy(2) - yf ^(z) 

^ _1 /yfVi(z)-^yf/2(2) yfVi(z)-#/2(z) 
n 2xy(3)-yf/3(z) 2xy(4)-yf/3(z) 



(3.127) 



with the choice c = 2 for the arbitrary index. Next we address solving (13.461) , which for 
charge 2 can be written as 

, _ ^ViUj Cjj 5cd - i{yx) {uA)cd (3128) 
reducing with the choice c = 1, = 2 to 

^ ^i/i-^-y2/2- ^ 

This equation was actually considered and solved in terms of {^-functions in the BPS 2- 
monopole context in 1 70 1 but we found it more convenient to express the solution in terms 
of elliptic integrals. One only has to use the properties (13.381 1, the facts that = 72/3 
and cyclically as well as ffi^z) — fjiz) — 4{ci — Cy) to show by direct substitution that the 
solution of (I3.129t with the initial condition xiO) =0 is 

zsign(z)^^-^ \nUf,\z),m) - n^hm) + \z\ 

4X1/ y3 L 

where the elliptic integral of the third kind is defined as 

nk(s,m)=/ , =, ni = 4:^^. (3.131) 

I (l-m^2)V(l-k2^2)(i_^2)' y2 

This solution for x determines cp through the definition (13.451) , which in turn specifies the 
solution for the bulk Green function equation, see ()3.41|) , and we obtain, 

Wa{z) = (p{z)M{z)a2=exp{zux + x{z))^^=^. (3.132) 

VM(z)i2 
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Finally combining our results we have the following exact solution for the bulk Green 
function equation. 



Wa{z) 



(x X 17)3 . . . X k' 1/1 V2 
exp I tz- ^ - zsign (z) 



4x1/ ys 



nk(/3-i(z),m) -nk(l,m) + |z| 



X 



-yi/i(z)-(-l)«zy2/2(z) 



These can be packaged into f^j^ (z) = w^a ^ and f^j^ (z) = wl"^^' where, again, the superscript 
(b) means that the explicit solutions ()3.126|l for y = Oy should be plugged into w. 



1/2 



2xy+(-l)W3(z) 

2xy-(-l)W3(z) 

ih+i) 



1/4 



(3.133) 



3.5.4 Gauge field 

Now that we have the exact Green function for topological charge 2 we are in a position 
to construct the gauge field, field strength and action density in a straightforward manner 
using the general formulae of section l33l 

To illustrate our exact solutions we have plotted the action density for two sets of pa- 
rameters for what was called the 'crossed' configuration. The plots are in the 1-3 plane at 
f = with holonomy OJ = 1/4. 

In figure 13.31 the parameters determining the solution are k = 0.997, D = 8.753, 
a = n/2, 6 = n/A, and (p = — 7r/4, see section 13.5.21 As k is sufficiently close the unity, 
corresponding to the point-like limit, the constituents are clearly visible. 

For an example of strongly overlapping constituents, see figure 13.41 This configuration 
has parameters k = 0.962, D = 3.894 and for the rest the same as in the previous case. 
Since monopole constituents of opposite charge are closer to each other then in the pre- 
vious case, the time dependence is also stronger. They are overlapping so much, that the 
action density does not even show 4 separate peaks. The plot is similar to the well-known 
doughnut shape [71J of coinciding static charge 2 monopoles, however, here we have a 
combination of 2 pairs of monopoles, and the caloron configuration is time dependent. 

We will see in section 14.21 that using the zero-modes, it is possible to isolate the con- 
stituents, even though the action density blurs them. 

The exact results from the previous section can be used for BPS monopoles too, this 
was shown in section 13.4.21 where we derived the monopole limit of calorons. In figure 
I3.5l we plotted the energy density (13.891 1 of a typical charge 2 configuration with parameter 
values k = 0.57 and D = 6.915, somewhere between the known doughnut shape 1171 1 
corresponding to k = and well-separeted monopoles at k = 1. 



3.5.5 Abelian limit 

In this section we make some comments on the gauge field in the abelian limit. We have 
seen in section 13.4.41 that upon dropping all exponential contributions coming from the 
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Figure 3.3: Field strength squared of an SU{2) caloron of charge 2. 



non-abelian cores of the constituent monopoles - essentially shrinking them to zero - the 
gauge field becomes abelian. Let us focus on the temporal component. It can be written in 
this limit as 

Ao{x) = ^^{x), (3.134) 

which is clearly abelian and defines the field <l>(x) (remember that in the present gauge 
the asymptotic value icOiTi is gauged away). Using our results on the Green function we 
can compute which we know in any case to depend algebraically on x. Let us recall the 
Green function evaluated at the jumping points in the abelian limit ()3.95|l , 

1 _ 1 

in terms of which we can calculate the quantities (p and (pj J3.68H3.70t to obtain. 



<p = -~ T, = 2^Tr^2^2^ (3-136) 
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Figure 3.4: Field strength squared of an SU{2) caloron of charge 2 with considerably overlapping con- 
stituents. 



where pj = exp(Ao a;) pj exp(— Aq co), see (13.321) . Now using (13.711 1 to express the gauge 
field by (p and (pj we arrive at 

O = 1 ^ , . (3.137) 

1 - 2 Tr 52^2^ 

Let us emphasize that in this formula for the asymptotic Higgs field of the caloron we know 
every term explicitly. Since for abelian fields linear superposition preserves self-duality 
we expect that O can be written as a difference <1> = Oi — O2, with each contribution 
coming from monopoles of type A. The minus sign is due to the sign change in magnetic 
charge. Such a factorization is not at all obvious from (13.1371 1. In [72J it was shown, using 
the twistor correspondence, that for BPS monopoles the algebraic tail of the Higgs field - 
in other words in our abelian limit - is harmonic almost everywhere with peculiar support 
on extended disks. If we indeed wish to identify the two parts (£>a with contributions of 
the two types of monopoles, then we should recover this result. We will come back to 
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Figure 3.5: Energy density of an SU{2) monopole of charge 2. 



this question in section 14.11 where we discuss the abelian limit of the fermion zero-mode 
densities. 
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Chapter 4 

Dirac operator and zero-modes 



The index theorem states that the anti-chiral Dirac operator D = —o'^{d^ + ) has k zero- 
modes for charge k calorons and D = o'ni^y. + A^) has none. The zero-modes of the former 
are used in this chapter as probes for constituent monopoles. 

The circle corresponding to finite temperature necessitates a choice of boundary con- 
dition. We seek normalizable solutions of D^^pz = which are perodic up to an arbitrary 
phase factor. 



in the gauge where the gauge field is periodic. The choice z = corresponds to peri- 
odic zero-modes most relevant for supersymmetric gauge theory, whereas z = 1 /2 is the 
canonical choice for the anti-periodic fermions of finite temperature field theory. 

It is well-known for instantons on that if the k lumps making up the instanton are 
sufficiently well separated then a zero-mode density ^p^^p will peak roughly at each lump. 
It is also well-known that for instantons both the gauge field and the zero-modes have 
algebraic dependence on the coordinate whereas for monopoles - due to their mass - the 
dependence in the vicinity of the non-abelian core is exponential. The same behaviour 
applies to calorons in the vicinity of each constituent monopole. Zero-modes localize ex- 
ponentially for z ^ }iA stated by the Callias index theorem [1251 and they delocalize - 
or spread - for z = jd^ with only algebraic dependence. Localization of the zero-modes 
for z ^ jd^ will be used as a tool to identify the constituents in a full caloron solution. 
In particular we will show to what extent can the constituents be identified as point-like 
objects. 

The k zero-modes may be packaged into a 2n x k matrix ipz and the straightforward 
generalization of the corresponding ADHM formula in ()2.19|l states that the densities can 
be computed from the Green function. 





(4.1) 



fz{x)Mx) 



1 



□/x(z,z). 



(4.2) 



47r2 
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with fx{z,z) given by formula ()3.61|l . Here □ is the 4-dimensional Laplacian. We have 
computed the Green function in several asymptotic regions - including the one relevant 
for exact BPS monopoles - in section|3]4|and now will use these results for the zero-modes. 
Exact solutions will be presented for SU{2) and charge 2. 



4.1 Abelian limit 

We have seen in section 13.4.41 that in the abelian limit the Green function becomes static 
and accordingly the zero-modes will be static too. The abelian limit can be thought of as 
the high temperature limit as well. Since the masses of the monopoles are proportional to 
the temperature, the non-abelian cores shrink to zero size in this limit. This gives rise to an 
infinite localization for the zero-modes, hence they can be used to trace the cores II66I . 

Quite remarkably the full zero-mode density - that is the sum over all k zero-modes 
- in the abelian limit will not depend on z as long as it is in the bulk of an interval to 
allow for exponential localization. It can then be expressed by the conserved quantities of 
the Nahm equation. For z = }Ia the zero-modes do not localize exponentially, but spread 
algebraically over both A — 1 and A types of monopoles to allow for the hopping between 
the two types. 

Let us recall the result for the Green function in the abelian limit for ji^ < z < }Ia+1/ 

/x(z,z) = ^^, (4.3) 

where Ra{z) = {R^{z) + R^{z)) /2 and the quantities K^(z) have been computed in sec- 
tion |3]2]T] and have been shown to have only algebraic dependence on x and A^, the solu- 
tion for Nahm's equation on {jiA, }Ia+i)- Thus for the zero-mode density we have 

^tMM'c) = ~,^J^. (4.4) 

where the 3 dimensional Laplacian A was used as we have shown that the Green function 
is static in this limit. A z-dependent dual gauge transformation may rotate the individual 
zero-modes, thus in order not to have this ambiguity let us sum over the k zero-modes 
and define "^(x) = TrRy^(z)^^. Now we may interpret Pa{x) = 2TzTrxpl{x)xpz{x) as an 
electrostatic charge distribution with potential 

Anix) = -4.npA{x). (4.5) 

We have seen for |x| ^ oo that Ra{z) \x\ which gives '^^(x) k/\x\ (k coming from 
the trace of the identity matrix) implying that the total integral of Pa{x) over 3-space is 
the topological charge k. This can also be seen from the normalization of the zero-modes. 
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/ d^xtpltpz = 1, which when summed over all k of them and remembering that /3 = 2n also 
gives / d^xpA{x) = k. 

Following our discussion on the infinite localization of the zero-modes we expect 
to be harmonic as a function of x almost everywhere. In addition we expect it not to depend 
on z since the cores on which |0^(x) is supported are determined by the gauge field, which 
clearly does not depend on z. 

The zero-modes of a class of axially symmetric solutions with arbitrary charge |58| and 
the most general charge 1 solution l6Tll62ll63ll can be obtained easily and will illustrate most 
of what we have said above. These cases are special in the sense that they can be described 
by a constant dual gauge field on each interval. Consequently the constant matrices A, 
must mutually commute on each interval and hence can be diagonalized simultaneously, 
yielding on 



A = zdiag 



(4.6) 



defining all together nk vectors y^. From these, k are associated to each each interval 
(^^, }iA+i)> representing the locations of the k monopoles of type A. In this case the matri- 
ces (z) solving the bulk Green function equation (|3.34|l are particularly simple. 



/J(z) =diag (exp (^±z 



x-Ya 



,exp 



resulting in similarly simple expressions for J^^ (z) = ' (z)/^ 



±z 



R±(z) =diag ( 



(4.7) 



(4.8) 



which are also constant, but nevertheless satisfy the Riccati equation H3.52|l as they should. 
Now taking K^(z) = {R\{z) + R-j^iz))/! yields. 



(4.9) 



We see that Ya is independent of z, also that it is harmonic almost everywhere and the 
localization of the zero-modes for A fixed is on the k monopoles of type A. This is a clear 
signal that the constituents can be described as point-like objects for arbitrary separations 
between them in the present case. 

For generic higher charge solutions it is not directly obvious that the zero-modes follow 
a similar point-like behaviour. Due to their infinite localization we expect in any case a 
singular behaviour that is almost everywhere zero with a distributional support on the 
cores. In general we will find an extended disk-like structure, which, however becomes 
point-like for large separation between the constituents. The singularity is of course a 
byproduct of the abelian limit, in a full non-abelian solution everything is smooth. We 
have calculated the exact Green function in section 13.51 for charge 2 so we will see in this 
case how the singularity structure is resolved by the exact solution. 
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4.1.1 Multipole expansion 

In this section we will study in detail - for generic caloron solutions - the potential i^Ai^) 
and the charge distribution Pa{^) it gives rise to. Since we will be focusing on one given 
type, the index A will be dropped. 

In principle our formula (I3.54|l can be used to calculate R(z) = (R^(z) + i^^(z))/2 ex- 
actly and hence 'f as well. However, this way of analysing Y is a bit awkward as the posi- 
tion X is hidden in the roots ^ or null vector y of the algebraic curve in a rather complicated 
fashion. It is much more practical to use the Riccati equation H3.52|l and the asymptotic 
conditions R^{z) — > r = \x\ for large r we have found earier, 

±K±'(z) +R±(z)2 = (zA(z) = -A{z)^ + 2ireA{z) +r^, (4.10) 

where we have introduced the unit vector = Xi/r. One can assume that Tr A, = since 
any trace part can be incorporated by a shift in x/. Solving H4.10|l for a general solution 
of Nahm's equation on the right hand side is very difficult but an efficient expansion is 
possible in powers of r, furthermore one does not need to solve any differential equation 
in doing so. The reason is that in terms of the variable Y^(z) = R^{z)/r the Riccati 
equation H4.10|l becomes 

y±(z)2 ± y±'(z) J = 1 + 2iAe^ - A(z)2l, (4.11) 

and the term containing the derivative is only sub leading in 1 /r, resulting in purely al- 
gebraic recursion relations. The asymptotic condition for the new variable is Y ^ 1 for 
large r. 

It is a simple excercise to extract the first few terms of from ()4.11|l . Let us introduce 

the useful combinations U = —iAe and V = —A , in terms of which the expansion of Y"^ 
becomes 

11 11 1 

Y^ = l-U- + -{-U^+ y±U')^ + -{-2U^-U"+{U,V}± 2U^'t V')^ + ■■■ (4.12) 

where {, } stands for the anti-commutator. Wherever the derivatives U' or V' appear one 
has to use the Nahm equation to substitute some polynomial of A, into the expressions. 
From the above expansion one can compute Y(z) = 2(Y+(z) + Y^(z)) and the expansion 
of its inverse, the trace of which finally becomes, 

Tr ^ = + ^Tr {3U^ " + ^Tr {5U^ - 3UV)^ + ... . (4.13) 

Two important observations are in order. First, all the potentially z-dependent terms are 
actually constants as a result of A, satisfying Nahm's equation. For example, 

Tr{3U^-V) = {-i)^Tr{A,Aj){3e^ej-3ij) 
Tr (51/3 _ 3UV) = i-if Tr ( A^A^AO {Se^Cje,, - e,6^^ - e^d,^ - e^b^^) , (4.14) 
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with tensors built out of in such a way that they are totally symmetric and traceless. 
We have seen in section 13 . 1 . 1 1 that if Tr Ai^ . . . Ai^^ is contracted with such a tensor then the 
result is conserved. Thus, with the help of the constant, totally symmetric and traceless 
Qi...!m tensors defined in eq. (13.91 1 the expansion for 'f = r^Tr can be written as, 

r(x) = - + 3^ + 15?^^^ + .... (4.15) 

The dipole term is absent because we have chosen Tr A, = by shifting x, appropriately. 
The second very important observation we wish to make is that Y is harmonic. It is well- 
known that if Ci^...i,„ is a totally symmetric and traceless constant tensor then the mono- 
mials f^jgfj . . . defined on are spherical harmonics with eigenvalues —m{m -\- 1) 
for the Laplacian. Therefore each term C/j...,„,e,j . . . e,,,, / r"'+^ is harmonic and so is at 
least to this order of the multipole expansion. 

In fact one can continue the expansion to any given order, compute first Y ^ (z) from the 
Riccati equation (this is quite tedious, but only amounts to algebraic recursion relations), 
then y (z) as well as its inverse (also a very tedious expansion) and finally its trace at which 
point everything simplifies considerably. The pattern observed above will continue to hold 
and for definiteness we give the term with coefficient of in Y , 

^4 = ^Tr (^351/4 + imu" - 30^1/2 ^^y^ - y" ^^u'^^ = 105 Cijueieje^ei . (4.16) 

As mentioned above, the pattern continues to hold for higher orders as well and the coef- 
ficient of r^"^^^ will be proportional to (— z)*" Tr (A/^ . . . A,^^), contracted with a tensor 
with first term . . . e^^^ and remaining terms that are exactly such that they make it totally 
symmetric and traceless. Hence instead of {—i)^ Tr (A;^ . . . Ai^^) one can use the conserved 
tensors Ci^,„i,„, leading to ~ Qi...!,„^fi • • • ^i„,- The constant of proportionality is a trivial 
symmetry factor, the number of ways 2m indices can be contracted pairwise, yielding the 
remarkable formula 

~ (2m)! J_ 

m=l 

Here we have included an arbitrary dipole term CiCilr^ that was assumed zero by an ap- 
propriate choice of the coordinate system. Thus we conclude that the multipole moments 
of the total zero-mode density in the abelian limit are the conserved quantities of the Nahm 
equation. 

Let us summarize what has been done. We are not even close to solving Nahm's equa- 
tion for arbitrary topological charge, but it turned out that only the conserved tensors are 
needed in the abelian limit. The more complicated remaining (fc — 1)^ constants of inte- 
gration we have described in section 13. 1.1 1 are solely responsible for the non-abelian cores. 
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We are not able to solve the Riccati equation ()4.10|l either, yet we were able to derive an 
exact multipole series for y{x) in terms of the conserved tensors. This series is manifestly 
harmonic almost everywhere and independent of z. 

Let us repeat that the relevance of y^x) is that it determines the total zero-mode density 
in the abelian limit through 

'Yx^pl{x)ip,{x) = --La^W . (4.18) 

on 

Now let us check our formula against the special case of constant Nahm data where we 
did not need to resort to the multipole expansion. In this case, the diagonal A, given in 
(l4!6t leads to 

^^TrAu...Ai =— y"yf...yf, (4.19) 

which means that the conserved tensors are just the above monomials made totally sym- 
metric and traceless. Using one of the standard definitions of the Legendre polynomials, 
Pm, leads to 

where 0" is the angle between y" and e. Substitution into formula ()4.17|l indeed gives the 
well-known expansion ofj^u |x — y"!^^ in terms of spherical harmonics. 

Naturally, one would wish to sum the series ()4.17|l into a closed expression in the gen- 
eral case of non-constant Nahm data as well. One would then analyse its singularity struc- 
ture in order to determine what set the zero-modes are supported on. We will perform 
this analysis for charge 2, although by a different method, rather then by summing the 
multipole series. 

A final note about formula (|4.17|l . For general charge k only the conserved tensors 
with at most k indices are independent, the others can be expressed in terms of these in 
a polynomial way. This relationship between the multipole coefficients of Y results in a 
polynomial equation, 

E«-W^" = 0' (4-21) 



m 



where the coefficients am{x) are polynomials in x, and the independent conserved quanti- 
ties. For example, in the particularly simple charge 1 case the only invariant is C; = y;, all 
further tensors are polynomials of y^ and we have 

(x-y)V2-l=0. (4.22) 

For charge 2 the only invariant is C/y if we assume C; = and all higher rank C tensors can 
be expressed by Qj. In the next section we will derive the polynomial equation satisfied 
by y in this case. 
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4.1.2 Charge 2 

First note that the shape parameter k = 1 corresponds to constant Nahm data on the two 
intervals, resulting in 

1 1 

yix) = - r + - r, p{x) = S{x + y) +S{x-y) , (4.23) 

\x + y I \x — y| 

with y = (0,0, D/2). Now observe that the interchange C2 ^ C3 induces D iD and 
k ^ k' as can be seen from their definition (13.108b . Thus if we interchange the axis X2 and 
X3 and at the same time transform both D and k then "V will stay invariant. Therefore, 
from eq. ()4.23|l we can immediately obtain the corresponding result for k = 0. Strangely, 
the constituents moved to complex locations, something that has been observed already in 
the BPS monopole context 1.59 J . Nevertheless a manifestly real expression for the analytic 
continuation of (I4.23| l is 



'4r2 - D2 + J(4r2 - d2)2 + 1602^2 

r{x) = VS^ . (4.24) 

X4r2 - d2)2 + 16D2x2 



This function is singular on the circle of radius D/2 in the 1-3 plane, furthermore the X2 
derivative is discontinous on the disk bounded by the circle. This can be seen by expanding 
around small X2, but for 4(^2 + x2^ < 

which indeed reveals a discontinuity in the X2 derivative. We conclude that the charge 
distribution p{x) is singular on the whole disk, elsewhere it is smooth. This behaviour 
is very far from the exactly point-like situation and in order to understand how that is 
approached as k — 1 we now analyse the intermediate cases < k < 1. 

We know that R'^{z) satisfies the Riccati equation (13.521 1. However, not every solution 
of (13.521 1 will have the property that Tr R(z)^^ is independent of z. This means that R^{z) 
are very special solutions and since the equation is first order the whole solution is de- 
termined by the initial condition, which in turn has to be very special. Expanding (z) 
and A, (z) in Taylor series around z = 0, eq. (13.521 1 will determine all Taylor coefficients of 
R^(z) as a polynomial function of R^(0) and the Taylor coefficients of A;(z), which are 
polynomials in D and k. Now one can compute the Taylor series of TrR(z)-i and impose 
the constraints that all coefficients should vanish except for the constant term. This gives 
an infinite system of equations for R^(0) - of coure it is infinitely redundant - which ac- 
tually fixes it as an algebraic function of x, D and k. We immediately obtain y{x) since, by 
construction, ^{x) = TrR(z)^^ = Tr R(0)^^. Note that the choice z = corresponds to the 
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interval {-co, to), hence to monopoles of type A = 1, but by shifting z appropriately one 
can do the same for the interval {co, 1 — co), that is for type A = 2. 

Let us illustrate this method for unit topological charge. In this case A(z) = iy, a 
constant. The Taylor expansion of -R^ (z) - which is now a scalar, not a matrix - is taken 
to be 

(z) = a± + b±z + c±z^ + .... (4.26) 

Substitution into the Riccati equation (13.521) determines the coefficients b±, c±, ... as func- 
tions of J^^(O) = a± and we obtain, 

&± = ±((x-y)2-4) , c±=fl±(4-(x-y)2) , ... (4.27) 
giving rise to the following expansion of R{z)~^, 

+ 2 — z + 2 ^ ^ — z^ + . . . . (4.28) 



The requirement of z-independence implies fl+ = fl_ and {x — y)^ = a^a^, resulting in all 
further coefficients to vanish as well. Clearly, this leads to the correct result as it should, 

y{x) ^ |x-y|"^ 

This rather cumbersome way of obtaining y should in principle work for any topolog- 
ical charge, however, in practice we could only make use of it for charge 2 (apart from the 
trivial charge 1 case above). Even in this case it is very complicated as the first 11 Taylor 
coefficients in z are needed to constrain R^{0) unambigously. Leaving technicalities aside 
we only present the result of this procedure. Let us parametrize the initial conditions for 
the Riccati equation (|3.52|) by dimensionless quantities and Y^, 

R+{0) + R-{0) = D(Xo + X,T,) , R+(0) - R-(0) = D{Yq + Y,t,) , (4.29) 

and measure Xi in units of D/2. Then we have 

y = . (4.30) 

The constraints from the first 11 terms in the Taylor expansion imply that Yq = Yi = Ys = 
X2 = and that the variable Q = Xq — — X3 — Y2 satisfies the polynomial equation 

Q^+(k^-2 + r^) q2 - (k* + 2k^{xl -x\- x\) + ^x\ + r*) Q - + (2 - k})k^ + 
+4k2(x2 _ ;c2) _ k*(3x2 -x\- x\) + {{2 - k2)(3x? - x\^ x\) - 4x?) = 0, (4.31) 
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whereas Vi = Xi / Xq, V3 = X3/X0 and Yi can be solved for in terms of Q, 

r4 - q2 + 2(k'2(Q - 2x2) j^x\ + x\- x\) + (2 - k2)k2 



4xik'k2 



r4-Q2 + 2(Q-2x2-k'2(x2_;^2 ^2))_(2_k2)k2 ^ ^ 

= ^ — 5 ^-^^2 — — ^ — (4-32) 



^2 



r4 - Q2 + 2k2(x2 _ ;c2 - x2) + 4;^2 ^ 



4x2k' 

Thus we obtain from (14.301) the explicit expression for 'Y as an algebraic function of x, 

"V = , ^ , (4.33) 

2^J{Q + Yl){l-Vl-Vi) 

and the proper root of the cubic equation for Q is fixed by the asymptotic condition y 
II r, singling out the one with 

Q = r2 + 2x2 + k2(x2 -x2-x2) ^ ^^^^ ^^^^^^ 



This completes our derivation of the electrostatic potential Y for the most general charge 
2 solution. We have checked that the first 21 orders in the multipole expansion ()4.17|l agree 
with the closed form given above. In addition, we have checked numerically that the 
evaluation of Y using formulae H3.53H3.54|) also agrees with ()4.33|l . 

An alternative form for Y can be given as the root of a 6'^ order polynomial with coeffi- 
cients themselves polynomials in x and k. This is because Q is a root of a cubic polynomial 
and ^2 ig a rational function of Q, hence also fulfills a cubic equation. The coefficients of 
this polynomial are rather lengthy and are given in an appendix. 

Yet another way of writing Y , which turns out to be the most useful, is the following. 
Expanding it for small X2 yields, 

Y{x) = ^^[^^""'^ ^^3^2 + Oixl) , (4.35) 



D2k'2-4(x2^k'2x2jy 



for D2k'^ — 4 ^x2 _|_ k'^Xgj > g, generalizing the corresponding expansion for k = in 

(14.25b . In particular the cirlce of diameter D is deformed into an ellipse with major axis D 
and minor axis k'D. The singularity structure is similar to the k = case, the X2 derivative 
jumps at the disk bounded by the ellipse. It is useful to introduce polar coordinates {p,cp) 
suited for the ellipse in the 1-3 plane by (xi, X3) = (k'p cos cp, p sin cp). Taking the Laplacian 
then gives the following distribution, 

Pix) = -^Anx) = -HX2)^^^^2=M=. (4.36) 
^ 47r ^ ^ ^ ^ nk'p dp y/D^-4p^ 
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We see that due to the step function, 'Y is harmonic almost everywhere with support on 
the disk, including its boundary ellipse. We also see that there is a radially symmetric - in 
an elliptical sense - charge distribution on the disk. Further inspection of ()4.36|l shows that 
the total charge on the disk excluding its boundary is — oo, the total charge on the boundary 
ellipse is +00 such that the overall total charge is 2 as it should be. 

In section 15. 5.21 we have found two families of exact solutions of the matching condi- 
tions. In figure l4n we show the arrangement of the two disks for these families. The light 
and dark shading corresponds to constituents of plus and minus magnetic charge. One 
example of what we called the 'rectangular' configuration is on the left, a 'crossed' config- 
uration is on the right, with the curves indicating the approximate constituent locations as 
a is varied from n to 0, see section l3.5.2l for more details. 



y 




Figure 4.1: Two examples to illustrate the location of the disk singularities of the zero-modes for 511(2) 
and charge 2. What we called 'rectangular' configuration is shown on the left, the 'crossed' configuration on 
the right. 



The expression (14.3611 allows for the integral representation. 



D 

In T 



r{x) = - + — I dcp ^ =, (4.37) 

{ V^^-^P^'^P J{xi-k'pcos(pf + xl + {x3-psm(pf 



r n 



which can be used to check the k ^ 1 or equivalently k' ^ limit. For an arbitrary test 
function f{x) we have 



2n § 



j d^xp{x)f(x) = 2/(0) + ^ I dcp I -j^==^f [k' p cos cpAp sin cp) , (4.38) 
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which is smooth in k' and for f{x) = 1 verifies the proper normalization, / d'^xp{x) = 2. 
Evaluation at k' = gives. 



J d^xp{x)f{x) = 2/(0) + ^ y dx3 



T V 4 ^3 



dxi X3a3/(0,0,X3) 



D 



1+^3 JD^-4{xl + xl] 



2 

= 2/(0) + J dx3sign(x3)a3/(0,0,X3) = / (o,0,y^ + / (o,0,-^^ , (4.39) 

D 

which is indeed the expected result from (|4.23|l . We conclude that the rather peculiar dis- 
tribution of zero-modes that for generic k values is supported on an exteded disk smoothly 
tends to the point-like limit for k 1. In this limit the minor axis of the ellipse shrinks to 
zero while the major axis stays D, eventually tending to a segment of length D, with sup- 
port for the zero-modes at its two ends. We recall from ()3.123|l that for large separation D 
between the constituents k tends to unity, i.e. for large separation the constituents indeed 
become point-like. 



4.1.3 Higgs field and zero-modes 

We have derived a formula for the Higgs field <I>(x) of charge 2 calorons in the abelian 
limit in section 13.5.51 Its form that naturally emerged from the Nahm transform entangles 
the contributions from the 2 types of monopoles. Nevertheless, since linear superposition 
preserves self-duality for abelian fields, we expect that it can be factorized into two parts, 
each carrying Nahm data specific to one type of monopole. Now it was shown in 1721 that 
the algebraic tail of the Higgs field is harmonic almost everywhere and has the same kind 
of extended support on disks as the one we have found for the zero-modes. This makes 
us conjecture that the potential ^(x) for monopoles of type A actually equals the Higgs 
field of these monopoles in the abelian limit. Since we did compute the Higgs field for the 
caloron, see (I3.134|l and (I3.137t , we have verified the relation 

^{x) = -Tiix) - y2{x) (4.40) 

numerically and found agreement to high precision. This results in the exact identity 

Trxpl{x)xp^{x) = --La<I>^(x), (4.41) 

oTT 

for ji/^ < z < jiA+i- In other words, the fermion zero-mode density equals the abelian 
charge distribution of the monopole gauge field. Such a relation is at the heart of using 
chiral fermion zero-modes as filters to isolate the underlying topological lumps from rough 
lattice Monte Carlo configurations Il73ll74l 175117611 . 
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4.2 Exact results 

In the previous section we have seen that in the abelian Hmit, where exponential con- 
tributions coming from the non-abelian cores of the massive constituent monopoles are 
dropped, the zero-modes of the Dirac operator develop a singularity structure. In a full 
non-abelian solution, naturally, everything is smooth. We have computed the Green func- 
tion for charge 2 exactly, and using ()4.2|l we can now inspect how the singularity is resolved 
in the full solution. Exact zero-modes for unit topological charge have been computed in 

EzllZHI- 

We have plotted the action densities in section 15.5.41 for two particular cases. Now we 
will compare these with the behaviour of the zero-modes. The parameter values for figure 
14.21 are the same as for 13.31 but the scale of the zero-modes is enhanced by a factor of 5 
relative to the action density. The zero-modes are seen to follow the action density, with 
approximately coinciding peaks. This is due to having a k value close to unity, however, 
the configuration is not static, monopoles of opposite charge are close enough to produce 
time dependence. 




Figure 4.2: An SLr(2) caloron of topological charge 2, with the 2 periodic zero-modes on the left, the 2 
anti-periodic zero-modes on the right and the action density in the middle. 
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It is instructive to plot the configuration with strongly overlapping constituents as well, 
where the action density does not show the 4 constituent peaks. Nevertheless, the zero- 
modes reveal the presence of constituents. Since monopoles of opposite magnetic charge 
are closer than in the previous example, this configuration shows a stronger time depen- 
dence. The parameter values for figure l43l are the same as for l3.4[ 




Figure 4.3: An SU{2) caloron of topological charge 2, with considerably overlapping constituents. Even 
though the action density in the middle does not show the constituents separately, the zero-modes do. The 2 
periodic zero-modes are on the left, the 2 anti-periodic are on the right. 

We have derived in section 13.4.21 the Green function for pure magnetic monopoles as 
well and have plotted the action density of a typical charge 2 monopole in figure 13.51 In 
figure 14.41 the 2 zero-modes are shown together with the action density for the same con- 
figuration. 
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Figure 4.4: A charge 2 monopole for SU{2) as obtained from a caloron solution in an appropriate limit. 
The action density is in the middle and the two zero-modes on the sides. 
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Chapter 5 

Twistors and moduli 



The moduli space of self-dual fields over various flat spaces (monopoles over instan- 
tons over R^, calorons over x R'^ or Nahm-solutions over S^, etc.) carry natural met- 
rics. These are inherited from the L^-norm of the corresponding gauge field, restricted to 
modes transverse to the moduli space. The geometrical properties of the spaces over which 
the equations are defined are reflected by geometrical properties of the moduli space. In 
particular, both R"^ and x R^ are hyperkahler and as a result the moduli spaces will 
carry a hyperkahler metric. In the following section the main ingredients of the appara- 
tus needed will be summarized, followed by an elementary account of the hyperkahler 
quotient method. For a review on hyperkahler geometry see [79, 80J. 

We will be concerned with the moduli space ^ oi SU{n) calorons of arbitrary topolog- 
ical charge k and vanishing over-all magnetic charge with maximal symmetry breaking, or 
equivalently with non-trivial holonomy at spatial infinity, and its twistor space. 

As to the algebraic geometry of the moduli space we will derive a correspondence with 
stable holomorphic bundles over the projective plane CP^ that are trivial on two lines as 
opposed to one line as in the case of instantons. The difference is due to the different 
way R"* and x R'^ compactify to CP^. This correspondence is very much along the 
lines of geometric invariant theory as applied to instantons 181]. We arrive at it from an 
explicit parametrization of the moduli space by a number of finite dimensional complex 
matrices subject to a constraint, similarly to the quadratic ADHM equations for instantons. 
As this correspondence both for instantons and calorons is independent of the twistor 
construction, it is perhaps more fundamental than the usual correspondence between the 
moduli spaces and stable holomorphic bundles over the twistor space of the base manifold 
(or compactifications thereof). 

In general the (framed) caloron moduli space^ will be of real dimension Ank. This mod- 
uli space, just by naive counting of dimensions, can incorporate nk BPS monopoles of unit 

iWe will be only concerned with the framed moduli space that includes the moduli corresponding to a 
constant gauge rotation. 
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charge as every such monopole has 4 parameters, 3 for the location and 1 for a phase. 
We show that this is indeed the case, as in the limit of large separations the moduli space 
factorizes into nk copies of x M.^. For calorons of unit charge the exact metric was first 
conjectured in |53| based on string theoretic arguments and later proved in 1 57 1 for gauge 
group SU{2) with extension to SU{n) in 116 711 . 

5.1 Hyperkahler geometry and twistor theory 

Through the prototypical example of IR^ one can motivate twistorial ideas as follows 1*401. 
In two dimensions one can introduce complex coordinates z = x + iy thereby identifying 

~ C. Since S0(2) ~ there is essentially a unique way to do this in a manner 

compatible with the metric. Studying two dimensional problems in a single complex coor- 
dinate is a useful thing and one wants to study four dimensional problems using two com- 
plex coordinates. However, the situation in four dimensions is different because SO (4) 
is much bigger than U{2), in fact S0{4:) /U{2) ~ S^, which we will parametrize by ^. 
There is a whole family of choices to identify with C^, one of which is zi = xi + ix2, 
Z2 = x^ + ZX4. However, one does not want to single this out, for example any other choice 
with the coordinates permuted is just as good. Singling out a particular choice might spoil 
some symmetry of the problem or hide some features which are manifest in the formula- 
tion with 4 real coordinates. 

Thus one is lead to study problems in four dimensions by 3 complex coordinates 
{zi{Q,Z2{C),0 where Zi and Z2 themselves depend on ^ because they are the complex 
coordinates obtained by the identification that corresponds to the choice ^. If the original 
problem really lives on then the above description gives a parametrization of its twistor 
space CP'^ as a fibration over CP^ with fibre S^. This will be the general pattern as we will 
see. 

The idea of twistors originates from Penrose II82II and was invented in the Lorentzian 
signature context but has proved useful for the study of self-dual Yang-Mills fields in Eu- 
clidean signature too, as developed by Atiyah and Hitchin. The general idea is that differ- 
ential geometric structures on a smooth manifold are encoded in holomorphic data on its 
twistor space. For example in the case of Yang-Mills theory, an instanton solution on or 
S is encoded in a stable holomorphic bundle over CP3 with some certain extra properties. 
In general for a hyperkahler manifold its twistor space will encode the hyperkahler met- 
ric itself. Here we will describe only this aspect of twistor theory and neglect the original 
motivation of Penrose - general relativity. 

A Riemannian manifold is said to be hyperkahler if it is Kahler with respect to 3 com- 
plex structures and these satisfy the multiplicative relations of the quaternions. Thus such 
manifolds come with a metric g, 3 integrable complex structures I, } and K and the cor- 
responding 3 Kahler-forms coi = g{I-,-),C02 = g{}-r) and coo, = g{K-,-). The complex 
structures I,} and K obey = p = = IJK = —1 and are covariantly constant. It is an 
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elementary fact that the dimension of such manifolds is a multiple of 4. 

An equivalent characterisation can be given in terms of holonomy groups. Since par- 
allel transport preserves I, } and K the holonomy group of a 4n dimensional hyperkahler 
manifold lies in the intersection of 0(4n) and GL(n, H). The maximal such intersection 
is Sp{n) the group of n x n quatemionic unitary matrices, which group is one of the pos- 
sible holonomy groups on Berger's list [83J. The group Sp{n) is also the intersection of 
U{2n) and Sp{2n, C) hence a hyperkahler manifold is naturally a complex manifold with a 
holomorphic symplectic form. This observation will be heavily used in the construction of 
the twistor space and in practice will mean that many quantities will have a holomorphic 
dependence simplifying their study considerably. 

If M is hyperkahler and (^1,^2/ ^3) is a unit vector in R^, then = eil + 62} + e^K 
is again an integrable complex structure, if G CP^ corresponds toe G S^. Thus a hy- 
perkahler manifold is endowed with a whole = CP^ family of complex structures. It 
is advantageous to study all members of the family at once just as in the introductory 
example of R^. 

The twistor space of a 4m real dimensional hyperkahler manifold is defined to be a 
2m + 1 complex dimensional manifold Z together with a projection Z CP^ such that 
the fiber above ^ is, as a complex manifold, M endowed with the complex structure 7^. 
Topologically Z = M x CP^ but its complex structure is non-trivial and encodes all 3 com- 
plex structures of M. At a point {p,C) G Z the tangent space decomposes as TpM © T^CP^ 
and the complex structure is © i, where i denotes the standard complex structure of 
CPI. The projection Z ^ CP^ is in fact holomorphic. The 3 Kahler-forms combine into 
CO = CO2 + i<^3 + — C^{<^2 — 20^3)/ a holomorphic symplectic form with respect to on 
each fibre of Z. Because of the quadratic dependence on ^ (compare with (l3.10l - l3lTlt ) it 
takes values in i^{2). Here we use the same symbol <^(2) for the bundle over CP^ and its 
pull back to Z. Following standard notation, we denote by ff{k) the bundle on CP^ with 
transition function 1/^*^. 

On CP^ there is the antipodal map ^ which induces an antiholomorphic 

involution, or real structure on Z. Let us see how points in M are represented in Z. One can 
think of any point p G M as a holomorphic section {p} x CP^ of the projection Z CP^. 
They are obviously real, that is are left invariant by the real structure. One can show that 
the normal bundle of such sections is C^"(l). Here and henceforth we use the standard 
notation .^{k) for any bundle ^ twisted by ^{k) and for the trivial bundle with fibre V 
twisted by k we also simply write V{k). 

The crucial fact is that the above construction can be reversed. Suppose that the follow- 
ing are given, 

• a holomorphic projection Z CP^ where Z has complex dimension 2m + 1, 

• a holomorphic symplectic form on each fibre with values in ^(2), 

• a family of holomorphic sections each with normal bundle C^*" (1), 
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• a real structure on Z which induces the antipodal map on CP^, 

and the projection, holomorphic symplectic form and the family of sections are compatible 
with the real structure. Then the parameter space of real sections of Z is a hyperkahler 
manifold of real dimension 4m whose twistor space is Z. This is the precise statement how 
holomorphic data on Z encodes the hyperkahler metric of M. 

We will make use of the following result that is a direct consequence of the definition. 
If O is a flat hyperkahler manifold, then its twistor space is Z = 0(1). 



The self-duality equations - which are 3 independent equations - can be written as the 
3 components of a hyperkahler moment map set to zero. In this section we first review 
briefly the hyperkahler quotient construction and then demonstrate how it applies to self- 
duality. 

If G is a compact Lie group acting freely on a hyperkahler manifold M preserving the 
3 complex structures I,} and K, then it also preserves the 3 Kahler forms coj. Thus it is 
possible to define 3 moments maps, mj : M ^ q*, for each in the standard way [84J. We 
will usually identify the dual g* with q. One can show that the induced metric on the 
quotient 



is then also hyperkahler. This manifold is called the hyperkahler quotient of M with re- 
spect to G. 

Focusing first on the complex structure I, let us combine m2 and mo, into the complex 
moment map m = m2 + im^- It can be shown to be a holomorphic function with respect 
to I and furthermore to be preserved by the complexified group Gq- The corresponding 
complex symplectic form, co = C02 + icos, can also be shown to be holomorphic with respect 
to I. In fact, m is the holomorphic moment map with respect to the action of Gc and 
holomorphic symplectic form co. Then the quotient in (I5.lt can equivalently be constructed 



in other words a hyperkahler quotient can be realized as an ordinary symplectic - or 
Marsden-Weinstein - quotient but in a holomorphic setting [79J. The equation mi — 
is called the real equation, whereas m = is called the complex equation. 

It is easy to see explicitly how the two descriptions are equivalent. Setting m = is 
the same as setting m2 = = 0. In the holomorphic setting we solve m = up to G^, 
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(5.1) 



as 




(5.2) 
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whereas G^/ G can be used to solve mi = as well, since nii is not invariant with respect 
to Gc only to G. Thus solving m = up to Gc is the same as solving nii = m2 = ~ 
up to G. 

Let us count the dimensions. Setting the 3 real moment maps to zero reduces the di- 
mension of M by 3dim G, factoring by G then reduces the dimension by dim G, leaving 
over- all a manifold of dimension dimM — 4dimG. In the holomorphic description set- 
ting m = reduced the dimension by dim Gc = 2dim G, factoring by Gc further reduces 
the dimesion by 2dim G, resulting in again a dim M — 4dim G dimensional manifold, as it 
should. 

What we have done above was specific to singling out I. In the previous section we 
have seen that it is advantageous to combine the 3 complex structures and study them 
all at once. Let us now describe the hyperkahler quotient construction for any choice of 
complex structure J^, in other words we will consider it on the level of the twistor space. 
We have mentioned that coz = C02 + 2<^3 + 2^a;i — {C02 — icos) is holomorphic with respect 
to on each fiber of Z and we now analogously define 

mz = m2 + ims + 2^mi - ^^{m2 - im^) : Z — > 0c(2) , (5.3) 

as a 0c(2)-valued function on the twistor space of M. This complex moment map is holo- 
morphic on Z and we see that the choice ^ = corresponds to the discussion above when 
we have singled out I. 

Just as before, ooz is preserved by the complexified group Gc and we are led to the 
following description of the twistor space of the hyperkahler quotient of M, 

mzi(0)/Gc. (5.4) 

Clearly, ()5.4|l is a generalization of ()5.2|l which was valid in the fiber of Z over ^ = and 
now we have a description for the whole of Z. 

This concludes our brief review and now we turn to the (anti)self-duality equations on 
R^. First note that is flat and hyperkahler with hyperkahler structure 

1 . 

= 2^!'t/^^f' ^ • (5.5) 

By the same token the set of gauge potentials is also hyperkahler and flat, although 
infinite dimensional. Gauge transformations with gauge group G act as given in H1.3|l 
and preserve the hyperkahler structure. Thus we can invoke the hyperkahler quotient 
construction. The 3 moment maps are easily seen to be 

1 1 

= ^vlivFpiv = -Foi - -CijkFjk . (5.6) 

Thus setting = is equivalent to the anti-self-duality equations. The reason for being 
interested in anti-self-duality now, rather then self-duality, is simply that - following the 
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literature - we want to arrive at holomorphic quantities, rather then anti-holomorphic. 
The choice is a matter of convention and in this chapter we will stick to anti-self-duality 
and calorons will be anti-self-dual gauge fields on x R^. Time reversal interchanges 
self-dual and anti-self-dual gauge fields and leaves the moduli space metric invariant. 

First let us introduce fixed complex coordinates zi = xq — ixi and Z2 = X2 — 2x3 as well 
as a; = Aq + iAi and /3 = A2 + z A3. The moment maps in terms of these coordinates can be 
written, 

9 9 

m = m2 + im^ = -- — B - \ol,B] -\- - — a 

dzi dz2 

2imr = + + + J-iS + 5^iS' + [16, iS+]. (5.7) 

dzi dZi aZ2 0Z2 

Clearly, m is holomorphic in oc and fi and the set m = is preserved by G^-valued gauge 
transformations. On the other hand mi is not holomorphic and the set nii = is only 
preserved by the original G-valued gauge transformations. 

This completes our review of the hyperkahler geometry of the anti-self-duality equa- 
tions on R"^. These ingredients will be used for the construction of the moduli space of 
calorons and its twistor space in the next section where we discuss more details specific to 
our application. 



5.3 Moduli of calorons 

We have seen that the gauge equivalence class of caloron solutions is in a one-to-one corre- 
spondence with Nahm data ( A^, A) satisfying Nahm's equation on the dual circle, modulo 
dual gauge transformations. Let us recall how these dual gauge transformations g = g{z) 
act on the Nahm data, 

gAog-'-g'g-\ A^gAg-' (5.8) 
g(?^A)A^, (5.9) 

where \^ is the A^'^ row of the n x k matrix of 2-component spinors A. Simply dimen- 
sionally reducing to one dimension the formulae in the previous section we obtain a flat 
hyperkahler description for the dual gauge field A^. The space of all A is C^"^ which is 
also hyperkahler and flat. It is easy to see that the jumps ip^, which are quadratic in A, are 
the hyperkahler moment maps with respect to the action (15.911 at z = 

From the previous section it then follows that the moduli space of calorons is the hy- 
perkahler quotient of the space of all Nahm data ( A^<, A) with respect to the action (|5.8H5.9|I . 
This moduli space ^ will be our main object of study. 

Let us set G = U{k), Gc = GL{k,C) for its complexification and g = u{k) and gc = 
gl{k,C) for the corresponding Lie algebras. 
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The above quotient construction can be performed in two stages. First quotienting out 
by gauge transformations which are the identity at the jumping points z = followed 
by quotienting out by the remnant gauge transformations located at the jumping points. 
In the first stage the jumping data does not play a role as it is left invariant, hence one can 
concentrate on the dual gauge field A alone. On each (^^, }iA+i) interval it satisfies the ho- 
mogeneous Nahm equation subject to the boundary condition that at the endpoints A(^^) 
and A(/i^^i) are finite. Henceforth we will restrict our attention to this fixed interval only. 

This situation was analyzed in [85], see also 1861 : we will follow the same argument. 
Using the complex coordinates a. and /3 introduced in the previous section, dimensional 
reduction to one dimension of (15.71 1 leads to 

(a + ^e+y+Ka+] + [^,^+] = (5.10) 
/3' + [^^,/3] = 0. (5.11) 

We have seen that the space of solutions of equations ()5.10ll5.1l"|l modulo gauge transfor- 
mations by G is the same as the space of solutions of the complex equation ()5.11|l alone, 
modulo Gc gauge transformations f87|. This general principle of hyperkahler geometry 
can be shown explicitly to hold in this setting as follows. The complex equation H5.11|l can 
be solved by 

^ = (5.12) 

with a constant B = fi{}iA) ^ 0c ^rid Gc-valued 7(2) on which we impose 7(/^a) = 1 
for definiteness. Gc gauge transformations act on 7 as 7(2) — > g{z)y{z), hence the only 
gauge invariant quantity ish = ydi^+i ) (remember that we only quotient at this stage with 
gauge transformations which are the identity at the endpoints of the interval). Plugging 
this form into the real equation (|5.10|l gives for V = 7^7, 

(V'V-^^ ' + [B\ VBV-'^] = , (5.13) 

subject to the boundary condition ^(^a) = 1 and ^(f/A+i) = h'^h. Such a solution^ for 
V on the interval [jia, ^a+i] is unique and we can reconstruct 7 from V up to a G gauge 
transformation. In fact, V can be thought of as being Gq / G-valued, in accordance with the 
general discussion in the previous section. Hence indeed the moduli space for eqs. (15.101 - 
I5.11| l with G invariance is the same as the moduli space for the complex equation alone 
with Gc invariance. 

^Note that eq. iSlSi is the Euler-Lagrange equation for S = j'^^^^ Tr [\ {V'V-'^f + B^VBV-'^^ dz, that 
is a 1-dimensional WZNW-like model with a potential. 
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Si B2 Bn Bn+i = Bi 

hi h2 hn-\ hn 

Ml "2 Un Un+l = Ml 

Figure 5.1: Reduced Nahm data on the dual circle. 

In addition we have also obtained that this moduli space can be parametrized by {h,B) G 
Gc X 0c ~ r*Gc, the cotangent bundle of Gc. In terms of a. and fi they are given by 

/ Fa+1 \ 



h = Pexp — J a.{z)dz 

\ FA j 

B = fi{}iA). (5.14) 

The above discussion applies to each of the n intervals separately, giving reduced Nahm 
data {hj{,Bj{) for each interval labelled by A. Now we have to incorporate the jumping 
data. At each z = fi^ these are 2-component spinors in the fundamental representation of 
G, which can be cast into a k dimensional column vector u^, and k dimensional row vector 
Wa, in such a way that p2 + ip§ = u^® iv^- This is because we have seen in section (|3.1.2|) 
that whenever y is null, yipf is of rank 1 and now we take y = (0, 1, i) . 

These vectors make up the flat quaternionic space H"*^. In figure l5Jl the reduced Nahm 
data is summarized showing that the matrices B^ and vectors Ua, Wa are located at the 
jumping points with Ha giving the propagation in between. 

Now we proceed to the second stage of the reduction, namely quotienting out with the 
constant gauge transformations at the jumping points. These remnant complexified gauge 
transformations gA G Gc act as 

Ba ^ gABAgA^^ , hA > gA+lhAgA^ , Ua > gA^A , ^A ' ^ Ag~A > (5-15) 

with associated complex moment map equations 

Ba+1 - ^aBa^Ta = ^A+i ® wa+1 . (5.16) 
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At this point the moduli space M consists of elements in (T*Gc)" x H"'^ subject to eq. 
H5.16|l modulo the action of GJ^ given by ()5.15|l . with 

(rGc)" = rGc X ■ ■ ■ X rGc 

n 

Gl = Gc X ■ ■ ■ X Gc . (5.17) 

^ V ' 

n 

However, we can explicitly solve for n — 1 of the Bj{ variables using (15.161) and can 
gauge away the corresponding variables using the gauge transformations at }i2, ^3, 
. . . , fin- The remaining variables are B = Bi, h = hnh^-i ■ ■ - hi and some gauge transforms 
of the original Ua and Wa vectors which we will continue to label by and Wa and 
assemble into a k x n and an n x A: matrix u = (UaA) and w = (wAa)- The remaining 
symmetry by Gc is 

B — > gBg^^ , h — > ghg^^ , u — y gu, zv — > wg^'^ , (5.18) 
where g = gi- This action gives the complex moment map equation 

B-hBh-^ = uw . (5.19) 

The result of the present section is that the moduli space oiSU{n) calorons of topolog- 
ical charge k and zero overall magnetic charge with maximal symmetry breaking can be 
parametrized by matrices [h, B, u, w) subject to eq. (|5.19|l modulo the Gc-action (|5.18|) . 

Now recall that the ADHM construction gives the moduli space of SU{n) instantons 
on R"* (or S^) of topological charge k as the space of four k x k hermitian matrices B^, and a 
2-component spinor of n x matrices A^, subject to the quadratic equations 

[Bo,Bi] + [B2,B3] = 1(a+A2 + A+Ai) 

-1 

[Bo,B2] + [B3,Bi] = -(A+A1-A+A2) (5.20) 

[Bo,B3] + [Bi,B2] = ^ (a+Ai - A+A2) , 

modulo a natural U{k) action. Note a sign change relative to (|2.11|) as we are now describ- 
ing anti-self-dual, rather then self-dual fields. 

It follows from 1)5. 7|l by dimensional reduction to zero dimension that the ADHM equa- 
tions H2.11|l can also be seen as setting a hyperkahler moment map to zero. In terms of the 
complex coordinates 

A| + A+ 

a = ^ , a = Bn -I- iB\ 

h^^^^, /3 = B2 + zB3, (5.21) 
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where a isk x n and b isn x k, they turn into, 

[oi,^]+ab = (5.22) 
+ [j6,j6+] = 0, (5.23) 

where again (15.221 1 is called the complex, (15.231 1 the real equation. Applying the general 
principle of hyperkahler geometry to this case implies that the moduli space of instantons 
can be identified with the space of matrices {a,fi,a,b) subject to the complex equation 
[a;,/3] + ab = alone, modulo the complexified group GL{k,C) IISTII . 

What we have found is that the moduli space of calorons on x is almost the same 
space upon the identification 

a = h, jS = B, a = u, b = luh (5.24) 

except for the condition that one of the two k xk matrices should be non-degenerate as 
det/z ^ 0. Hence we have found an embedding of the caloron moduli space into the 
instanton moduli space as an open subset. The identification ()5.24|l gives a dictionary how 
to translate the 4n/c instanton moduli consisting of k scales and 4-dimensional locations 
plus some gauge orientations into the same number of caloron moduli consisting of nk 
3-dimensional locations and phases characterizing the nk constituent monopoles. 

5.3.1 Stable bundles on the projective plane 

Based on the above result a description in terms of stable holomorphic bundles on CP^ is 
possible. The relationship between stable bundles on projective space and instantons goes 
back to I3TII32. 88. 891. The original construction gave a correspondence between instan- 
tons and bundles on the twistor space of S^, which is CP^. The analogous correspondence 
for calorons was described in |56i giving bundles over the twistor space of x R'^. 

Later it was found in 118111 that compactifying the Euclidean 4-space to the projective 
plane CP^ instead of gives the instanton moduli spaces an interpretation in terms of 
bundles on CP^ without reference to twistor methods. In some sense this relationship is 
more fundamental than the twistor theoretic one. Our construction is a simple extension 
of fSTl and we summarize its essential ingredients below, for more details see 1901 . 

As said above, R"^ is compactified to CP^ by adding a "line at infinity". For any X G CP^ 
let [z\ : 22 : 23] denote its homogeneous coordinates. Then a monad is a sequence, 

Cx^ ^2k+n Dx^ (5_25) 

where the linear maps Cx and Dx depend linearly on z,, Cx is injective, Dx is surjective 
and in addition DxCx = 0. Because of the linear dependence on the coordinates we can 
write Dx = Z/D/, Cx = Z/C/ and the DxCx = condition gives six quadratic equations 
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DiCj + DjCi = 0. Every such monad defines a rank-n holomorphic bundle on CP^ by 
assigning to each point X G CP^ the fibre kerDx/imCx and this bundle will have second 
Chem class k. Conversely, it can be shown that any bundle on CP^, that is trivial on a line 
comes from a monad. 

The condition of triviality on a fixed line, say [zi : Z2 : 0], means for the matrices C/ and 
Di that D1C2 = — D2C1 is non-degenerate. In this case by an appropriate choice of bases 
for the three vector spaces in (15.2511 one can achieve that D1C2 = 1 and also that 




(5.26) 

Di = ( 1 ) D2 = ( -1 ) D3 = ( -|6 a ) 

for some a, /3, a and b matrices where the first two components of the C, and Dj are k x k, 
the last component of Cj is n x k and the last component of Dj is k x n. From the six 
quadratic constraints only one remains, namely [oi,fi] + ab = 0. This is the same as the 
complex ADHM equation H5.22|l hence showing the aforementioned correspondence be- 
tween instantons on and holomorphic bundles on the projective plane. 

Now we are in a position to determine what the extra condition deth 7^ means in 
terms of holomorphic bundles. Upon the identification (15.241 1 we obtain the following 
monad data 

C3 = ( I] 
\ wh J 

(5.27) 

D3 = { -B h u ) 

with Ci 2 and Di 2 as before. Using that D3C2 — — D2C3 = h and the same argument f Slll 
that led to the conclusion that triviality on the line [zi : Z2 : 0] means that D1C2 = — D2C1 
is non-degenerate, we conclude that the bundle has to be trivial on the line [0 : Z2 : Z3] as 
well. 

This is the result of the present section; there is a one-to-one correspondence between 
the moduli space of SU{n) calorons of charge k with maximal symmetry breaking and 
zero magnetic charge and the moduli space of stable rank-n holomorphic bundles on the 
projective plane having second Chern class k which are trivial on two distinct lines. 

A simple geometric picture clarifies how triviality on two distinct lines in the holo- 
morphic language and the base manifold x R-^ in the gauge theory language are re- 
lated. Two lines in the projective plane intersect in a single point. Our holomorphic bun- 
dle is trivial on both hence we can decompactify by removing them. The first removal 
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gives CP^ — CP^ ~ just as in the original instanton construction but now we have 
to remove the other CP^ as well. They intersect in a point that has been removed al- 
ready by the first CP^ leaving a CP^ — {*} — C behind from the second line. Now if 
we parametrize by (21,22) and remove the remaining plane (0,22) as well, we obtain 

- C ~ C* X C ~ X R3. 

Thus, what we have shown is essentially that the gauge potential A^(x), defined on 
S^xR^~C*xC, extends holomorphically to CP^. An analogous extension was proved 
for calorons in the twistor theoretic correspondence in |56|. We would like to emphasize 
again, that our construction - just as Donaldson's construction for instantons - does not 
rely on the twistor correspondence. 

5.3.2 Twistor space and spectral data 

So far we have identified the caloron moduli space M as a complex manifold but we have 
not said anything about the induced hyperkahler metric on it. This is encoded in its twistor 
space and has a convenient description in terms of spectral data, similarly to S\1{T) BPS 
monopoles as mentioned in secion lZ2.2l The spectral data will be derived from the twistor 
construction [79J which will be the subject of the present section. 

In order to find the twistor space J° one has to redo most of the first part of section l53l 
with an arbitrary choice of complex structure labelled by ^ G CP^ and trace the dependence 
on ^. This dependence can be expressed by transition functions from the patch LT = G 
CP^I^ ^ 00} to the patch y = g CP^I^ 7^ 0} under I, — ' t ^ \ since there exist 
holomorphic trivializations over both M and V. Quantities defined over V will be denoted 
by a tilde. 

It turns out that for this purpose the parametrization (15.18b is not very useful, it is 
better to go back to eq. (I5.16t with variables in (r*Gc)" x H"^ and symmetry GJ^. First 
we will describe the twistor space of r*Gc that corresponds to a given interval, then take 
n copies together with the jumping data H"'^ and carry out the reduction by GJ^. The 
hyperkahler structure on r*Gc and its twistor space have been described in 1.85 J . below 
we will reproduce the ingredients we need. 

To identify the twistor space of r*Gc recall that it came from an infinite dimensional 
hyperkahler reduction of Nahm data on an interval with regular endpoints. The vector 
space n of (a, j6) pairs is flat, hence its twistor space is 0(1) CP^. 

The fact that the twistor space is 0(1) means that one can pick the complex structure 
corresponding to ^ in such a way that (a, jS) change on the overlap of LZ and V according 
to 

Oi >• ^ ~ ^ 

/3 ^ ^ = ^. (5.28) 
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Complexified gauge transformations (which are the identity at the endpoints of the inter- 
val) act over U as 

^ — gh'^-^g'g'^ (5-29) 

and over V as 

S — > gH'^-tg'g'^ 

^ ^ gk'^-g'g'^- (5-30) 

The corresponding complex moment map equation (generalizations of eq. dS.llIt ) are easily 
computed over U and V to be 

/3' + K/3]+Ca' = 

a' +[di,^]+C^' = 0. (5.31) 



Their solution (generalization of eq. ()5.12|l ) over U is 

a. = — 7^7^"^ 

/3 = lBj-' + ^j'j-\ (5.32) 

and over V 

a. = 757^^ — 1^7 '7^"^ 

^ = 7'r\ (5.33) 

with constant matrices B,B and Gc-valued functions 7(2), 7(2) subject to the initial con- 
ditions "yijiA) = liV-A) = 1- From the latter we define the gauge invariant variables 
h = 7(^^_|_i), h = 'y{}iA+i), both in Gq, analogously to section l53l 

To specify the twistor space of T*Gc the transition function {h, B) {h, B) is needed. 
Substituting (a,^) from (l5!^ and {a, ^) from (l532t into (l533t gives 

7(z) = j{z)exp (^{z-fiA)^^ (5.34) 

where the factor of }Ia appeared in order to maintain 'y{}iA) = 1/ once j{}iA) = 1- From 
H5.34|l we obtain the desired transition function for {h,B), 

(KB) = (^hexp(^VA^y-^^ (5.35) 
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where = /^a+i — i^A is the length of the interval and is related to the masses of the 
constitutent monopoles of type A. Now we are able to identify the twistor space of T*Gc 
as 

Z 0c(2) CF\ (5.36) 

where the first arrow is a principal Gc bundle over the total space of 0c (2) with transition 
function exp {vj{B/Q and ^ is a section of this principal bundle, while B is a section of 
0c(2). 

The holomorphic symplectic form on each fibre of Z is the natural invariant 2-form 
CO — dTv (^Bh^^ dh) . Its transformation rule is 

CO = dTr (^Bh-'^dh^ = dTr (^^^Bexp (^-v^r^B^ h'^d (hexp (v^r^B^^^ = 

= r^dJr ^Bh-^dh^ + T^rfTr (exp (-v^r^s) ^^exp (v^r^B) b) = (5-37) 
= r^dJr [Bh-^dh^ +l,-\AdTr (BdB) = r^dTr (sh-^dh^ = r^cv , 

hence it is a globally defined ^(2) valued 2-form along the fibres, as it should be. 

The above discussion applies to every interval labelled by A separately and we obtain 
twistor spaces for each. Note that they are not identical, they have different transition 
functions as the masses of the constituent monopoles may vary. 

Now, parallel to section l53l we can incorporate the and Wa jumping data that makes 
up the flat H"'^ space. Again, because of flatness its twistor space is ]H"'^(1) — > CP^. 

The action of on the reduced Nahm data is the same as in (|5.15|l giving Q£ (2) valued 
moment maps which when set to zero give eq. ()5.16|l just as before. 

Upon quotienting with we recover the well-known spectral data of calorons IBSIISUII 
in the following way. Part of the gauge invariant quantities are the spectral curves Sa in 
^(2) associated to each interval defined by = {(C,7a) ^ ^(2) | det (7^ - B^(^)) = 0}, 
where Ba is a section of gc (2) and t]A is a section of ^ (2) . These curves have genus (fc — 1)^. 
Also gauge invariant are the sections over ^(2) 

= WA+ihAad] {f]A - Ba) ua ■ (5.38) 

Since {{hA, Ba), {ua, iva)) is a section of Za x ]H'^(1) with known transition functions, the 
transformation of tj^A under ^ — > | is 

rpA = r^''exp(^VA'^^xpA, (5.39) 

once the xpA is restricted to the curve Sa- In other words the xpA are sections of the line 
bundle L^^{2k)\s^ first introduced in L44i in the context of magnetic monopoles, see also 
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II60II . The transition function in H5.39|l can actually be taken as the definition of the line 
bundle U^{2k) over ^(2). 

The invariants ipA do not exhaust the list of gauge invariant quantities, the combina- 
tions 

^AC = WA+i^Aadj {f]A -Ba) - ■■ head] {tjc - Be) "C/ A > C (5.40) 

are all gauge invariant and obviously xpA = ^aa- One can think of the xp^c being 
associated to the interval {}ic, ^a+i) which for A — C > 1 contains several jumping points, 
while the xpA invariants are associated to a basic interval {jia, ^a+i)- Again, the transition 
functions for the invariants xpAC follow from eq. ()5.35|l . 

once they are restricted to Sq in each variable j/d for D = C, . . . , A. Here and throughout 
A — C stands for the difference modulo n. 

The jumping conditions H5.16|l impose constraints on the sections \pA- Specifically, 

^A+i^A = (Va+i - Va) ^A+hA , (5.42) 

where we have used the fact that i/?^ is a section over Sa, thus det {yia — Ba) = and 
similarly for A + 1. This means that over the intersection of neighbouring spectral curves 
Sa and Sa+i, where r]A = Va+i^ either xpA or xpA+i vanishes. Analogously to eq. (15.421) the 
jumping conditions (15.1611 give 

^A,D+i^DC = (Vd+i - Vd) ^ac (no sum) . (5.43) 
for A — C > 1. From the relations (I5.42ll5.43t it follows directly that 

A A-1 

Yl^D = ^pAC n (Vd+1 -Vd) , (5.44) 

D=C D=C 

in other words, if the interval (fiQ, }iA+l) contains several jumping points and hence can be 
broken into A — C + 1 basic intervals {}ic, ?^c+i)/ (?^c+i/ V-c+i)r ■ ■ ■ > {}^a> Fa+i) then there is 
a relation between the invariant xpAC and its "constituents" xpo associated to }iD+i) for 
D = C, . . . , A. In fact once ()5.44|) holds the relations ()5. 42115. 43)l follow. It is clear from eq. 
()5.44|l that locally the invariants xpAC can be expressed by the xpA> which implies that the 
spectral curves Sa and the sections xpA over them can be used as local coordinates for the 
twistor space 

More concretely, each spectral curve Sa is given by a polynomial of order k with leading 
coefficient 1, which may be parametrized by its roots rjAa- The sections xpA are restricted to 
these curves, thus can always be written as a polynomial of order k — 1, 

k 

^AaVA^ , (5.45) 

fl=l 
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with nk coefficents ^Aa- Then the parameters {^,f]Aa^^Aa) are local coordinates for the 
2nk + 1 complex dimensional twistor space iF. 

These serve only as local coordinates because there is an additional set of invariants, 
specific to the fact that for calorons the Nahm equation has periodic boundary conditions, 
namely the k invariants of the dual holonomy h = hnhn^i ■ ■ - hi. For these invariants there 
are additional constraints, which can be solved locally in terms of {i]Aa> ^Aa) hut not globally. 

The spectral data we have obtained from the Nahm transform recovers the spectral 
data of magnetic monopoles, as described in section 12.2.21 The n rational functions are 
simply 

'^(Va) = , w^^'^^R , = ^A+MVA - Ba)-'ua . (5.46) 
det {rjA - Ba) 

These rational functions involve all independent local coordinates, thus invoking the same 
argument as in section 12.2.21 for large separations between the constituent monopoles, we 
obtain a parametrization of in terms of nk 3-dimensional locations and nk phases, de- 
scribing all together nk monopoles. 

Although above we have obtained only a local parametrization of the twistor space, 
the spectral data that emerged is sufficient to calculate the exact hyperkahler metric on 
This is done using the generalized Legendre transform |l79l and is work in progress, close 
to completion. 
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Lattice gauge theory provides a non-perturbative framework to investigate the relevance 
of calorons in a dynamical context. This is necessary, because the mere existence of our 
new solutions does not guarantee that they play any role dynamically Nevertheless, the 
non-trivial constituent nature of the solutions is intimately tied to the average Polyakov 
loop, which is set dynamically as discussed in section [L3l Also a semi-classical calculation, 
based on the 1-loop determinant in a caloron background, provides evidence in favour of 
their dynamical relevance |19|. 

There are several lattice methods to investigate the low-energy behaviour of a gauge 
ensemble. Two main methods are cooling and the study of the zero-modes of the improved 
chiral Dirac operator. The first, cooling, only deals with the bosonic sector and probes the 
long range correlation in the gauge fields, their topological content in particular I92ll93li94t . 
Our exact results on the caloron gauge field may be compared with the lattice findings. 

The second method probes the fermionic sector and uses the same philosophy as we 
have for motivating our study of the fermion zero-modes. That is, the low lying spec- 
trum of the Dirac operator carries information on the topological content of the underly- 
ing gauge field and is responsible for several non-perturbative aspects of QCD. We have 
shown how continously changing the boundary condition in the compact time direction 
for a zero-mode makes it hop between constituent monopoles. The surprising fact is that 
on the background of Monte Carlo generated configurations - which are usually very 
rough - the zero-modes of the improved chiral lattice Dirac operator show very similar 
behaviour. This property can be used to probe the monopole content of the gauge fields 
iHITTOl. 

The results of our exploratory lattice investigations - which nicely show how numer- 
ical simulations complement the study of the formal structures - is summarized in this 
chapter. The computations were done with modest computer resources and are sufficiently 
accurate for our purposes. For dedicated large scale Monte Carlo simulations we refer to 
the works mentioned above. 



91 



6.1. Lattice gauge theory 



6.1 Lattice gauge theory 

In order to simulate SU{n) Yang-Mills theory on a computer first and foremost spacetime 
and the gauge fields must be discretized. This is usually done by replacing spacetime by 
a 4-dimensional square lattice with 4 group valued fields defined on each oriented link. 
These Lr^(x) link variables represent parallel transport - or path ordered exponential - 
from X to X + fl where x is a site and x + fl is the next site in the }i direction. The distance 
between x and x + fiis the lattice spaceing a. The gauge potential of the continuum theory 
is encoded in the link variables. In order to make the continuum limit more explicit one 
can define a lattice gauge potential by 

Lr^(x) = exp(flgYMA^(x)) . (6.1) 

Gauge transformations g{x) are defined at each site and act on the link variables by 

Uj,ix) ^ gix + fl)U^{x)gix), (6.2) 

the form one would expect from a path ordered exponential from x to x + p of a continuum 
gauge field Aj,(x). The assignment that the inverse is associated to a link with opposite 
orientation is also quite natural. 

From the transformation property above it is clear that the trace of any ordered product 
of link variables along a closed loop is gauge invariant in exactly the same way as a Wilson 
or Polyakov loop is gauge invariant in the continuum. In particular the trace of a basic 
plaquette of sides a, 

U^y{x) = U~'^{x)U~'^{x + v)Uy{x + ]l)U^{x) , (6.3) 
is also gauge invariant. The action proposed by Wilson IIIOIII is then defined to be 

S = i^E (l - ^ReTrLr^^.(x)^ , (6.4) 

where the sum is taken over every possible basic plaquette in the lattice and /3 is the bare 
coupling, not to be confused with 1 /fcg T as it was used in previous chapters. One can show 
that in the classical continuum limit, — > 0, the Wilson action agrees with the continuum 
action corresponding to (|l.l|l if one sets /3 = In/gy-^. Here gyu is the bare coupling, 
which is subject to renormalization in the quantum theory. The above matching is purely 
a classical one. 

Since we are dealing with finitely many degrees of freedom the Feynman integral is an 
ordinary - yet multiple - integral over SU{n). The number of integration variables grows 
with the volume and quickly becomes too large for explicit evaluation. Instead, a Monte 
Carlo method is employed to compute observables. 
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Simulations at zero temperature are done in a 4-dimensional box with the same number 
of lattice points on each side, provided the lattice spacing is the same in every direction. 
Non-zero temperature is implemented by choosing less number of points in the time direc- 
tion than the spatial directions (assuming a homogeneous coupling in all 4 dimensions). 
Once the number of lattice sites is fixed, the lattice spacing and hence the physical tem- 
perature is controlled by /3 = 2n/gYM- Increasing /3, i.e. lowering gyu, reduces the lattice 
spacing as dictated by the running of the coupling, 

at 1-loop order. 

We have already met interesting lattice observables, the Wilson or Polyakov loops, but 
one can easily construct many more. For our purposes an important one, apart from the 
Lagrangian, is the topological charge density. A widely used lattice version is 

1 

q{x) = -^g-2£fivpaT^rUjiy{x)Upaix) . (6.6) 

6.2 Phase transition 

The order parameter of the confinement - deconfinement phase transition is (p(x)), the 
vacuum expectation value of the trace of the Polyakov loop ()1.25|l . For SU{2) it behaves as 

and the phase transition is second order. For high temperatures j6 > /3c, in the deconfined 
phase, the center Z2 symmetry which interchanges the vacua p = 1 and p = —1 is spon- 
taneously broken. In the confined phase, j6 < /3c, this symmetry is restored; see section 

We will demonstrate the phase transition ona4xl6xl6xl6 lattice with periodic 
boundary conditions using a simple heatbath Monte Carlo algorithm [1021. The appropri- 
ate quantity to be measured is the average (p) of {p{x)) over the spatial lattice. This is 
plotted in figure |6H against /3. We see that the phase transition occurs around the critical 
point /3c ~ 2.3. 



6.3 Cooling 

In the process of cooling, the link variables are altered such that the total action is lowered. 
Updating every link variable over the lattice is called a sweep. Sweeping through the lat- 
tice many times, one is necessarily ending up with self-dual solutions as these are minima 
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Figure 6.1: Phase diagram of SU{2) gauge theory. The vacuum expectation value (p), half the trace of the 
average Polyakov loop, is plotted as a function of /5. 



of the action. However, due to the finite discretization, instantons can "fall through" the 
lattice if their size becomes smaller then the lattice spacing. Following our remark on the 
classical continuum limit, a ^ 0, this effect can only be 0{a^). Concretely, the discretized 
1-instanton solution of size p has lattice action 1911 



S(instanton) = 2n^^ 1-^ - - +••• (6-8) 




for gauge group SU{2), which will be assumed throughout. Indeed, the leading term is 
the continuum 1-instanton action STT^/gyM correction is such that a decreasing size 

decreases the action. Thus in the process of cooling, where the action is always lowered, 
instantons will tend to shrink. Their size will eventually reach the lattice spacing and 
assuming a location in between lattice sites, they disappear. 

This is illustrated in figure \62\ for a 16 x 16 x 16 x 16 lattice with periodic boundary 
conditions. An initially random - or hot - configuration was cooled to reach a configura- 
tion with roughly 6 instantons, after which we have plotted both the action (measured in 
units of S/r^/gYM) ^^'^ topological charge against the number of cooling sweeps. Each 
sharp drop by one unit corresponds to one instanton falling through the lattice. Between 
the falls the topological stability is reflected by long plateaus, where the action and the 
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Figure 6.2: Instantons falling through the lattice one by one. The action in units of Stt^/^ym ^^'^ 
toplogical charge is plotted against the number of cooling sweeps. The action is always slightly higher then 
the topological charge. 



topological charge roughly agree, indicating a self-dual configuration. At the final stage 
the last remaining instanton disappears as well, leaving the trivial vacuum behind with 
each plaquette being the identity of SU{2). 

If one is to investigate stable self-dual solutions on the lattice, the above shrinking phe- 
nomenon must be eliminated. The following action serves this purpose. 



48 



El 



:Tr 



(6.9) 



where the first term represents the Wilson action (16.41) , while the second means a sum over 
every 2x2 plaquette. Note that for SU{2) all traces are real. Clearly, Si is the Wilson 
action. The usefulness of Se is that the discretized 1-instanton will have an action 



(instanton) =2n^^\l 



5 \p 



+ 



1 2 - 15£ 
210 



+ ... 



(6.10) 



which shows that once a negative value is picked for e, an instanton will be stable against 
shrinking I9H . The value e = minimizes the lattice artefacts to this order - thus Sq is 
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called the improved action - and also leads to stable instantons, but the stabilizing force is 
greater with a negative value of e. The choice e = —1 is called over-improvement. 

Over-improved cooling stabilizes instantons so well, that starting from a random initial 
configuration with very high action, one can easily get stuck at multi-instanton configu- 
rations with high topological charge. Thus in order to investigate self-dual configurations 
with charge 1, 2 or 3, it is most useful to first apply Wilson cooling and only switch to 
£ = — 1 when the action is already below 2, 3 or 4, in units of STT^/gyM- 



Figure 6.3: An SU{2) caloron of unit charge on the lattice, obtained by over-improved cooling. The topo- 
logical charge density is on the left, half the trace of the Polyakov loop on the right, showing that the lump 
to the left is a monopole of charge -1 and the lump to the right is a monopole of charge +1. 

At finite temperature the situation is the same as at zero temperature. Over-improved 
cooling, however, will have an additional advantageous feature, it will tend to separate the 
constituent monopoles of opposite charge. Thus not only the calorons will be stable against 
falling through the lattice, but also their constituent nature will be revealed. We have 
performed simulations on a4x 16x16x16 lattice, with periodic boundary conditions in 
all 4 directions. 

In figure 16.31 a charge 1 caloron is shown. The action density is summed in one di- 
rection and plotted in the remaining plain at a fixed time slice. It was obtained by first 
Wilson cooling a random initial configuration for 665 sweeps until the action decreased 
below 1.8 X Szr^/gYM- Th^ri it was followed by 1200 over-improved cooling sweeps and 
the improved action Sq became 1.014 in units of ^^^/gyw whereas the topological charge 
0.877. The reason for the topological charge only being an integer up to 12% is that we used 
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Figure 6.4: An SU{2) caloron of charge 2 with well-separated constituents on the lattice, obtained by over- 
improved cooling. The topological charge density is on the left, half the trace of the Polyakov loop on the 
right, which shows that the two smaller lumps are monopoles of charge -1 and the two higher lumps are 
monopoles of charge +1. 

the simple operator 1)6. 6|l which has O(fl^) lattice artefacts, whereas the improved action So 
was deliberately chosen such, that only O(fl^) artefacts are present. 

Higher charge calorons - discussed in great detail in the previous chapters - can be 
illustrated as well. Figure 16.41 shows a charge 2 caloron. The topological charge of this 
configuration is 1.82, and has improved action 2.001 in units of STT^/gyM- obtained 
by first Wilson cooling for 580 sweeps until the action dropped below 2.8 x Sn^/gyf^ and 
then switching to over-improved cooling for 2150 sweeps. 

This concludes the presentation of our numerical studies, complementing the exact 
results of the preceding chapters. For more extensive investigations, see I95ll96ll98ll . 
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We have discussed the explicit construction of multiply charged caloron solutions with 
non- trivial holonomy. For gauge group SU{n) and topological charge k, the non- trivial 
Polyakov loop causes the caloron to dissociate into nk massive magnetic monopoles. Thus 
a caloron of charge k should not be seen as an approximate superposition of k charge 1 
calorons, but rather as the approximate superposition of nk constituent monopoles, each 
carrying fractional topological charge. The constituents have their own identity, they do 
not know to which charge 1 caloron they belong, the non-linear superposition of all of 
them gives rise to the full 4 dimensional gauge configuration. The action density clearly 
shows nk lumps, once the constituents are separated far enough from each other. The nk 
monopoles come in n distinct types, each being charged under a different U{1) subgroup. 

Along with the bosonic sector we have investigated the fermionic sector as well, and 
have found the zero-modes of the Dirac operator in the caloron background. The compact 
time direction corresponding to finite temperature necessitates a choice of boundary con- 
dition for the fermions and the index theorem dictates the existence of k zero-modes for 
every such choice. In order to describe the correct partition function, physical fermions 
are required to be anti-periodic, however, for diagnostic purposes one may impose for the 
zero-modes periodicity up to an arbitrary phase. The usefulness of employing zero-modes 
to identify the constituent monopoles lies in the fact that by fixing this phase, the zero- 
modes localize to only one type of monopole and are blind to the remaining n — 1 types. 
By continously changing this phase one can detect all the nk constituent monopoles, where 
the k zero-modes jump from one type to the other whenever the phase passes through one 
of the eigenvalues of the asymptotic Polyakov loop. In this way one can analyse each 
monopole separately, effectively deconstructing the caloron into its bare constituents. 

An important tool we have employed to clearly see to what extent the constituents can 
be described as point-like objects, is the abelian limit. In this limit all exponential contri- 
butions originating from the non-abelian cores are neglected, only algebraic tails survive, 
giving rise to abelian fields. The masses of the monopoles in this limit are effectively infi- 
nite, and the zero-modes localize to the non-abelian cores. These collapse to zero size, thus 
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the zero-mode density exhibits a singular behaviour with singularity structure tracing the 
cores. The remnants of the cores in the abelian limit are in general extended as seen by the 
zero-modes, however, we have established that once the separation between constituents 
is large, they become point-like. 

The moduli space of calorons factorizes into nk copies of x R"^ for large separation, 
confirming the right description in terms of nk phases and nk 3-dimensional locations. 
This is in contrast with instantons on M.^ where the correct description is in terms of k scale 
parameters, k 4-dimensional locations and relative gauge orientations. As to the algebraic 
geometry of the caloron moduli space we have established a correspondence with stable 
holomorphic bundles on the projective plane that are trivial on two complex lines. This 
was derived from a similar correspondence between instantons on and bundles that are 
trivial on one complex line, which ultimately is a result of the compactification of ~ 
to CP^. This implies that the appropriate way of approaching the problem is not through 
the factorization R^/Z ^ x R'^, but rather through the embeddings S^xR^^C*xCc 

C CP^, at least from a holomorphic point of view. 

This observation may turn out to be relevant for the study of doubly periodic instantons 
as well L46,, 47J. They are defined on x R^ and traditionally - especially in the Nahm 
context - are studied through the factorization R'^/Z^ ~ x R^. Now since x R^ ~ 
C* X C* which again can be compactified to CP^ by the sequence of embeddings C* x C* C 
C* X C C C CP^, we suspect the existence of a correspondence between the moduli 
space of doubly periodic instantons and stable holomorphic bundles over CP2 with some 
additional triviality constraints. 

We believe that our results on the moduli space of calorons may be useful for com- 
pactified supersymmetric gauge theories as well IITSI . Similarly to the formulation on R^, 
a number of correlation functions are saturated by instantons - or calorons in the com- 
pactified case. The computation of these correlation functions involve integration over 
the moduli space of instantons, which in fact can be performed by localization techniques 
| |103( . Understanding the metric on the caloron moduli space is a necessary ingredient to 
compute these correlation functions for the compactified theories. If one wishes to com- 
pute the contribution from every topological sector, a summation over charge is necessary, 
thus the metric is needed for the moduli space of calorons with arbitrary charge. This is 
where we believe our results on the twistor description of the moduli space may turn out 
to be useful. 

We have saved our primary motivation to the end. Understanding permanent quark 
confinement in QCD remains to be a challenging task. What we have set out to investigate 
was the nature of the fundamental topological excitations in the confined phase. We have 
found that instantons are composites and are made up of monopoles due to the non-trivial 
Polyakov loop background that is present in the confined phase. The instanton liquid 
model 111041 has been very successful in demonstrating chiral symmetry breaking and we 
believe that a similar model, but now with the constituent monopoles taking the place of 
instantons, holds considerable promise for the future. Especially, because the presence of 
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magnetic monopoles in the QCD vacuum is in fact required for the dual superconductor 
picture to emerge [20, 21J. In the 3 dimensional Yang-Mills-Higgs model monopoles have 
been shown to be responsible for confinement long ago [105J. 

It is well-known that the dilute gas assumption of the instanton liquid model breaks 
down for large instantons and also that small instantons are not sufficient for explaining 
confinement. In our constituent monopole picture the scale parameter of large instantons 
is converted to the separation between monopoles. It is not unlikely that the density of 
constituents at low temperature is so high that they form a coherent background and as 
such will no longer be recognized as separate lumps. With high quark density leading to 
deconfinement, it may perhaps be that a high constituent monopole density will result in 
confinement. 

Whether or not constituent monopoles play a role in a dynamical context, remains a 
difficult non-perturbative question. The natural testing ground is the lattice, which shows 
growing evidence in favour. Simulations were performed very much along the lines of 
our investigations, first identifying the semi-classical gauge structure |93|, followed by 
studying the localization properties of chiral fermion zero-modes 11731 174 1. These studies 
- some of which we have also touched in an exploratory fashion - seem to imply that 
constituent monopoles are the relevant degrees of freedom in the confined phase. 

Could it be that the non-perturbative dynamics of constituent monopoles is such, that 
it causes the quarks to confine? At present, we do not have an answer to this question. 

It should be noted that the idea of instantons having more fundamental building blocks 
is by no means new. The moduli space of fc-instantons in the 2-dimensional CP"^ model 
can be parametrized by nk complex numbers (the remaining moduli are not relevant for 
what follows). All of these have an interpretation as 2-dimensional locations for so-called 
instanton quarks |106|, or in our terminology, constituents. The 1-loop determinant in the 
instanton background defines a statistical ensemble for these nk instanton quarks which 
may be mapped to a 2-dimensional Coulomb gas 1 107, .108 1. Since the topological charge 
is k and there are nk constituents, the topological charge carried by each Coulomb particle 
is fractional. 

Despite the fact that the 2-dimensional CP"^ model shares some features with QCD, 
such as confinement, a similar picture is difficult to arrive at in gauge theories. Never- 
theless, the results presented in this thesis further suggest that the description in terms of 
instanton quarks is indeed the relevant one. It has been argued recently that in the con- 
fined phase of QCD the effective low energy theory is of a Sine-Gordon type |1(H1 11101 . 
which also has a Coulomb gas representation lll 1 1 Ml 121 . Just as for the CP"^ model, but 
now in four dimensions, the Coulomb particles carry fractional topological charge 11131 . 

We are hopeful that our results and related work mentioned above will ultimately tie 
together into a consistent picture of a confining QCD vacuum. 
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We present below the polynomial whose root is 'Y for charge 2. It follows from our dis- 
cussion that this funcion is harmonic almost everywhere, see section l4.1.2l for details. 
We have 

flo + flil^^ + a^f^ + ^3^^ = , 

with 



4x2x2k6 - 4k2x4x2 + exjk^ - 2x|x2i^4 _ 21,6 + ^4 ^ _ 5^41,2 ^ 2x2^2 _ 6k4x2+ 

4k6x2 + 4x6k2 + 4x2k4 - 6x4k4 + xfk^ + xfk4 - 2k6xf + 4x^k4 + 2x2k8 - 6x|k6- 
6k6x| + 2xf x2 + 2x2x6k2 + 2xfk6x| + 6x^xfk4 + ex^xfk^ + x^xfk* + 2x§k4x2- 
6^2^2-^ 4k x-j^x^ —f" ^x^k X-J^ ~|~ 8k x^x^ ~1~ 6k x^x^ 2x2x2k 6x2k x^ 6x2X'^k 
lOx^x^xfk^ + 6x^x^xfk4 + 6x|x2xfk4 - 2x^x2 xfk^ + 2x^x^xfk4 - xf - 2x6 ^ ^8_ 
6x^xf k2 + 2x6x2k2 - 2xf k^x^ + lOk^x^xf + 6x|xf k^ + 2xf x^ + xf x| - 6x^x4 - 2x^k2+ 

2x2k4 + 2x4k2 + 2k6x2 + 2k2x6 + 2k4x§ - 2x2k8 - Ix^^^ - ex^k^ + X^k^ + k4x§ + 

2xjk6 - 2k2x8 + 2x2x^ - 6xf x^ + 4x§x2 + 4x^x^ + 6xf x^ + 2x^x4 + 2x^x§ + x\x\-^ 

Ay2y2y2-L.6 0^2^4^21,2 /:i,2^2^2^2 _ /:^2^2^2i,4 9^6 _i_ ^8 1 r4 Ry'^y^V^ 

2x§x^k2 - 6xf x|k2 + lOxf x|k2 + 4x§xfk2 + 2xf x^k^ - ek^x^xf + lOx^x^k^- 

2 2i 2 1 4 2 2 4 4i 2 242 422 222 42 2\^ 

6xo — 6x2x2k + 6x2 
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/:y2y2i,6 1 4^6^21,4 _ £,^4^21,6 _ AylyZ^^e _ AxJlyAyl 1 £,y4i,4 _ jyAyl^A^ _ 91^61 

+ k« - 6xfk2 + 2x2k2 - 6k^xl + 4k6x2 + 4;cfk2 + 4x^k4 + 6x|k4 + x|kS- 
2k X2 H~ 4x2k ~|~ 2x2k 6x2k 6k X2 ~\~ ~\~ 2x2X'2k ~|~ 2x2x2k ~|~ 6x2x2k ~|~ 

f.y.ly.l^.6 I Y^Y^v4: I 9^6^21,4 _ (:y.2y.4]^4 _ Av4y2y4 , Oy^y^VS _i_ qi<-4y2y2 , (£11,4^4^4 _i_ 

2x2X^k 6x2X^k 6x2X^k 10x2x2x2k ~\~ 2X2X2 ~|~ 2X2X2 H~ 6x2x2x2k ~|~ 

f.yAy2y2-\^4 _ 9y2y2y41<-2 i 0^2^2^41,4 i ^4 9y6 _l y^ _l 8y4y21<-2 i i:,Y2Y2L-2 _i_ 9y6y2l.2_ 
\j.A-2-^2 1 z../\.2-^2 2 z.A.2-^2 2 1 ■^-^2 2 l'-a-2-^2 ^ja.2-^2'^ ^^2 "^2 

2x2x2k 10k X2X2 ~f~ 6x2x2k 2X2X2 ~[~ X2X2 6X2X2 2x2k ~|~ 2x2k ~|~ 

2xjk^ + 2k^xl + 2k^xl + 2k4x6 _ 2xlk^ - Ix^k^ - Sxjk'^ + xjk^ + k^x^ + 2xjk^- 

9V,2--8 I 9-y.2-^2 f.^4:^2 1 a ^6^2 1 -^-6-^.2 1 /:-y.4--4 _i_ -^.4-^4 1 -^2-^2-^2i,6 9-^.2 --4 -.2i,2 

x_ IN. A.2 z.A.2-A'2 vjA.2-^2 ".A-2-^2 ".A-2-^2 ^j.a-2-^2 -A-2-^2 "-A-2-^2 3 z..a-2'^2 3 

/:i,2 ^2 -.2 -.2 _ /:^-.2^2^2i,4 _ ^^6 \ -L — Ir'^r^V'^ — Ar^r^lr^ 1 1 r)-.4-.2i,2_ 
U-N. -1-2-^2 2 0^2-^3-^2 3 3 3 ^A.2-^2 uA.2-^^2^ ' ±ua.2-^^2^ 

4x2x2k 2x2x2k " 6k X2X2 10x2x2k " 6x2x2k ^ 8k X2X2 6x2x2k ~\~ 

1 9 ^4 ^2 ^2 1 J_ -.2 ^2 ^2 I (0^2^4:].4 _ ^ y.1 y.l\.6 _ A].! y.^ yl 1 ^4^2l,4 1 

UA2-^2 2 ^-^2-^3'^2 ^^.2^2-^2 ^ ua.2-^^2'^3 / ^*-?a.2-^2 z-X'2^X'^r< rtis. ^2-^2 ^ a2-^2 ^ 

9Y4k2 _ Y2k2 -L k4r2 -L Y2k4 J- 9r4k4 -U r6k4 _ r4k6 -V(>y2 , yAy2 , ^4:^2^4 \.4y.2y.2_ 
^A.'^rs. a.'^js. is. ^2 ~r ^2^^ \ ^■\-'2^js. I A-f2^ry A.'2^rs. js. ^2 ^2-^2 ~'~ A2A2 ^2-^2 

^4 _|_ ^6 _|_ ^2^2],2 I 9-.4-.2i,2 i i,4-.2-.2 i 9^4^2],2 i 0-.2-.4 i y2\.l _ 9r2v4 i ^4i,2_|_ 
^2 A2 I A2A2 ' ^A.'^A.'2^]\. I is. A2A2 ^ ^^2-^2 ' »_'A2A2 i A2-*n- z.a2-»n- i -A,2-r^ \^ 

1,6^2 _ 9V2-.6 I V4-.6 i y4L'4 _ r^L-^ _ 9r2r^ _u 'Xr^r'^ -U r^r^ _ 9Lr2-y.2--2^2i 
JS. ^2 Z-JN. ^2 -^2 ^ 3 3 ^a.-^a.'^ \ ^^^2-^2 \ ^2-^2 ■^js. ^2-^2-^2^^ 

9-^2-^2-^21,4 I v6 _ v4 _ Av4v2v2 _ 9-^4-^21^2 i riylylA^l _ jir2y.2y.A i Ov4v2v4_ 
ii.A-2'^3'^2 3 3 TCiA'2'^2 ^.A^2*'*'3-'^ \^ ^a^^a^-^js. ^js. ^'^^'2. *-'t^2*'*'3-^ 

—^2 —^21^2 I ^il^4 —^2 —^2 I —^2 —^2 —^2^ 
2 3 '-'J^ 3 ~'~ ^'^2 2 3/ 
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ai = -48xf x^k^ + 16xf x^k^ - 144x|x§k4 - 96x|xf k* - 48x|x|k4 + 48x§xjk6+ 
192x2'^^-k 96k x^^^ 16k ^2^3 48x2x2k ~|~ 48k •^2*^3 96x2x2k 
48x2x2k H" 48x2x2k H" 48k ■^2*^3 16x2x2k 144x2x2k H" 16x2x2k 
48x2x2k 16x^X2 16X2X2 48x2x2k ~|~ 96x2x2k ~|~ 32x2x2k 96k dC'yX''2^ 
48x2x2k ~|~ 96k 48x2x2k ~|~ 192x2x2k 48x2x2k 16x2x2k 

Uk^^xlxl - 48ki°x|x2 - ASx^x^k"^ - 48k8x2x4 + 96x^xjk^ - ASx\x\k^- 
16x2x2k 48x2x2k 16x2x2k ~|~ 96x2x2k ~|~ 32k X'2^X-^ ~\~ 32x2x2k ~|~ 
96xfx2k« - 48x^xfk^ - ASxjxp^'^^ - 16xf xfk^ + gex^x^k^ + 48ki°x2x2- 
96x4xfk4 - 48ki0x^xf - 144x§xf k^ + 144x§xf k^ - 96xfxpi^ + 144kSxf x^- 
96k8xf x^ - 48x6x2k4 + lek^^x^x^ + 48xf x^k^ + 48x2xfk4 - 16xlxlk^°+ 
48xf x^k^ + 48k4x2x| + 48x6x2k*^ - 48x8x2k*^ - UAx^xlk^ + 16x2x^k4+ 
48ki0x2x6 + 48kSx2xf - ASxfxpi^ + lex^x^k^ - lek^x^xf + 48x8x2kS- 
48k8x^x2 - 48xf x§k4 + 48x^x^k2 - 96k4x|x| - 48x^x^k4 - 96x^x|k8+ 
48xf x|k6 - ASx^xlk^ + 48x|x6k6 - 48x^x§k8 - 96x|x|kS + 144x|x|k*^+ 

48x2x2k ~\~ 48x2x2k ~|~ 48x2x2k ~l~ 192x2x2k ~\~ 144x2x2k 64X2X2X2 

64x^ 

^3'^! 32x2x2x2k 96x2x2x2k 96x2x2x2k ~|~ 
240x2x2x2k 48x2x2x2k ~\~ 144x2x2x2k ~\~ 176x2x2x2k 96x2x2x2k 

1 1 9-^2-^6-^21,6 _ ^^-^2-^2-^61,10 _ q^t^4-^4-^2i,6 i OAOv^v^v^V^ -U QAv4v2v4i'4_ 
xxz.A2'^2 2 wriA2'^2 3 7L/A2'^2 3^ ^^^-''^2 x 3 7L/A2'^2'^2 

112x2x2x2k ~|~ 32x2x2x2k 64x2x2x2k 48x2x2x2k ~|~ 176x2x2x2k 

ARr^Y^Y^k^ _i_ M^Y^r^r^k^ 4- ^QJY^r^Y^k^ -\- Af^Y^Y^Y^k'^ — QftY^Y^Y^k^A- 

Tct^./V2 1^2 2 xTcTc./V2 ''*'2 2 y^A^'2^A^'Y A^^js. \^ Tct^./V2*'*'2 ]^ ^w^2 ''*'2 2 

ARr'^Y^Y^k^ — ^AAY^Y'^Y^k^ -\- A9,k^^ y^ y^ y^ -\- ARY^Y^Y^k^ — ^AAk^Y^Y^Y^ — 
^tj,A-2-^2 ^'^'2_ \ A.2-^2-'^2 ^c/.A-2-^2 3 ± -t jv A.2-^2 2 

96r2r4y4l<.6 i Q/:,y4^4^2i,6 _ 40^4^4^21,8 _ 0^:^4^2^21,6 1 4.oi,8y4^2^2_ 
17^/^2-^2 3 ^ yiJA2-^2 3 rroA2A2A2J^ i7ua2-^2 3 1^ Ttojs. a2-^2'^3 

48ki0x? X2 ^2 

+ 96k^x? X2X2 ~l~ 48x2x2x2k 48x2x2x2k ~|~ 96x2x2x2k 

flo = 64x^x§x5k4(k - 1)2 (k + 1)2 . 
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